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Chapter 1 

Introduction. 



In this thesis we study sets of points in the plane and their Voronoi diagrams, in 
particular when the points coincide. We bring together two ways of studying point 
sets that have received a lot of attention in recent years: Voronoi diagrams and 
compactifications of configuration spaces. We study moving and colliding points 
and this enables us to introduce 'limit Voronoi diagrams'. We define several com- 
pactifications by considering geometric properties of pairs and triples of points. In 
this way we are able to define a smooth, real version of the Fulton-MacPherson 
compactification. We show how to define Voronoi diagrams on elements of these 
compactifications and describe the connection with the limit Voronoi diagrams. 

Voronoi diagrams and supermarkets. 

Consider all supermarkets in a city. We divide the city in sectors by considering 
the closest supermarket: all people in the sector of some supermarket are closer 
to this supermarket than to any other supermarket. The Voronoi diagram of the 
supermarkets is this subdivision of the city into sectors. The Voronoi cell of one 
supermarket is just the sector of the supermarket. An example of a Voronoi diagram 
of six supermarkets, numbered 1 to 6, is given in Figure 1.1. 




Figure 1.1: The Voronoi diagram of six supermarkets, labeled 1, 2, 3, 4, 5, 6. 

In the figure we see some boundaries, for example the boundary between super- 
market 1 and supermarket 3. People living on this boundary are at equal distance 
between those supermarkets. The boundary is part of a line, the dotted line in the 
picture, which is called the bisector of 1 and 3, because on one side of the line people 
are closer to supermarket 1 and on the other side people are closer to supermarket 
3. So, the bisector divides the plane into two sectors or half-planes. The sector 
containing supermarket 1 is called Voronoi half-plane vh(l, 3) and the sector on the 
other side of the bisector is Voronoi half-plane vh(3, 1). We can express the sector 
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V(l) of supermarket 1 in terms of these Voronoi half-planes: 

V(l) = vh(l,2)nvh(l,3)nvh(l,4)nvh(l,5)nvh(l,6) (1.1) 

Of course this formula just states that you arc in the sector of supermarket 1 if you 
are closer to supermarket 1 than to 2, closer to supermarket 1 than to 3, etcetera. 
If we want to code the bisector of 1 and 3, we need only two ingredients: one point 
on the bisector and the angle that the bisector makes with the horizontal line. This 
will be important later on. 

Driving supermarkets: changing Voronoi diagrams. 

In the Dutch countryside food is supplied by little supermarket lorries, or supcrvans. 
These vans drive through the countryside looking for customers. In the part of the 
countryside we are considering there are five supervans. If we know their positions 
at a certain moment, we can determine the Voronoi diagram of the five vans at 
that moment. But, if the vans drive continuously the Voronoi diagram of the vans 
changes continuously as well. In Chapter ^ we model these driving supervans and 
their changing Voronoi diagrams. For every supervan there is a curve that gives at 
any time t the position of the supervan. If we want to know the Voronoi diagram at 
time t, we specify t in the curves describing the position of the vans. Using positions 
given by the curves we compute the Voronoi diagram. 




a b c 

Figure 1.2: a. The positions of five driving supervans before collision, b. The Voronoi 
diagram of the vans, just before collision, and c. at collision. 



There is only one problem. At one day a very unfortunate accident happens: all 
vans run into each other at t = 0. How should we define the Voronoi diagram of the 
vans at t = 0? In this particular example we decide to analyze the situation just 
before the accident happens. That is, we let time run backwards and investigate 
what is going on for small negative t. The route covered by the five vans before the 
accident is shown in Figure 1.2. a. A Voronoi diagram of the positions of the vans 



just before the accident is given in Figure 1.2.b. 

Our strategy is to define a Voronoi diagram at t = that is consistent with the 
situation just before t = 0. We call such a diagram a 'limit Voronoi diagram' and 
an example of this that resembles the Voronoi diagram in Figure 1.2.b is shown in 
Figure 1.2 .c. Compare the two diagrams. Note that the directions of the outgoing 
edges in the two diagrams are very similar. The Voronoi cell of van 3, in the middle 
of the diagram of Figure Figure 1.2.b seems to have completely disappeared in 
Figure |l.2j .c. These two remarks address questions that we answer in Chapter ^ 
and the chapters after that: what information is needed to to create a diagram as 
shown in figure 1.2, c? And, can we still say something about the Voronoi cells that 
seem to disappear in the limit, like the cell of van 3? 
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Voronoi diagrams for coinciding points. 



1 . 




2 • 



a b c 

Figure 1.3: The Voronoi diagram of two coinciding points. 



We restrict ourselves for a moment to the case of two points, p\ andj>2- The position 
of each point at time t is given by a curve. Suppose that the two points meet at 



t = 0. Such situation is shown in Figure 1.3.a-c. In Figures 1.3. a and b, the two 
points arc still distinct. Therefore we can draw the line l\2 that passes through both 
points. This line makes some angle a±2 with the horizontal axis. So this angle «i2 
is in some sense the angle of the points p\ and p2- As the positions of the points 
depend on time t, the angle ai2 also depends on time, that is ai2 — a\2{t). The 
Voronoi diagram of points pi and p2 is determined by the line at equal distance 
from the points p\ and p2- That line is exactly the bisector of p\ and p2'. the line 
perpendicular to line 1x2 passing through the middle of the line segment 12. If p\ 
and p2 coincide, we define the middle of the line segment 12 as the point p\ = p2 
itself. And we define the angle 0:12(0) as the limit for small negative t of ai2(i). 
Now we define the bisector at t = in terms of this angle and this middle point. 
That is, the bisector of pi = p 2 is the line passing through pi — p 2 perpendicular 
to the direction 0:12(0). But this implies that we have created a Voronoi diagram 
for the two coinciding points p\ and P2 ! 

We drop this particular example but conclude the following: we can define a limit 
Voronoi diagram for two coinciding points p\ and P2 if we know the following infor- 
mation: 

- the position pi = p 2 of the coinciding points. 

- an angle ot\2 mod 2tt. 

In this way we can define a limit Voronoi diagram for an arbitrary number of co- 
inciding points as well: as long as we have for every pair of coinciding points pi 
and pj a position pi = pj and an angle onj , we can define the bisector of pi and pj . 
And, using the bisector we can determine the two Voronoi half-planes vh(jpi,pj) and 
vh(pj,pi). Now we are done, as any Voronoi cell can be expressed as an intersection 
of half-planes, as we saw in Equation [l.l| . 

We work along these lines in Chapter ^. For example, the curv es d escribing the 



positions of the two points p\ = pi(t) and P2 — p%{t) in Figure F3 are given by 
Pi{t) = (t 2 ,t + t 2 - M 3 ) and p 2 (t) = (t, -t - 3< 2 + 2t 3 ). We only allow curves given 
by pairs of polynomials in t. Wc call points described by such curves polynomial 
sites. After having defined Voronoi diagrams for polynomial sites we show how to 
determine the Voronoi diagrams without having to compute all bisectors. 

First we assume that we have a set of n polynomial sites that all coincide at t = 0. 



So, this is the situation in Figure 1.2. a, where five polynomial sites cluster together. 



Leaving out polynomial site p% (t) in the beginning will not change the shape of the 



limit diagram as we see in figure 1.2.c. We show in Section 4.5 which sites exactly 



8 



Chapter 1. Introduction. 



can be left out without changing the boundary of the diagram. So, we want to 
characterize those sites that are somehow in the interior of the cluster. Of course 
we could as well characterize the sites that are on the exterior of the cluster. This is 
done in terms of the convex hull of the sites in the cluster for small enough positive t. 




Figure 1.4: Plugging Voronoi diagrams. 



For a general set S(t) of n polynomial sites, the boundary of the Voronoi diagram at 
t = is determined by plugging. We demonstrate plugging in Figure 1.4 where we 
compute the Voronoi diagram of S(t) — {p\, . . . ,p%\, for pe = (—3, 2), p~j = (0, —3) 
and ps — (3,2). The other polynomial sites form a cluster in th e o rigin. In fact 
we used this cluster already for the supervan example in Figure F2. The cluster 
locations of S(t) at t = are the distinct positions of the sites at i = 0. So, we 
have fo ur c luster locations. The Voronoi diagram of the cluster locations is given in 
Figure IT. a. Next, the cell of every cluster location is filled, if necessary. Three cells 
in our example correspond to a single point, so no filling is needed there. The limit 
Voronoi diagram of the points clustering in the origin, is shown in Figure 1.4.b. We 
simply plug this diagram in the appropriate cell in order to obtain the limit Voronoi 
diagram of S(t). The result is shown in Figure 1.4 .c. 



The limit Voronoi diagrams are really new diagrams: in gen eral they can not be 
realized as classic Voronoi diagrams, for example Figure F4.c. In classic Voronoi 
diagrams, every Voronoi cell has positive area, for example. This is not true for 
limit Voronoi diagrams. 



Compactifications of configuration spaces. 

A collection of n points in M 2 is often called a configuration. The configuration space 
of n distinct points in R 2 is just the space that consists of all possible configurations 
of n distinct, ordered points. Suppose for example that n = 3. Any element of the 
configuration space consists of three lebeled points Pi,j»2,P3 G K 2 such that: p\ 
is distinct from p2\ p\ is distinct from p^; and p2 is distinct from p%. A natural 
way of describing the configuration {pi,P2,P3) is by listing the coordinates of pi, 
P2, and pz- But if we list all six coordinates in single file, we obtain an element 
c = (pi,P2,P3) hi K 6 . That is, the configuration space of three distinct point in 
the plane is part of a six dimensional space. In fact it is six dimensional, as almost 
all elements of R 6 can be seen as some configuration (^1,^2,^3) of three distinct 
points. 

We are interested in the Voronoi diagram of n points in the plane. That is, we want 
to define a Voronoi diagram for every configuration in a configuration space of n 
points. We are especially interested in possible Voronoi diagrams for point sets that 
contain coinciding points. We call such point sets degenerate configurations. The 
idea of a compactification of a configuration space is as follows: we want to construct 
a space that encodes all possible configurations of n distinct points, both non- 
degenerate and degenerate. Degenerate configurations should be on the boundary of 
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this space. By adding this boundary the degenerate configurations are incorporated. 
For a bounded space, adding the boundary is the same as compactifying, which 
explains the name 'compactification'. One important reason to compactify is the 
hope to be able to extend some definition, in our case the definition of Voronoi 
diagram, to the degenerate configurations. This extension would give access to 
limit objects which in our case are limits of Voronoi diagrams. 

We denote the set of ordered n tuples of all pairwise distinct points in the plane 
by CONF n . Although this is not very useful for us, an easy example is the com- 
pactification (K 2 ) n of the configuration space CONF n itself. The problem with 
this compactification is that it gives very little information on degenerate con- 
figurations: for two coinciding points p\ — P2, the only point that is added by 
compactifying is the point (pi,p2). Consider the two points pi(t) and pi{t) where 
P2(t) — Pi(t) + t(cosa,sinQ;). If t goes to zero, P2(t) will coincide with pi(t). De- 
scribe the Voronoi diagram of pi(0) and p 2 (0) as in Figure |L~3| . Then every distinct 
value of a corresponds to a distinct direction of the bisector of the points pi(0) 
and ^2(0). But this means that there is not one Voronoi diagram corresponding to 
Pi(0) — p 2 (0) but a complete collection of diagrams, parameterized by a mod 2tt. 



Compactifying using angles. 

By now, the attentive reader should be convinced that it is not a strange idea to 
use angles for a suitable compactification. In Chapter we may consider all angles 
between n points: we start with n distinct, labeled points in the plane. For every 
two points pi and pj with i < j we consider the directed line Uj passing through 
the two points. This line makes some angle onj mod 2tt with the horizontal axis. 
The angle map ipDA n maps a configuration of n distinct points to the («) directed 
angles with i < j. Let T be a triangle with vertices 1, 2 and 3. If we know 
the angles ai2, CH13 and 0123, we know in fact the shape of the triangle T. So, 
fpDA 3 (CONF 3 ) describes in a way all distinct triangles on three points. We describe 
explicitly this image ipvA 3 {CONF 3 ) in Chapter ||. We explain the consequences of 
the action of the symmetric group 53 on the labels of the three vertices. Finally 
we give geometric interpretations of boundary points of ipDA 3 (CONF3) in terms of 
degenerate configurations. 

The first compactification that we encounter is CDA n , the compactification of di- 
rected angles. This space is defined as the closure of the graph of the angle map 
ipDA^ , so a point 7„ € CDA n consists of n points and (™) angles between those 
points. If 7„ is a boundary point of CDA n , then not all points need to be distinct. 
In any case we know for every two points Pi,Pj G j n an angle c%ij and the position(s) 
Pi Pj or pi = pj. This implies that for 7„ a Voronoi diagram V(p/ n *) is defined. 

In Chapter |^ we prove a continuity theorem for Voronoi diagrams of data sets from 
CDA n . It states essentially that two data sets j n ,r] n £ CDA n that are close in 
the Euclidean metric, have Voronoi diagrams whose boundaries are close in the 
Hausdorff metric. The Hausdorff metric is very suitable to compare images. Two 
sets A and B in the plane are within Hausdorff distance r iff r is the smallest number 
such that any point of A is within distance r from some point in B and vice versa. 



An algebraic description: from angles to slopes. 

The line through the points (xi,yi) and (x 2 ,y 2 ) is given algebraically by the equa- 
tion 

y-yi = (x-xi) (1.2) 

x 2 - Xx 
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The quantity is of course the slope of the line. Instead of writing down the 

angle for every two points Pi and Pj , we consider all slopes . It turns out that 
there exists an algebraic relation between the six possible slopes for four distinct 
points po, Pi, P2 and p 3 . This relation is called the six-slopes formula, compare 
Figure |1.5| . The formula is given by A0123 = 0, where 

A0123 = («oi - ai 2 ){a 02 - a 23 )(a 3 - ai 3 ) - (a i - ai 3 )(a 02 - ai 2 )(a 3 ~ a 23 ). 

We prove in Chapter ^ that there also exists a triangle formula. This is a relation 
tijk = between slopes and x-coordinates of the points pi, pj, and pk- So by now 
we have two formulas involving points and slopes: one holds for any three points 
and the other holds for any four points. Suppose we have some configuration c of 
points and slopes such that all triangle formulas and all six-slopes formulas hold. 
A question that we consider in Chapter is the following: is such configuration c 
always a limit of non-degenerate configurations. To that end we introduce a variety 
that consists of exactly those configurations such that all triangle formulas and all 
six-slopes formulas hold. That is: the algebraic variety T n is just the set of zeroes 
common to all triangle formulas and all six-slopes formulas. It serves as an algebraic 
counterpart of the compactification CUA n of undirected angles. We prove that T n 
is not a smooth variety and give a geometric interpretation for the non-singular 
configurations. 




Figure 1.5: For four distinct points the six-slopes formula holds. 



A smooth and clickable compactification. 



Look at Figure 1.2 again. On the right, a limit Voronoi diagram of the points 
Pi, . . . ,ps is displayed. For every point pi a Voronoi cell is visible, except for p 3 . 
One could say that the Voronoi cell V(p^) is so small that we cannot see it. A 
solution in this case could be to rescale or magnify the picture somehow until the 
cell V(p 3 ) becomes visible. After magnifying enough we would get a picture that is 



very similar to Figure 1.2.b In more complicated configurations we might need to 
rescale or zoom in several times at several positions in order to distinguish every 
cell. 

This shows the need of a 'clickable' compactification on which we can define 'click- 
able' Voronoi diagrams. An example is given in Figure |L~6| . In the top 'screen', six 
points arc visible. But the two points in the middle arc in fact clusters of points. 
One cluster consists of the points p 3 , p-j and p\ 2 . This can be 'seen' by clicking 
on the cluster: The bottom left screen appears, displaying these points separately. 
The bottom right screen appears after clicking on the other cluster and displays the 
points p4, p 8 , pg, pw and pi\. 

The Fulton-MacPhcrson compactification is a well-known compactification of the 
configuration space of n distinct labeled points that has such a 'clickable' structure. 
To analyze Voronoi diagrams we do not need the full generality of this compactifica- 
tion. That is why we construct a real version of this compactification, incorporating 



11 





Figure 1.6: A clickable Voronoi diagram in a clickable configuration. 



some ideas of Kontsevich-Soibelman. Let c S CONF n be some configuration of n 
distinct points. Kontsevich and Soibelman write down both the angles a,j be- 
tween any two points, and, for every ordered triple of points (pi,Pj,Pk), the ratio 
Plj — |p*Zp fc | • They take the closure of the image of CONF n under the map that 
assigns all angles and ratios to a set of distinct points. This results in a manifold 
with corners FM2(n). 




Figure 1.7: The hook h\j — (/3^ , af* ) hinged at pi from pj to pt- The little arrow 
indicates the positive direction. 



In Chapter ^ we show how to adapt their approach so that the resulting compact- 
ification is a smooth manifold. Instead of just a ratio, we write down a hook for 
every triple of points. A hook hf- consists of this ratio /3jjr together with an angle 
afj between the legs of the hook, cf. Figure 1.7. We define a space XAH n as the 
closure of the image of CONF n under the map 

i>AH n '■ CONF n -> (R/ttZ)(2) x ((-oo,oo] x R/2nZ)/ ~ k ) 6 (°\ 

(p 1 ,...,p n ) i-> ((a y ) (PI* , a\j )), i, j, k pairwise distinct 

where ~fe denotes the identification of (/3y , afj ) with {—Pf- , a\j + n). 

Any point x G XAff n that is added to CONF n by taking the closure has some 
ratio Plj (x) equal to zero. We think of these points as corresponding to degenerate 
configurations, that is, configurations that have some coinciding points. We exploit 
these zero ratios to define a series of screens corresponding to x. A screen is just a 



copy of the plane R . In Figure 1.6 for example, there are three screens. We 'fill' 
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these screens with the degenerate configuration. We can do it so that any two points 
in the configuration are separated in at least one screen. Moreover, the process is 
such that it shows how to write XAH n locally as the graph of a function. This 
proves that XAH n is indeed a smooth manifold. 

The close connection between XAH n and the manifol d wit h corners FM2 has two 
immediate applications. First of all there is Theorem 7.65 that reveals the corner 
structure of FM2: the corners appear automatically if we describe the natural map 
from FM 2 to XAH n . 

Moreover, we can apply this same method of defining and filling screens to the 
manifold with corners FMi{n). That is, we add a family of screens to a degenerate 
configurations x £ FMi{n). But FM2(n) contains all angles mod 2ir between two 
points. Therefore we c an, without ambiguities, define a clickable Voronoi diagram 
as shown in Figure |l.6| for every x £ FMiin). We conclude Chapter [j] by showing 
how to jump between several descriptions and models in this thesis. Suppose we 
start with some x £ XAH n , for example. We can always construct a set S(t) of 
polynomial sites {p%(t), . . . ,p n (t)} with the following property: if we write down all 
angles and ratio's for the sites in S(t), we exactly obtain x. 



Comparing prices: higher order Voronoi diagrams. 

In what preceded we have extended Voronoi diagrams by limit Voronoi diagrams. 
Another extension or generalization of the classic Voronoi diagrams are the higher 
order Voronoi diagrams. We discuss some results on higher order Voronoi diagrams 
in Chapter [|. 

People in the Netherlands like to compare prices. That is, they have no fixed 
supermarket, but frequent the, say, five closest supermarkets to pick the bargains. 
For those people we have to divide the city into different sectors that we call 5- 
sectors. In one 5-sector, the five closest supermarkets are the same. So, if people 
know which 5-sector they live in, they know the five closest supermarkets. An 



example is given in Figure [La. In Figure 1.8. a the 6 supermarkets of above are 
displayed again. In Figure [Lq.b the four different 5-sectors are drawn. So, in 
sector 13456, the five closets supermarkets are supermarkets 1, 3, 4, 5 and 6. Or, 
if we put it the other way around: in 5-sector 13456, supermarket 2 is the most 
far away supermarket. This diagram is the fifth order Voronoi diagram of the six 
supermarkets. 




a b 

Figure 1.8: a. The Voronoi diagram of six supermarkets 1, ... ,6, and b. the fifth 
order Voronoi diagram of the same supermarkets. 



Note that in Figure 1.8.b only four distinct 5-sectors are present. This implies 
that two supermarkets (which ones?) are nowhere the most far away supermarket. 
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We can change however the position of the supermarkets in such a way that every 
supermarket is somewhere the most far away supermarket. (How to do this? Put the 
supermarkets on the vertices of a convex 6-gon.) This shows that if we change the 
position of the supermarkets, the set of 5-sectors may change as well. In Chapter [| 
we write down all 5-sectors. Moreover, we write down all 4-sectors, all 3-sectors and 
all 2-sectors. And we add all 1-sectors and all 6-sectors, although they are trivial 
(why?). If we collect everything we get a collection of subsets of {1, 2, 3, 4, 5, 6}. Or 
more general: we start with a set S of n distinct points in the plane. The Voronoi 
poset 11(5) consists exactly of those subsets of {1,. . . ,n} that have a non-empty 
Voronoi cell. 

Although the 5-sectors may change if we change the position of the supermarkets, 
there are certain invariants. An invariant is just a number that is independent of the 



position of the supermarkets. In Theorem 3.24 we prove for example the following: 
suppose the number of supermarkets is odd. Then the number of odd sectors equals 
the number of even sectors. Or more formal: for almost every configuration 5 the 
total number of cells in all even order Voronoi diagrams of 5 equals the total number 
of cells in all odd order Voronoi diagrams of 5. 



Chapter 1. Introduction. 



Chapter 2 

Voronoi diagrams 



In this chapter we give a short introduction to Voronoi diagrams. Most material 
covered in this chapter can be found in good books on computational and discrete 
geometry. A very readable general introduction to computatio nal geo metry is the 
book by De Berg, Kr evel d, Overmars and Schwarzkopf, see | BKO£ |. The book 
by Edelsbrunner, see |Ed|, is more in-depth, but less suited as an introduction. 
Okabe, Boots and Sugihara wrote a monograph on Voronoi diagrams for a broad 
audience, see | OBS | . Most recent is the chapter on Voronoi diagrams writ ten b y 
Aurenhammer and Klein in the Handbook of Computational Geometry (| AK |), 
Some concepts that are used in the rest of this thesis are shortly discussed in an 
overview paper on computational topology, cf . | DEG | . 

Apparently, Dirichlet in 1850 and Voronoi in 1908, cf. [Vo , were the first that used 
a notion of Voronoi diagram. They considered Voronoi diagrams of regular point 
sets, associated to quadratic forms. Consult [CF] for a recent and recommendable 
book on Voronoi diagrams of quadratic forms. 



2.1 Convex hull. 



A subset A of the plane is convex if for any two points p,q €z A the line segment pq 
is contained in A as well. The convex hull CH(A) of a set A is the smallest convex 
set containing A. Any two non-coinciding points p — (p xi p y ) and q = (q x , q y ) define 
two hull half-planes hh pq and hh qp where: 



hh v 



■ = {r = (r X) r v ) e M 2 | Bet(p,q,r) > 0}. 



The determinant Det(p, q,r) is given by: 



Det(p, q, r) 



1 Px Py 

1 q x q y 

1 "f X Tqi 



Let S be a finite set of distinct points in M. 2 . We can write CH(S) as an intersection 
of hull half-planes. More precisely, CH{S) equals the intersection of those half- 
planes defined by points in S that contain all points of S 



ch(s) = n 



,Pj£S 



As a non-empty intersection of half-planes, CH(S) is a convex polygon. We can 
represent this polygon by listing its consecutive vertices in clockwise order. In this 
way CH(S) corresponds to a cyclically ordered list of points from S. 
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Chapter 2. Voronoi diagrams. 




Figure 2.1: The convex hull of the points 1, . . . , 7. 



Example 2.1 

In Figure ^3), CH{S) <-> {2, 4, 5, 6, 7}. 

Suppose the points in S start moving around in the plane. Assume that at the 
start of the motion, no three points in S are collinear. A small enough disturbance 
of the points in S does not change the ordered list of points of CH(S). Only 
when a point pk passes a convex hull edge PiPj, the list changes (simultaneously 
Det(pi,pj,pk) — 0). We call this change in CH(S) a convex hull event. 



2.2 Voronoi diagram. 

Start again with a set S = {p\, . . . ,p n } of n distinct points in the plane. The 
Voronoi cell V(pi) of a point pi € S is defined as 

V{pi) := {<? e R 2 : dfa, q) < d( Pj ,q), i ? j}. 

Here, d(p, q) denotes the ordinary Euclidean distance between p and q. Note that 
we define a Voronoi cell as a closed subset of R 2 , in contrast to the choice made 
| [BKOSl- Thc Voronoi diagram V(S) of S is the family of subsets of R 2 



consisting of the Voronoi cells and all of their intersections. The boundary of a 
Voronoi cell consists of Voronoi edges and Voronoi vertices. A point q <E R 2 is 
on the Voronoi edge e{pi,pj) if d(q,Pi) — d(q,Pj) and d(q,Pk) > d(q,Pi) if k 7^ 
A point q £ K 2 is a Voronoi vertex if it is present on at least two Voronoi edges. 
A circle C is an empty circle with respect to S if there are no points of S inside 
the circle. For any three points pi, Pj, and pk, that are not collinear, there exists a 
unique circle Cijk passing through p,-, pj, and pk- A circle Cijk is a Voronoi circle 
if it is an empty circle. 




Figure 2.2: The Voronoi diagram of some points 1, . . . , 7. 



Lemma 2.2 

q € R 2 is a Voronoi vertex if and only if q is the center of a Voronoi circle. 



2.3. Topological changes. 
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Proof. |BK0S| 1, Theorem 7.4. ■ 

The bisector B(pi,pj) of two points pi and pj is the line equi- distant to Pi and pj. 
It is perpendicular to the line segment piPj passing through \{pi + Pj)- A point q 
is in the Voronoi half-plane vh(pi,Pj) if it is not closer to pj than to p^. 

vh(p t ,pj) = {q £l 2 | d(q,pi) < d(q,pj)}. 

As a consequence, the Voronoi half-plane vh(j>i,pj) is bounded by the bisector 
B(pi,pj). Any Voronoi cell is an intersection of half-planes 

V (Pi) = f] vh (phPj)- 

Therefore, any Voronoi cell V(j>i) is convex and is either bounded or unbounded. 
Lemma 2.3 

V{pi) is unbounded if and only if pi G SCH(S), where SCH(S) denotes the boundary 
of the convex hull CH(S). 



Proof. |OBS| , Property V2. ■ 
Lemma 2.4 

Let pi, pj and pk be three distinct points that are not collinear. Let c be the center 
of the circle through pi, pj and p^. Then: 

c = B{pi,pj) n B{pi,p k ). 
Proof. The point c is equi-distant to Pi,Pj and pk- B 



2.3 Topological changes. 




Figure 2.3: A diagram of type {134, 142, 245, 354}. 



Definition 2.5 

Let V(S) be a Voronoi diagram. Represent every vertex x of V(S) as an ordered 
list of labels of the points of S on the Voronoi circle C(x) of x. The order of the 
labels corresponds to the cyclic, clockwise order of the points on C(x). The set of 
all these lists is the type of the Voronoi diagram V(S). 



Example 2.6 

A diagram of type {134, 142, 245, 354} is depicted in Figure 12^3 
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A point set S is in general position iff no three points are collinear and no four 
points are cocircular. If S is in general position, then every Voronoi circle has 
exactly three points on its boundary. When the points in S start to move around in 
the plane, the type of V(S) changes exactly when the configuration of empty circles 
changes. This can happen generically in two ways. 




Figure 2.4: A circle event. 



The first is when two empty circles C a bd and Cbcd coincide. This is a circle event, 
see Figure 2.4. Before the event, a, b, c, and d define two empty circles C a bd and 
Cbcd- If c moves to the left in the leftmost figure, a, b, c and d become cocircular. If 
c continues moving left, one arrives in the situation of the rightmost figure, where 
Cabc and Cacd are the empty circles. The corresponding change in type is given by: 



{. . . , abd, bed, . . . } 



{. . . , abed, . . . } 



{. . . , abc, acd, . . .}. 






Figure 2.5: A convex hull event. 



The other way by which the type of a Voronoi diagram can change is by means of a 
convex hull event, see Figure |2.5| . Consider the circle defined by the points a, b and 
c. In the figure on the left, b 6 SCH(S). The circle C a bc contains all other points of 
S in its interior. Suppose b moves to the left. At some stage, b passes through the 
line segment ab. At this moment, the circle C a bc swaps over, and becomes empty, 
as in the picture on the right. This means that in the type of V(S), an extra term 
acb appears. 



2.4 Delaunay triangulation. 

Starting from a Voronoi diagram V(S) we define the Delaunay graph, DG(S), of 
S. The vertices of DG(S) are exactly the points in S. Two vertices p% and pj are 
connected by an edge in DG(S) exactly if there exists an edge e(pi,pj) of positive 
length in the Voronoi diagram V(S). Let D(S) be the straight line embedding of 
the Delaunay graph DG(S). For point sets S that are in general position, D(S) 
triangulates the convex hull CH(S). For this reason, the straight line embedding 
D(S) of the Delaunay graph DG(S) is called the Delaunay triangulation of S. 
This terminology is also used for point sets S that are not in general position. In 
the latter case D(S) is not necessarily a triangulation. An example of a Delaunay 



triangulation is presented in Figure 2.6 



2.5. Geometric transformations. 
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Figure 2.6: The Delaunay triangulation of the points introduced in Figure 2.2 



The Delaunay triangulation D(S) is dual to the Voronoi diagram V(S) in the 
following sense: vertices in the Voronoi diagram correspond to faces in the De- 
launay triangulation, while Voronoi cells correspond with vertices of D(S). As a 
consequence, the effect of a convex hull event and a circle event on the Delaunay 
triangulation D(S) is easily found. At the convex hull event depicted in Figure 



2.5 a triangle T a t, c appears in the Delaunay triangulation D(S). Moreover, for any 
triple of points pi, Pj, and pk in a point set S it follows that pi, pj, and p^ gives 
a triangle in D(S) if and only if the circle Cijk is a Voronoi circle. Indeed, in the 



circle event of Figure 2.4 an edge bd of the Delaunay triangulation flips over in an 



edge ac. This proves the following lemma. 
Lemma 2.7 

Let S be a set of distinct points in M. 2 in general position. The type of the Voronoi 
diagram V(S) is the list of triangles in the Delaunay triangulation D(S). 

The following characterization of Delaunay triangulations explains why Delaunay 
triangulations are often used in the generation of altitude maps and contour plots. 
The minimal angle of a triangulation is the smallest angle that occurs in any of 
the triangles of the triangulation. In Figure [2.6| , the minimal angle is the angle 
Z526. 

Lemma 2.8 

Let S be a point set in general position. The Delaunay triangulation D(S) is that 
triangulation that maximizes the minimal angle over all triangulations of S. 



Proof. [BROS], Theorem 9.9 



2.5 Geometric transformations. 
2.5.1 The lifting transformation. 

There is an easy correspondence between Delaunay triangulations of point sets in 
M. 2 and convex hulls in K 3 : the lifting transformation. This transformation is 
established by the map 

ip : R 2 ->■ R 3 , 

(x,y) i-> (x,y,x 2 + y 2 ). 



It maps a point in the plane onto a point of the paraboloid P defined by z = x 2 + y 2 , 
see Figure 2.7. The map ip has two properties that are important to us. 



20 



Chapter 2. Voronoi diagrams. 




Lemma 2.9 

Let C be a circle in the xy-plane. 

(i) ip(C) = P n He, where He is a plane in R 3 . 
(ii) A point q is inside C if and only if ip(q) below He- 

Proof. The proof follows directly from combining the equation of a circle in the 
plane and the equation of the paraboloid P. ■ 

Let S be a set of points in the a y-pla ne and S' the set of images of the points in S 



on the paraboloid P. By Lemma 2J] , ip maps empty circles defined by S onto faces 
of the lower convex hull of S'. Therefore, the projection of the lower hull faces of 
CH(S') onto the zy-planes gives exactly the Delaunay triangulation D(S). 

Let S be a set of n points in general position. A subset A of k points of S is a 
k-set if it can be separated from the complementary set B — S \ A of n — k points 
by a plane Va- Here we say that Va separates A from B if Va can be oriented so 
that all points in A are on the positive side of A, while all points of B are on the 
negative side. Lemma |2.9| shows that any set of k points contained in a circle in R 2 
can be mapped to a fc-set in R 3 . 



2.5.2 Mapping points to planes. 

A set of n planes V defines a subdivision of R 3 into connected pieces of dimension 
0, 1,2 or 3. This subdivi sion is the arrangement A(V) of V. Arrangements form 
a topic of their own, cf. []OT| . We assume that any point in any plane V, E V for 
i = 1, ... ,n can be written as (x, y, fa {x, y)), for some linear function fa : R 2 — ► M. 
That is, V consists of non- vertical planes only. We say that p — (p x ,p y7 p z ) 6 R 3 
is above plane Vi iff p z > fa(Px,Py)', similarly for below. The k-level of the 
arrangement -4(V) consists of those points in R 3 above or in k — 1 planes and below 
or in n— k planes of V. The n-level of an arrangement A(V) is also called the upper 
envelope of A(V). 

Using the map ip again, we show that we can determine the Voronoi diagram of 
a set S of points in R 2 by computing the boundary of an upper envelope in R 3 . 



Consider Figure 2.8 that demonstrates the principle one dimension lower. To any 
point p in a fixed point set 5* we associate a plane h p in R 3 . That is, we map p to 
the unique plane h p tangent to the paraboloid P at the point ip{p)- Lemma 2.1C 



shows that the set of planes {h p \p 6 S} completely encodes the relative distances 
of points q G R 2 to points in S. 



2.5. Geometric transformations. 
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Lemma 2.10 

Let q G K 2 and let h p (q) be the intersection of the vertical line through q and the 
plane h p . Then 

d{q,pj 2 = ip(q)-h p (q). 



Proof. An easy computation or [Ed], Observation 13.3 



! - 



Applying this result, we compute the Voronoi diagram V(S) as follows. Let q £ 
Suppose that q is in the interior V(pi)° of the Voronoi cell V(pi). By Lemma 2.10 
this is equivalent to h Pi being the first hyperplane that one encounters from the set 
{h Pj | pj G S} if one goes downwards in negative vertical direction, starting from 
ijj(q). But this means that V(S) is exactly the projection of the boundary of the 
upper envelope of the arrangement {h Pj | pj £ S} on R 2 . 

We can connect the Delaunay triangulation D(S) of a point set S with the lower 
convex hull of a set S'. Meanwhile the Voronoi diagram V(S) can be obtained from 
the upper envelope of an arrangement. Remember that the Delaunay triangulation 
D(S) is dual to the Voronoi diagram V(S). By now it may not come as a surprise, 
that there exists a duality transfer from lower convex hulls to upper envelopes and 



vice versa. For more information consult [Ed, BKOS 
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Chapter 3 

A Voronoi poset. 



Given a set S of n points in general position, we consider all fc-th order Voronoi 
diagrams on S, for k = 1, . . . , n, simultaneously. We recall symmetry relations for 
the number of cells, number of vertices and number of circles of certain orders. 
We introduce a poset II(S') that consists of the fc-th order Voronoi cells for all 
fc = 1, . . . , n, that occur for some set S. We prove that there exists a rank function 
on 11(5) and moreover that the number of elements of odd rank equals the number 
of elements of even rank of n(S'), provided that n is odd. 



3.1 Introduction. 

The dynamics of Voronoi diagrams in the plane is well understood. When n — 1 
points are fixed and one point is moving continuously somewhere inside the convex 
hull, combinatorial changes of the Voronoi diagram correspond to changes in the 



configuration of empty circles, see Chapter 1 and AGMR | for example. Changes 
in the configuration of non-empty circles correspond to combinatorial changes of 
higher order Voronoi diagrams. Here the fc-th order Voronoi diagram associates 
to each subset of size fc of generating sites that region in the plane that consist of 
points closest to these k sites. 

We consider all fc-th order Voronoi diagrams simultaneously for fc between 1 and n. 
We do so by introducing the Voronoi poset II (S) of a set S of n distinct sites in the 
plane. The poset consists of all sets of labels that correspond to a subset of sites 
that defines some non-empty Voronoi cell in some fc-th order Voronoi diagram. 

Higher order Voronoi diagrams have been investigated by numerous people. Many 



results are published in an article by D.T. Lee, see [Le|. A survey is given in 
Edelsbrunners book on algorithms in combinatorial geometry, sec [Ed|. In Chapter 1 
we have introduced fc-sets and we have shown that there exists a map ip that changes 
the point inside circle relation in R 2 into a point below hyperplane relation in R 3 . 
It turns out that these circles containing points serve as a 'b uildi ng block' for higher 
order Voronoi diagrams as we discuss in full detail in Section [3^. As a consequence, 
formulas counting fc-sets in R 3 can be applied in the counting of vertices, edges and 
cells of higher order Voronoi diagrams. Instead of considering circles that contain 
a fixed number of, say, fc points, one can also consider circles that contain at most 



fc points. This is done in |GHK] 



Let T be a set of n points in R in general position that are the vertices of a convex 



polytope. Sharir, [3h|, Lemma 4.4 and Clarkson and Shor, |Cq|, Theorem 3.5 prove 



that the number of fc-sets of T is given by 2(fc + l)(n — fc — 2). They prove this 
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formula using probabilistic methods that we do not discuss here. 




This formula can also be derived in the context of fc-th order Voronoi diagrams from 
Lee's results as has been observed b y several people, see again Clarkson and Shor, 
flCSfl or A ndrz ejak et al. , [AAHSWJ. We give this derivation explicitly and state in 
Theorem |3.18| that 



C-i ~\~ C ri 



2(< + l)(n-t-2), 



where Cj denotes the number of circles defined by a set S oi n points in general 
position in the plane, containing exactly i points of S. For an illustration, see 
Figure 3.1. Moreover we explicitly derive similar formulas for the number of cells 
fk in the fc-th order Voronoi diagram Vk(S), see Lemma 3.16, and the number of 
vertices Vk, see Lemma 3.17, in Vk(S). 



fk + Jn-fc+l 
v k + w„_ fe 



2k(n - k + 1) + 1 - 7i, 
4fc(n - k). 



These 'symmetry relations' are independent of the particular position of the sites 
in S, provided S is in general position: while the number of cells in some fc-th 
order Voronoi diagram may change, depending on the configuration, the sum of the 
number of cells in the fc-th order diagram and the number of cells in the (n— fc+l)-th 
diagram remains constant. 



In Section 3.3 we introduce the Voronoi poset mentioned above and prove that 



n(S') has a rank function. As an application of the symmetry relations we prove in 
Theorem 3.24 that the number of elements of odd rank in T1(S) equals the number 
of even rank, provided that n is odd. 

The Voronoi poset of a set S of n moving points seems a natural object to study 
as changes of the poset occur exactly at those moments where S is not in general 
position. As there are tight connections between higher order Voronoi diagrams, 
fc- levels in certain arrangements in R 3 and certain fc-sets in R 3 , the study of the 
Voronoi poset may have applications in these areas as well. 



3.2 Higher order Voronoi diagrams. 

3.2.1 Definition of k-th order Voronoi diagram. 

Let S = {pi, . . . ,Pn\ be a set of n distinct points in R 2 in general position. Let 
< fc < n. For every point p in the plane we ask for the fc nearest points from S. 



3.2. Higher order Voronoi diagrams. 
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That is, we look for a subset A C S, such that 

\A\=k, VieA, Vy€S-A: d{p,x) < d{p,y). 
For two points in R 2 , we define a half-plane 

h(x,y) := {p£R 2 \d(x,p)<d(y,p)}. 
The Voronoi cell of A C S of order |^4| is the intersection of half-planes 

V(A) := f| h(x,y), 

x£A, y£S-A 

whenever this intersection is not empty. As an intersection of half-planes, V(A) is 
a convex polygon. 

Remark 3.1 

It is left as an exercise to the reader to show that assuming general position implies 
that a Voronoi cell is not a line segment or a single point. 

The k-th order Voronoi diagram is the subdivision of M 2 , induced by the set 
of Voronoi cells of order k. For later purposes, we identify the fc-th order Voronoi 
diagram with the set of non empty /c-th order Voronoi cells. 

V k (S) := {V(A) | A c S, \A\ = k, V(A) ^ 0}. 






a. b. c. 

Figure 3.2: A first, second and third order Voronoi diagram. 



Example 3.2 

Let S = {(45,86), (76,40), (40,42), (1,9)}. Figure |J[a. shows the first order Vo- 
ronoi diagram. The second order diagram is shown in Figure 3.2.b. and the third 
order diagram in Figure 3.2.C. All non-empty Voronoi cells are indicated by their 
generators. Note that not all possible triples occur. 

Remark 3.3 

A planar graph that represents a point-face dual of the fc-th order Voronoi diagram 
can be constructed as follows, cf. [AS]. Write down for every A C S with \A\ = k 
and Vk(A) ^ its centroid c(A), defined by c(A) = (1/fc) J2 p£ aP- Two centroids 
C(A) and c{B) are connected by an edge exactly iff Vk{A) and Vu{B) share an edge. 



3.2.2 Circles and higher order Voronoi diagrams. 



In this section, we state some elementary properties of higher order Voronoi dia- 
grams. Every edge in Vk(S) is part of some bisector B(a,b), with a, b E S. The 
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Voronoi vertices are exactly those points that are in the centers of the circles deter- 
mined by three points from S. Therefore, under our general position assumption, 
every Voronoi vertex has valency three. The following theorem describes the local 
situation around a Voronoi vertex. The symbol Q a b denotes the circle passing 
through the points a, b, and c. 

Theorem 3.4 

Let x be the center of Q a b , for a, b, c S S, let 

H = { z e S | d(x, z) < d(x, a) }, 

and let k = \H\. Then x is a Voronoi vertex of Vk+i (S) and (S) . The Voronoi 
edges and cells that contain x are given in Figure \3.3( . Moreover, all Voronoi vertices 
are of this form. 

b a 



V(H U {a,b» 




Voronoi edge of V k+] (S) 



V(H U (c)) Voronoi edge of V k+2 (S) 

I 
I 

b'a 



Figure 3.3: The Voronoi diagram around x. 



Proof. [Dc|, Theorem 1 and Theorem 2. 



Let a,b,c and H be as defined in Theorem 3.4. We define the order of a circle 
O a ,b,c as 1^1- Notation: \Q abc \ := \H\. An order k Voronoi circle Q abc 
is a circle through three points a, b and c from S that contains exactly k points 
from S — {a,b,c}. In fact, from all (™) Voronoi circles Q abc and all sets H. a ,b,c, 
compare Theorem |3.4|, almost enough information is provided to construct all fc-th 



order Voronoi diagrams VhiS) for k = 1, . . . , n — 1. 

Algorithm 3.5 

Voronoi diagrams of all orders. 

Input: set S of n points in general position. 

Output: all fc-th order Voronoi diagrams Vf-(S) for k = 1, . . . , n — 1. 
l: Compute all circles © a b c defined by S. 
2: Compute all sets H a ,b,c defined by S. 

3: Take all circles of order k — 1 and order k — 2. The centers of these circles are 
exactly the vertices of Vk(S). 



Theorem 3.4 gives for every vertex the three incident edges and the three inci- 
dent cells. 



3.2. Higher order Voronoi diagrams. 
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Two vertices are connected by an edge iff the two vertices have two incident 
cells in common. Skip the edge if it is used. 

Edges that are not skipped are unbounded edges. Their direction and orien- 
tation still have to be computed. The direction is simply the direction of the 
bisector containing the edge. The orientation follows from Figure |3.3|. 



Remark 3.6 

Denote the number of circles of order k by c k and the number of vertices in a k-th 



order Voronoi diagram by v k . As a consequence of Theorem 3.4 we get 

Vk = c fe _i + c fe _2- (3.1) 



3.2.3 Counting vertices, edges and cells. 

The following theorem shows that the total number of vertices, edges and Voronoi 
cells does not depend on the positions of the points in S, assuming general position. 

Theorem 3.7 

Let Vk, ek, and f k denote the number of vertices, edges and cells in V k (S) for some 
set S of size n in general position. The total number of vertices, edges and cells in 
the Voronoi diagram of all orders are as follows. 

ft) Efe=i y k = \ n ( n ~ - 2 )- 
(H) Efc=i e * = \n{n- If. 
(ni) ELiA = \n(n* + b). 

Proof. We prove the three claims. 

(i) Every circle center defined by three distinct sites from S is a Voronoi vertex 
in some k-th and (fc + l)-th order Voronoi diagram. As there are Q) distinct circles, 
the first claim follows. 

(ii) Consider the arrangement of bisectors A(S). Fix one bisector B(a, b). As S 
is in general position, we may assume that the bisector B(a,b) is divided into n — 1 
line segments by the Voronoi circle centers abx^,abx^, . . . ,abx n , where we write 
S = {a, 6, £3, . . . , x n }. Every line segment is an edge in some k-th order Voronoi 
diagram. As there are distinct bisectors, claim (ii) follows. 

(iii) The Euler formula, Vk — e& + //- = 1, holds for every order. Therefore 

n n n 

k=l k=l k=l 

which completes the proof. ■ 

The number of vertices, edges and cells in Vk(S) depends on the configuration of 
S as the ordinary Voronoi diagram shows. The following theorem gives expressions 
for these numbers. Let fjf denote the number of unbounded cells in the fc-th order 
Voronoi diagram. By definition /q° := 0. 

Theorem 3.8 

Let S be in general position. Then the number of vertices, edges and cells in the 
k-th order Voronoi diagram can be expressed as follows. 
(1) v k = 2(/ fc - 1) - 

(ii) e k = 3(/ fc -l)-/ fc ~ 

(iii) f k = (2k - l)n - (P - 1) - E-=i fi-v 



Proof. |Ed|, Q. 
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Note that /„ = 1. Substituting k = n in the expression for in Theorem 3.8 yields 
the following equation for the total number of unbounded cells: 

n 

J2f™i = n(n-l). (3.2) 

The unbounded cells in the fe-th order Voronoi diagram can be characterized as 
follows: let pq denote the line segment with endpoints p and q and l pq the line 
through p and q. 

Property 3.9 

A cell V(A) of the k-th order Voronoi diagram Vk(S) is unbounded if and only if 
one of the following two conditions holds, 
(i) There exists a line I that separates A from S — A. 

(ii) There exist two consecutive points p and q, withp,q € S — A, on 5CH(S — A) 
such that the points in A — pq are in the open half plane defined by l pq opposite to 
CH{S-A). 

Proof. flOBSfl , Property OK4. ■ 
Under the general position assumption, we only need to consider condition (i) in 



Property 3.9. It is clear that in this case the following symmetry holds: 

fk = fn-k- (3-3) 



3.2.4 Circle events and hull events. 



In Section 2.3 we have discussed the type of a Voronoi diagram in connection with 
circle events and convex hull events. We generalize these notions to fe-th order 
Voronoi diagrams. Let S be a set of points in general position and fix some k € 



1, . . . , n. As in [Le we call a Voronoi vertex of Vk(S) old if its corresponding circle 
has order k — 2 and new if its corresponding circle has order k — 1. We will also use 
the words old and new to indicate order k — 2 and order k — 1 circles with respect 
toV k (S). 

Example 3.10 

In the classic or first order diagram, vertices correspond with empty circles, so all 
vertices are new. 

Represent every vertex x of 14 (S) as an ordered list of labels of the points of S on 
the k — 1 or k — 2 order Voronoi circle C(x) corresponding to x. Order the labels 



as in Definition 2.1. The type of a fc-th order Voronoi diagram Vk(S) consists of 
two lists: a list of ordered labels of old vertices and a list of ordered labels of new 
vertices. 

Example 3.11 

Consider the configuration in Figure |3.5| .a. The circles 124 and 143 are empty, while 
the circles 123 and 243 both contain one point. Therefore, the type of the second 
order Voronoi diagram is given by: {{124, 143}, {123, 243}}. 

Suppose that the points in S start moving. S is in general position, as long a no 
four points are cocircular and no three points are collinear. If four points become 
cocircular, then two circles C a td and Cbcd, both surrounding a subset TcSof size 



3.3. The Voronoi poset. 
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\T\ = fc, coincide and change into two circles C a b c and C ac d, see Figure 2.4. This 
corresponds to a change of the Voronoi vertices abd and bed into Voronoi vertices 
abc and acd in both Vfc + i(5 l ) and V)~±2(S). This is the generalization of a circle 
event and is called a k-th order circle event. So, a zero order circle event is just 
a circle event. 



k a 4 a ^ a 

i b n k 3 k T b n-k-3 k b * n k 3 

i i 1 

/ c + c ; c 



Figure 3.4: A fc-th order convex hull event. 



If three points become collinear , an order fc circle abc changes into an order n — k — 3 



circle abc, see Figure |3.4| . This implies that the vertex abc disappears in Vk+i{S) 
and Vk+2(S) and appears in V n -k-i(S) and V n -k-2(S). This is the generalization 
of a hull event and is called a k-th order hull event. So, a zero order hull event 
is just a hull event. At a fc-th order circle event the type of Vk+i(S) and 14+2 (>S) 
changes. At a fc-th order hull event, the type of Vk+i(S), Vk+2(S), V n -k-i(S) and 
V n -k-2(S) changes. 



3.3 The Voronoi poset. 
3.3.1 Definition and examples. 

Fix a labeling of the sites in S and identify a set of sites A C S that defines a 
non-empty Voronoi cell V(A) with the set of labels L(A) C [n] of the sites in A. A 
subset L of [n] may or may not correspond to some Voronoi cell V(Al). For fc = 1 
we retain the ordinary Voronoi diagram, which implies the correspondence 

Vi{S) - {{l},{2},...,{n}}. 

We define V (S) = {0}. The set {{1, . . . ,n}} corresponds to V n (S). We consider 
the set of all Voronoi cells that appear for a given set S of points and call the set 
of corresponding labels the Voronoi poset 11(5) of S: 

U(S) := |J { L(A) | V(A) e V k (S) }. 

k 

This definition also makes sense when we drop the general position assumption. 

We order the elements in the poset by set inclusion of the sets L{A). This yields 
a partially ordered set. For more on partially ordered sets consult |ZJ. The poset 
is bounded since we have the empty set as 0, the unique minimal element, and the 
set [n] as 1, the unique maximal element. In general, a poset is called graded if 
it is bounded and if every maximal chain has equal length. We show that II(S') is 
graded. Below we give an example showing that 11(5') is in general not a lattice. 

Property 3.12 

II(S') is graded. 
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Proof. We show that r(L(A)) = \L(A)\ is a rank function for 11(5). A rank 
function maps an element x from a poset to a unique level in such a way that the 
level corresponds to the length of any maximal chain from x to 0. Let L(A) £ 11(5), 
with |L(j4)| = k. Every point x € V(A) has all k points from A as its k nearest 
neighbors. Order those points with respect to their distance to x. As we assumed 
general position it is always possible to change the choice of x in such a way that 
this order is strict. By removing at each step the furthest point still available, we 
obtain a chain of length k that descends to 0. ■ 

We analyze the two smallest cases, assuming general position. 
Example 3.13 

For n — 3 there is one poset, the full poset on [3]. That is, 
n 3 (5) = {0,1,2,3,12,13,23,123}. 

Example 3.14 

For n = 4 there are two essentially distinct posets as is evident from the circles 
defined by four points. Since n = 4, Voronoi circles are of order one or two. Eulcrs 
formula rules out four circles of order one. At the same time less than two circles 
would not yield enough cells in the first order diagram. 




Figure 3. 5. a has two empty Voronoi circles. Clearly, subset 23 is the only subset of 
[4] missing, leaving us with : 



n 4 (5i) = {0, 1, 2, 3, 4, 12, 13, 14, 24, 34, 123, 124, 134, 234, 1234}. 

This example shows that the Voronoi poset is in general not a lattice. A lattice 
requires that every two elements of the poset have a unique minimal upper bound. 
In this example, the elements 2 and 3 have two minimal upper bounds, namely 123 
and 234. 



In Figure 3.5.b. there are three empty Voronoi circles. The cell 123 cannot appear 
in the third order diagram, but all other subsets of [4] do appear, thus: 



n 4 (5 2 ) = {0, 1, 2, 3, 4, 12, 13, 14, 23, 24, 34, 124, 134, 234, 1234}. 



3.3.2 The order complex of the Voronoi poset. 

The standard way to associate a topological space to a finite poset (P, <) is by 
means of the order complex A(P) of the poset, see El, Wa|. The order complex 



3.4. Symmetry relations. 



31 



is the simplicial complex of all nonempty chains of P. A chain of P of length k is 
a totally ordered subset 

X < Xi < X 2 < ... < Xfc, 

of elements x E P. The well-known geometric realization associates a topological 
space with a simplicial complex. 

As a Voronoi poset on a set of n points p\ , . . . , p n always has a unique maximal 
element {1, . . . , n}, the geometric realization of the order complex is a cone and 
therefore contractible. This shows that the topological space that we have associated 
with S is homotopy equivalent with a point, and therefore not very interesting. 

More promising is to consider the complement, that is the anti Voronoi poset 
aP(S), consisting of those subsets of {1, . . . , n} that are not in the Voronoi poset. 
Another possibility is to consider the arrangement of bisectors. 



3.4 Symmetry relations. 

Given a set S of sites, we count for every order k the number of vertices i>/-, the 
number of edges e& and the number of non empty Voronoi cells /&. The f-vector 
of 11(5) is the vector {/i, f 2 , ■ ■ . , f n }- The c- and e-vector are defined analogously. 
Note that the /-vector of n(S') may change if the position of the sites in S changes. 

Example 3.15 

Consider the two configurations Si and £2 on four points, presented in Example 
p\l4 The /-vector of II(Si) equals (1,4, 5, 4, 1), while the /-vector of U(S 2 ) equals 



(1,4,6,3,1). 



3.4.1 Symmetry in the number of cells. 

It turns out that a symmetry exists in the /-vectors. 
Lemma 3.16 

Consider the f-vector of H{S), where \S\ = n. Then fk + fn-k+i is a constant 
independent of the position of the points in S . More precisely, 

fk + fn-k+i = 2fc(n-fc + l) + l-n. (3.4) 



Proof. We apply Theorem 3.8 to fk and f n -k+i'- 

k 



fk + fn-k+i = (2k-l)n-k 2 + l-J2f^i 

i=l 

+ (2( ? i- £; + !)- l) ? i- (n-fc + 1) 2 + 1- ^ f£i, 



n-k+l 



= 1 



= 2kn-2k 2 + 2k + l-n + n(n-l)-(Y,f^i+ E /~l )• 

i=l i=l 

We join the two sums by applying Symmetry Equation and evaluate the result 
by using Equation p\2| . 

k n—k+X n 

E/-i + E /~i = E/-i = n(n-l). 

i—1 i—1 i—1 

The lemma follows from combining the two equations above. ■ 
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3.4.2 Symmetry in the number of vertices. 

A similar equation holds for the number of vertices of a collection of Voronoi dia- 
grams Vfc(S), for k = 1, . . . ,n — 1. 

Lemma 3.17 

Let S be a set of points in general position with \S\ = n. Let Vk denote the number 
of vertices in the k-th order Voronoi diagram. Then: 

Vk + ^n-fc = 4fc(n — k) — 2n. (3-5) 



Proof. Using Theorem 3.8 we write v k + v n -k m terms of numbers of cell s. N ext 
we regroup and apply Symmetry Equation 3.3. After applying Theorem 3.S we 
combine using symmetry again. Finally, using X)T=i fi^i = n ( n ~ -0 completes the 
proof. 



v k + v n - k = 2(/ fc - 1) - A?° + 2(/ n _ fc - 1) - /~ k , 
= 2(/fe + fn-k - 2 - 

/c n—k 

= 2(n 2 -2n + 2kn 2k 2 + ^ + / fe °°)), 



4=1 



= 2(n 2 -2n + 2kn -2k 2 /i^l)> 



— 2n + 4fcn - 4fc 2 , 
4fc(n - fc) - 2n. 



3.4.3 Symmetry in the number of Voronoi circles. 

Recall that the order of a Voronoi circle equals the number of points of S contained 
in its interior. We define the c- vector of S as the vector 



c(S) = {c ,ci, . . . ,c„_ 3 }, 



where Cj denotes the number of circles of order i. The following theorem states that 
for n arbitrary points in general position, the number of circles containing exactly i 
points on their inside plus the number of circles containing exactly i points on their 
outside is constant. We prove this by applying the above results. 

Theorem 3.18 

Consider the c-vector of II(S') ; where \S\ = n. Then Cj + c„_i_3 is a constant 
independent of the position of the points in S . More precisely, 

Ci + Cn-i-3 = 2(i + l)(n-2-i), (3.6) 
= 2i(n-i-3) + 2(n-2). 



Proof. We prove the theorem by induction. 

[i=0]. We use the lifting transformation. This transformation changes the 
point-inside-circle relation in 2-dimensional space in a point-below-plane relation in 
3-dimensional space, 
defined by 

A, ■ TB2 , Tn>3 



See also Section 2.5.1. The lifting transformation map <j) is 



(x,y,x 2 +y 2 ). 
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It lifts points in the plane to the unit paraboloid in three-space. As every circle 
defined by S in the plane contains only three points from S, every hyperplanc 
defined by 4>(S) contains only three points from <f>(S) as well. The number cq of 
empty circles of S in the plane equals the number of facets of the lower hull of 4>(S) 
in three dimensions. At the same time, the number c n -3 of circles that contain all 
other points of S equals the number of facets of the upper hull of 4>(S). All images 
of points in S under <f> are part of the convex hull of <fi(S). Since the convex hull 
of a point set of n points consists of 2n — 4 facets, if every facet is a triangle, see 
BKOSfl , Theorem 11.1, the claim follows. 

[induction step] . We deduce the expression for c% + c n -k-3 by applying Equation 
JO, followed by combining Lemma 3.17 and the induction hypothesis: 



Cfc + C„_fc_ 3 = C fc _i + Cfe + C„_fe_ 3 + C„_fe_ 2 - (Cfc_i + c„_ fc _ 2 ), 

= Vk+i + t'n-o+i) — (cfc-i + c„_ fe _ 2 ), 
= 2 (2 (ft + 1) - l)(n - (ft + 1)) - 2(ft + 1) 

-(2(fc-l + l)(n-2-(ft-l))), 
= 2(ft + l)(n-2-fc). 



Let fk ■= fk + fn-k+i an d c~i :— Ci + c„_j_ 3 . By the reduced f-vector, denoted 
/, we mean the vector of fk's for all distinct ft. That is 

/ : = {/o, fl, ■ ■ ■ j /[n^ij}- 

c is defined similarly. As a consequence of Lemma |3.16 and Theorem 3.18 , / and c 
are only dependent on n. 

Example 3.19 

As an example we present the reduced /- and c- vectors for n G {3, . . . , 12}. 



n 


/ 


c 


3 


(4, 6) 


(2) 


4 


(5, 9) 


(4) 


5 


(6, 12, 14) 


(6, 8) 


6 


(7, 15, 19) 


(8, 12) 


7 


(8, 18, 24, 26) 


(10, 16, 18) 


8 


(9, 21, 29, 33) 


(12, 20, 24) 


9 


(10, 24, 34, 40, 42) 


(14, 24, 30, 32) 


10 


(11, 27, 39, 47, 51) 


(16, 28, 36, 40) 


11 


(12, 30, 44, 54, 60, 62) 


(18, 32, 42, 48, 50) 


12 


(13, 33, 49, 61, 69, 73) 


(20, 36, 48, 56, 60) 



Remark 3.20 

Computer calculations did not suggest any similar symmetry relation for the number 
of edges. 



3.4.4 Relations between cells and circles. 

Corollary 3.21 

fi= fo + Cf-i = Cf-i + n + 1. 



Proof. This follows directly from Lemma 



3.16 



and Theorem 13.18 
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Property 3.22 

Let f°° denote the number of unbounded cells in the i-th order diagram and let Cj 
denote the number of circles of order i: 

f?° + (ci_! - Ci _ 2 ) = 2(n-i). (3.7) 

Proof. We prove the property by induction. 

[i=l] c_i is zero by definition. The number of vertices v\ in the first order 
Voronoi diagram equals the number of circles of order zero, cq. The claim follows 
from applying Theorem |3.8| : 

/r° + ( Co - c -i) = /r+«i, 

= /r + 2(/!-i)-/r, 

= 2(n-l). 

[induction step] Assume we have proved that 

/f + (ci-i - Ci_ 2 ) = 2(n-i)« 
We rewrite this, by using induction, as 

»+i 

c»_i = 2ni-i(i + l)-5^A_i. (3.8) 

fc=i 

Evaluate Cj — Cj_i: 

c, - Cj_i = (cj + Ci_i) - 2ci_i, (3.9) 
= - 2cj_i, 
= 2(/ <+1 -l)-/£ 1 -2c j _i. 

Substituting this expression for Cj — Cj_i and applying Theorem [3.^ and Equation 
3.8 proves the claim: 

+ = 2(/i+! - 1 - <*_!), 

= 2(n-i-l). ■ 

Corollary 3.23 

T/ie c-vector totally determines the f -vector. The correspondence is given by 

fk = n-fc+l + Cfc_ 2 . 



Proof. Applying Equation 3.7 we get 

k 



Y^f^i = (A-l)(2n-fc)-Pfc_ 2 . 



i=i 



The claim follows from evaluating Theorem 3.S using the expression above. 



3.5 Even versus odd order cells. 

Given a grading on a set of objects, it is common to consider the Poincare polyno- 
mial P(t) of the grading. The i-th coefficient of this polynomial equals the number 
of objects of grade i. In our case, the objects are the elements of the Voronoi poset 
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II(S'), while the grading is given by the rank function on the poset. Recall that the 
rank of an element x in 11(5) is just the order k of the Voronoi diagram in which 
x occurs as a cell. The i-th coefficient of the Poincare polynomial P(t) is given by 
fi, as fi gives the number of cells in the i-th order diagram Vi(S). So, the Poincare 
polynomial P(t) of H(S) with respect to our rank function is given by 

P(t) = f + flt + f2t 2 + ---+fnt n . 

As an application of the symmetry relations we compare the number of cells in the 
even order Voronoi diagrams with the number of cells in the odd order diagrams. 
In terms of the Poincare polynomial P(t) of above, the following result can also be 
formulated as 

P(-l) = 0. 

Theorem 3.24 

Let S be a set of points in general position with \S\ = n > 3. Assume n is odd. In 
this case, the number of cells in the even order Voronoi diagrams equals the number 
of cells in the odd order Voronoi diagrams. 

Proof. Write fa = fi + f n -i+i- We show that ^4 = 0, where: 
A = -fo + fi-h + •••-/„-! + /„. 

So A is the number of cells in the odd order diagrams minus the number of cells in 
the even order diagrams: 

A = -f + fi + -/«±i + t n , 

2 2 

where 



Clearly, /o = 1, as fo counts the empty set. /i is the number of points in S plus 
the number of cells in V n (S), so /i = n + 1. Applying Equation |3.4| gives: 

/a±l = -(-1) 2 —7— 



n 2 + 3 



Straightforward calculations show that: 

tn = (-1) 

from which it follows that: 

1 / lX i±i"- 2 + 3 . , , n+i n 2 + 3 
A = + n + l -(-!)— — - — + (-!)— — n = 0. 



The claim of Theorem 3.24 does not hold when n is even. However, the following 
result does hold. 
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Lemma 3.25 

Let S be a set of points in general position, with \S\ = n > 3. Assume n is even. 
Let A(S) denote the number of cells in the odd order Voronoi diagrams minus the 
number of cells in the even order diagram. Then: 

n = 0(4) => A(S) odd. 
n = 2(4) => A(S) even. 



Proof. Similar computations as in the proof of Theorem 3.24. ■ 
Note that as Vk = Ck-i + Cfe-2 it follows immediately that: 

n-l 

£(-l) fe +V = 0, 
fc=l 

for all n, where Vk denotes the number of vertices in the fc-th order Voronoi diagram. 



Chapter 4 

Limits of Voronoi diagrams. 



The Voronoi diagram of a set S of n distinct points in M 2 associates to a point 
p € S that part of the plane that is closer to p than to any other point in S. In this 
chapter we assume that the position of any of the n points in S = S(t) is given by a 
pair of polynomials in one parameter t, such that no two points are represented by 
the same pair of polynomials. We call the elements of S(t) polynomial sites. In this 
setup it is possible that sites coincide at, say, t — 0. We define a Voronoi diagram 
V(S(0)) := lim m V(S(t)). That is, V(S(0)) is defined as a limit diagram in a 
particular sense of the ordinary Voronoi diagram of the positions, at small positive 
t, of the polynomial sites. We show how to extend the notion of type to polynomial 
sites. This enables us to determine the combinatorics of V(S(0)). It turns out that 
in general some sites can be omitted without changing the boundary of the Voronoi 
diagram V(S(0j). We characterize those sites that can not be omitted and present 
an efficient algorithm to determine these sites together with the boundary of the 
Voronoi diagram. 



4.1 Introduction. 

Whenever Voronoi diagrams are studied, it is assumed that all points defining the 
diagram are distinct. We call these points sites. In the case of dynamic Voronoi 
diagrams, where the sites are moving continuously over time, the assumption above 
means that sites are not allowed to coincide at any moment. In this chapter we 
investigate in one particular setting what happens when we do allow sites to coincide 
in the plane. 

We consider a set 

S(t) = { Pl (t),..., Pn (t)}, 

of n sites in the plane such that the position of site pi at time t is given by a pair 
of polynomials in t, one for every coordinate. That is, both the movement of the x- 
and the y-component are described by a polynomial in R[t]. 

In fact we do not have to restrict ourselves to sites described by pairs of polynomials: 
the theory developed in this chapter works whenever the movement of both the 
x- and y-componcnt of all sites is described by functions that can be expanded 
as a convergent Taylor series. For the sake of simplicity we restrict ourselves to 
polynomials however. 

In Chapter 1 we have described the Voronoi diagram of a set of distinct sites both 
in terms of half-planes and of empty circles. To check whether some point q is 
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in a given half-plane h(pi,p2) or inside a circle c{p\ 1 p2 1 p 3 ) for p\,p2,pz G R 2 
boils down to evaluating the sign of an easy polynomial expression //, (pi ,P2,q) or 
fc(PiiP2,P3, q)- The main idea of this chapter is to replace pi for i = 1, 2, 3 and q 
by Pi(t), q(t) € M[<] 2 and to evaluate the sign of the coefficient of the lowest degree 
term in f h (pi(t) , p 2 {t) , q(t)) and f c (p 1 (t),p 2 (t),p 3 (t),q(t)), at t = 0. In this way 
we define a Voronoi diagram V(S(t)) at i = 0, even if sites do coincide at t = 0. 
For a set of sites that coincide at t = we regard V(S(0)) as the limit diagram 
lim^o V(S(e)). We develop this arithmetic for polynomial sites in all details in 



Section 4.2 



Here we add a warning: we can define a limit Voronoi diagram at t = that is 
consistent with small positive t or small negative t. In general these two approaches 
give distinct results. Examples will be presented in the main text of this chapter. 

There is a similarity between our approach and a well-known technique in compu- 
tational geometry that is used to avoid computations involving degenerate input. 
In this technique, a set of distinct points that is not in general position is perturbed 
slightly in such a way that it is in general position after perturbation. For example, 
three points that are collinear before perturbation will not be collinear anymore 
after perturbation. As a consequence it is possible for a generic algorithm, that is, 
an algorithm that can only handle points in general position, to handle the input 
and compute the wished geometric structure, e.g. a Voronoi diagram. Such generic 
algorithms are in general much easier than algorithms for arbitrary input as a lot 
of degenerated cases can be ignored. In the so-called Simulation of Simplicity tech- 



nique, cf. [EM, Ed2], polynomials in the variable e are added to points of a set of 
distinct, degenerated points in such a way that evaluation of the new points for 
small enough e produces a point set that is in general position. An overview on 



robust geometric computation is given in [ Ya | . 

In Section |4~3| we extend the definition of type, see Definition 2.5, to a set of poly- 
nomial sites S(t) at t = 0, provided that S(0) fulfills some general position assump- 
tions. This type gives us a complete combinatorial structure that matches with the 
combinatorial structure of the Voronoi diagram V(S(t)) for small enough positive t. 

In Section |4.4j we define a Voronoi diagram for n not necessarily distinct points in 
the plane and (™) angles between those points. We apply this definition to introduce 
Voronoi diagrams for a set of polynomial sites S(t) at t = 0. The resulting Voronoi 
diagram or polynomial sites matches with the ordinary Voronoi diagram of the 
positions of the sites at small positive t. In this section we also give the connection 
between half-planes and Voronoi circles defined by S(t) at t = 0, thereby connecting 
the notion of type to polynomial sites diagrams. 

The shape of a Voronoi diagram is defined as the union of the boundaries of the 
Voronoi cells. It turns out that in general some sites in a set S(t) of polynomial 
sites can be omitted without changing the shape of the Voronoi diagram of S(t). 
A question that we pose is how to determine this shape efficiently and how to 
characterize those sites that do determine the shape. We answer this question by 
splitting up the problem into two parts. 



First we consider in Section 4.5 the shape of the Voronoi diagram of one cluster of 



polynomial sites at t = 0. By an /-cluster we mean a set of polynomial sites such 



that the positions of all sites coincide at one location I at t = 0. Lemma [4.35| fully 
classifies the sites Pi(t) in the cluster such that &re&(V (pi(0))) > at t = 0. 



The second part of the question is solved by Lemma 4.41. This lemma states that 
the shape of the Voronoi diagram at t = of an arbitrary set of polynomial sites 
can be found as follows. First, compute the ordinary Voronoi diagram of all distinct 
locations at t = 0. Second, plug in the cell of location I, the shape of the Voronoi 
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diagram of the ^-cluster. 

Section 4.6 shows that a lot of combinatorics can be hidden in the edges of the 
shape of the Voronoi diagram. In Section 4.7 we demonstrate the theory developed 
in the former section by a somewhat bigger example. We conclude this chapter in 
Section 4.9 with some remarks on generalizations to fc-th order Voronoi diagrams 
and shortcomings of this setup. 



4.2 Preliminaries. 

We define half-planes and circles for polynomial sites. A polynomial site p(t) = 
(p(t)x,p(t) y ) consists of a pair of polynomials p{t) x ,p(t) y € R[t]. Two polynomial 
sites u(t) and v(t) are called distinct if they represent distinct elements in l[f]xM[f], 
Throughout this section assume that S(t) = {pi(t), . . . ,p n (t)}, where every pi(t) a 
polynomial site and Pi(t) distinct from Pj(t) as a polynomial site whenever i ^ j. 

Remark 4.1 

The notions in this section are introduced so that they match the situation for 
ordinary points in the plane obtained by substituting a very small positive value of 
t in the set S(t) of polynomial sites. 



4.2.1 Polynomial lines and their directions. 



Let u(t) and v(t) be two distinct polynomial sites. The polynomial line l uv (t) is 
defined by 

1 U(t) x U(t)y 

1 V(t) x V(t)y 

1 x y 



The ruling coefficient rc(/(t)) of a polynomial f(t) S M[t] is the coefficient in its 
lowest degree term. The ruling sign rs(/(i)) is the sign of the ruling coefficient 
rc(/(t)). We define the direction tp uv of the polynomial line l uv (t) at t = as the 
argument of the point 



(cos (f) u 



, sin ( 



■) 



hm( , —. . — 

U0 V | d(t) x 



d{t\ 



r) 



where d(t) = (d(t) x , d(t) y ) = v(t) — u(t). Note that <f) uv is determined up to multiples 
of 2ir. We often use that value of 4> U v that lives in (— vr, tt]. 

Example 4.2 

Let u(t) = (t, t) and v(t) = (-t, t 2 ). Then l uv (t) = t 2 + 1 3 + tx - t 2 x - 2ty. Putting 
luv{t) equal to zero yields an ordinary line for every t ^ 0. As d(t) = v(t) — u(t) = 
(—2t, —t + t 2 ), we have lim t ^ = \- The direction of l uv {t) at t — is given 
by arctan(i) — 7r, however, as the ruling sign rs(d(i)a;) equals —1. 



Remark 4.3 

Note that this definition of <p uv indeed matches with the direction of the directed 
line that passes first through u(t) and then through v(t) for small positive t. For a 
definition that would match with small negative t, we should take into account the 
odness or evenness of the power of the lowest degree term of l uv (t). 



Example 4.4 

Let u and v be as in Example 4.2. For small negative t, the direction of the line 



^uv that passes first through u(t) and then through v(t) is close to arctan(i). If we 
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multiply both u(t) and v(t) by t, that is change u(t) into u(t) = (t 2 ,t 2 ) and v(t) 
into v(t) = (— t 2 ,f 3 ), then the direction of l uv for small negative t and for small 
positive t both are close to arctan(i) — 7r. 

4.2.2 Collinearity. 

Let u(t), v(t) and w(t) be three polynomial sites. As long as the sites do not coincide 
we can analyze whether they are collinear. Consider the determinant D(t) given by 



D(t) = D m (t) 



1 U{t) X U{t)y 
1 V(t) x V(t)y 

1 W(t) X w(t)y 



We call the three sites collinear if D(t) is the null polynomial. Otherwise let 
rs(Z?(i)) denote the ruling sign of D(t). If rs(D(t)) = 1, we say that w(t) is on the 
left of the polynomial line l U v{t) at t = 0, and if rs(D(t)) = —1, we say that w(t) 
is on the right. 

Example 4.5 

u{t) = (0,-t), v(t) = (-t,0), w(t) = (t,-2t). Then D(t) = 0, so u, v, and w are 
collinear at t = 0. 

Example 4.6 

Change w to w(t) = (t - t 2 ,-2t). Now D(t) = t 3 , so w is on the left of the 
polynomial line l uv {t) at t = 0. 

Example 4.7 

For small negative t, the site w(t) is on the right of l uv . If we change w(t) into 
w(t) = (t — t 3 , —2t), then w(t) is on the left of l uv (t) for small negative t while it is 
still on the right for small positive t. 



4.2.3 The center of a circle. 

Let u(t), v(t), and w(t) be polynomial sites such that u(t), v(t), and w(t) are not 
collinear at t = 0. Let D = D(t) be the determinant of above and let d = d(t), and 
e — e(t) be the determinants given by 



u (t) y u (~t)y 
v(tf x +v(t) 2 y v{t) y 

™(t)l+w(t) 2 y W(t)y 



■u(t) y u(t) x 

v{tf y v(t) x 
^2 



We define the circle center of u(t), v(t), and w(t) at t = as the point c given by 
the coordinates 

c T = — lim , c„ = — lim — — . 

We allow a circle center to be located at infinity. When leaving out the time 
dependency, these are of course the ordinary formulas describing a circle center. 
The circle C(u(t),v(t),w(t)) is oriented clockwise at t = iff w(t) is on the right 
of l uv (t) at t = 0. If w(t) is on the left of l uv (t) at t = 0, then C(u(t), v(t), w(t)) is 
oriented counterclockwise. 

If w(0) = v(0) = w(0) = (0, 0), then we say that the circle defined by u(t), v(t) and 
w(t) has positive radius at t = if and only if 



maximum ( \c x \ , \c y \ ) > 0. 



(4.1) 



4.2. Preliminaries. 
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Example 4.8 

Let again u(t) = (0, -t), v(t) = (-t, 0), and w(t) = (t - t 2 , -2t). Then a = t 3 and 
d = e = 4t 3 — 2t 4 + 1 5 . So c = (—2, —2) and we can conclude that u, v and w define 
a circle of positive radius at t = 0. As is on the left of l U v(t) at t = 0, the circle 
C(u(t) 7 v(t),w(t)) is oriented counterclockwise at i = 0. 

Example 4.9 

Let it(t) = (-i 3 ,2i), = (t 3 ,-2t), and = (-t 4 ,3i 2 ). Then C(u,v,w) is 

oriented clockwise at t = 0, while the circle center c of C{u,v,w) is situated at 
infinity. 

4.2.4 Cocircularity. 

Let u(t), v(t), and w(t) be distinct polynomial sites. Assume that at t = the 
sites u(t), u(t), and are not collinear and that the circle C(u(t),v(t),w(t)) is 
oriented clockwise. As long as the sites do not coincide, we can analyze if a fourth 
polynomial site q(t) is inside or outside of the circle defined by u(t), v(t), and w(t). 
Consider the polynomial I(t) £ M[t] defined by 



I(t) = 



u(t) x 


U(t)y 






1 


v(t) x 


V(t)y 


v(t)t - 




1 


w(t) x 


W(t)y 




Vw(t) 2 y 


1 


q(t)x 


q(t) y 


q{t)l - 




1 



We call the four sites cocircular at t = if I(t) is the null polynomial. Let 
otherwise rs(J(i)) denote the ruling sign of I(t). If rs(/(i)) = 1, we say that q(t) 
is outside of the polynomial circle c MVW (t) at t = 0, and if rs(D(t)) = —1, we say 
that q(t) is inside the circle. A set of polynomial sites S(t) is in general position 
at t = 0, iff no three sites are collinear at t = and no four sites are cocircular at 
t = 0. 

Example 4.10 

Let u(t), v(t), and w(t) be as in Example [4.8[ Let q(t) = (0,0). We check whether 
q(t) is inside or outside of the clockwise oriented circle C(u(t), w(t), v(t)) at t = 0. 
As I(t) = At 4 + 0(t 5 ), we conclude that q(t) is outside of C(u(t),w(t),v(t)) at t = 0. 

Remark 4.11 

For small negative t, the circle C(u(t),w(t),v(t)) is oriented counterclockwise. So 
the orientation swaps at t = 0. The site q(t) is inside the clockwise oriented circle 
C(u(t),v(t),w(t)) for small negative t. 

4.2.5 Ordering the sites. 

The lexicographic ordering for polynomial sites is as follows. Let u(t) and v(t) 
be two polynomial sites. First consider the polynomials u x and v x that give the 
x-component. If u x ^ v x , order u(t) and v(t) according to the coefficients of the 
term of lowest degree of u x and v x that are distinct. If u x — v x , then compare the 
coefficients of lowest degree that are distinct of u y and v y . We denote this ordering 
byO<. 

Example 4.12 

Let u(t) = (2t 2 ,t) and v(t) = (t,t). The first order coefficient of u x equals while 
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the first order coefficient of v x equals f . As this is the lowest order coefficient of u. 
and v x that is distinct, it follows that u(t) < v(t) with respect to < . 



4.3 The type of a set of polynomial sites. 



Let S(t) — {pi(t), . . . ,p n {t)} be a list of distinct polynomial sites that is in general 
position at t = 0. A clockwise polynomial circle c uvw (t) with u(t), v(t), and w(t) 
in S(t) is a Voronoi circle at t = iff there are no polynomial sites from S(t) 
inside c U vw(t) at t = 0. Represent a Voronoi circle c(t) by the ordered list of labels 
of the polynomial sites u(t), v(t), and w(t) that define c(t). The order of the labels 
corresponds to the c yclic, clockwise order of the defining sites on c(t) at t = 0, 
compare with Section 4.2.3. The set of all these lists is the type of S(t) at t = 0. 




Figure 4.1: The Delaunay triangulation of the sites of Example 4.13 at t = 0.1. 



Example 4.13 

Let S(t) be the set of polynomial sites given by 

qi = (2i, 2< 3 + t 4 ), 

q 2 = (-2i-i 3 , 2i), 

q 3 = (-2t 2 -2i 4 , -3t + 2i 2 -i 3 + f 4 ), 

94 = (i 2 - 2£ 3 + i 5 , -i 4 ). 



The relative position of the sites at t = 0.1 is shown in Figure 4.1. The set S(t) de- 
fines at t = the clockwise oriented circles 134, 142, 132 and 243, where we list only 
the labels. The polynomial site q\ is inside the circle 132 at t = 0, while the other 
circles are empty at t = 0. Therefore the type of S(t) is given by {142, 134, 243}. 

Remark 4.14 

Note that this definition of type again represents the situation for small positive t. 
This is because the type is defined in terms of the 'inside circle relation', compare 



Section 4.2.4 



4.3.1 Abstract Delaunay graph. 

The notion of type enables us to define an abstract Delaunay graph aD(S(t)) of 
S(t) at t = 0. In fact, we use Lemma 2.7 as a definition for aD(S(t)). The vertices 
of aD(S(t)) at t = are the labels occurring in the type. Two labels are connected 
by an edge if both occur in the same Voronoi circle at t = 0. Note that the closed 
paths of length three in aD(S(Q)) correspond with the Voronoi circles at t = 0. The 
multiplicity of an edge is the number of distinct Voronoi circles that contain both 
the vertices incident to the edge. 



Example 4.15 

Consider the type at t = of the set of polynomial sites from Example 4.13. It 



4.3. The type of a set of polynomial sites. 



43 



defines a abstract Delaunay graph aD(S(t)) with vertex set {1,2,3,4} and edge 
set {12, 13, 14, 23, 24, 34}. The multiplicities of the edges 12, 23, and 13 equal one, 
those of the other edges equal two. Compare this to the Delaunay triangulation 



shown in Figure 4.1 



Property 4.16 

Let p a (t) and pb(t) be two distinct sites of S(t). 

(i) If p a (t)pb(t) forms an edge of the convex hull CH(S(t)) for small enough 
positive t, then the multiplicity of the edge ab in aD(S(0)) equals one at t = 0. 

(ii) If p a (t)pb(t) forms an edge of the Delaunay triangulation D(S(t)) but not 
of CH(S(t)) for small enough positive t, then the multiplicity of the edge ab in 
aD(S(0)) equals two at t = 0. 



Proof. Combine the definition of type and Lemma 2.3 



4.3.2 Combinatorial convex hull. 

The convex hull of a set of distinct points in the plane can be defined as an inter- 



section of half-planes, see Section 2.1. But if a set of sites shrinks to one point, 
the convex hull of the sites also shrinks to one point. However, we can define a 
combinatorial convex hull of a set of polynomial sites at t = by making use of 
half-planes. 

Let S(t) be a set of distinct sites in general position. It is not necessary that all 
sites coincide in (0,0) at t = 0. Some i is a vertex of the combinatorial convex 
hull cCH(S(t)) of S(t) at t = if there exists a polynomial site Pj(t) ^ Pi(t) such 
that all sites Pk(t) with k ^ i,j are on the right of the polynomial line hj(t) at 
t = 0. We call such a line a bounding line. Two vertices i and j are connected by 
a directed edge, notation i — > j, iff hj(t) is a bounding line of cCH(S(t)) at t = 0. 

Property 4.17 

Let S(t) be a set of distinct polynomial sites S(t) in general position. Cq denotes 
the combinatorial convex hull cCH(S(t)) at t = 0. Let a, b € Cq and p c (t) £ S(t). 
(i) Cq is an abstract directed circuit graph. 

(ii) The edges of Co are exactly the edges of multiplicity 1 of the abstract Delaunay 
graph aD(S(t)) at t = 0. 

(Hi) The direction of an edge ab of Co equals a —> b iff the Voronoi circle contain- 
ing both a and b is oriented like abc. 

Proof. We prove the claims separately. 

(i) A circuit graph is a connected graph that is regular of degree two. For t 
positive, small enough, the boundary of CH(S{t)) is an oriented polygon where the 
vertices have the same labels as the vertices of cCH(S(t)) at t = 0. The polygon is 
oriented so that all points of S(t) not lying on an edge pi(t)pj(t) of CH(S(t)) are 
on the right of the directed line defined by pi (t) and pj (t) that passes first through 
Pi(t). This orientation and the co nnect edness are inherited when t vanishes. 



(ii) This follows from Property 4.16 



(iii) If abc is the clockwise orientation of the Voronoi circle defined by p a (t), Pb(t), 
and p c (t), then p c (t)(t) is on the right of l a b(t) at t = 0. ■ 

Example 4.18 

Consider the ordered set S(t) = {u(t), v(t), w(t)} of polynomial sites introduced in 
Example [Dj. As w(t) is on the left of l uv (t) at t = 0, the combinatorial convex hull 
cCH(S(t)) of S(t) at t = is given by 132. 
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Example 4.19 

Let S(t) = {qi(t),q 2 (t),q 3 (t),qi(t)} as in Example fOg. Then cCH{S{t)) at t = 
equals 132. 



4.4 The Voronoi diagram of a set of points and 
angles between the points. 

In this section we define the Voronoi diagram of a set of points and angles between 
the points. This enables us to introduce a Voronoi diagram for a set of polynomial 
sites. Moreover, we connect the notions of polynomial lines and polynomial circles 
in the context of Voronoi diagrams. 

4.4.1 Definition. 

Suppose S — {pi, . . . ,p n } is a set of n not-necessarily distinct points in the plane. 
To every pair of points (pi,pj) we add an angle otij such that the following rule 
holds: if two points pi and pj are distinct then there is a unique line though pi and 
Pj that makes some angle onj £ R/27rZ with the x-axis, directed from — oo to oo; if 
Pi and pj do coincide then we allow every value in R/27rZ for a,j. 

Let 7„ be a set of n points pi(j n ), ■ ■ ■ ,Pn(jn) an( i (<>) angles ai2, ■ ■ • , ai n ^x)n t na t 
obey aboves rule. Fix two points pi = Pi("f n ) and pj = Pji^n)- The bisection 
point b(pi,pj) is the point \{Pi + Pj)- If Pi 7^ Pj, the bisection point is just the 
middle of the line segment PiPj. If Pi — Pj then b(pi,pj) coincides with the double 
point pi — pj. The perpendicular bisector B(pi,pj) is the line through b(pi,pj) 
perpendicular to the angle otij — 0^(7™)- Let n be any non-zero vector, pointing in 
the direction aij. The Voronoi half-plane vh(pi,pj) is the half- plane defined by 

n-(x-b x ,y-b y ) < 0. (4.2) 

The Voronoi cell V(pi) is defined as 

V{Pi) = P| vh{pi,pj). 

j^i 

A point x is on the Voronoi edge e(pi,pj) iff it is on the intersection of the Voronoi 
cells V(pi) and V(pj), that is 

xeeipupj) ^> x G V(pi) n V(pj). 

The Voronoi diagram is the family of subsets of K 2 consisting of the Voronoi 
cells V(pi) and all of their intersections. The shape or boundary of the Voronoi 
diagram is the union of the boundaries of the Voronoi cells. 

Remark 4.20 

We will always use this definition in a geometric context which imposes restrictions 
on the angles atj . 

Remark 4.21 

We show how to get Equation 4.2. The line given by n • (x, y) = defines the line 
through the origin that is perpendicular to n. We have to translate this line over 
the bisection point in order to get the bisector B(pi,pj). As n • n > 0, the Voronoi 
half-plane vh(pi,pj) is the half-plane bounded by the line B(pi,pj) in the direction 
of the tail of n. 
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Remark 4.22 

Taking for n the vector (cos a,j , sin a,j ) of length 1 in the direction gives the 
inequality (y — b y ) sma^ < (b x — x) cosaijj for the Voronoi half-plane vh(pi,pj). 

4.4.2 The Voronoi diagram of a set of polynomial sites. 

Let S(t) = {pi(t), . . . ,Pn(t)} be a set of polynomial sites that is in general position 
at t = 0. The angle aij(0) for any two sites Pi(t) and Pj(t) is just the direction <pij 
as defined in Section [4.2[ We define the Voronoi diagram of the set of polynomial 
sites at t = as the Voronoi diagram of the points (pi(0), . . . ,p n (0)) and the angles 
«12, ■ ■ ■ , a(„_i)„. 

Remark 4.23 

At this stage, we do not define Voronoi vertices: if a Voronoi vertex is defined as an 
intersecti on of three Voronoi half-planes, it can become a complete line, compare 
Example [4.26| . 



Example 4.24 

Let S(t) be the set of polynomial sites, introduced in Example 4.13 . In Table 4.1 
we list the Voronoi half-planes, the directions of the bisectors that bound the half- 
planes, and the inequalities defining the half-planes. This results in the Voronoi 



diagram V(S(0)), depicted in Figure 4.2. The Voronoi regions V{q\),V{q2) and 



V(qz) all have positive area, while V{qi) = {(0,0)}. This shows that the shape of 
the Voronoi diagram of q±, q2, 93, and #4 equals the shape of the diagram of q±, q2, 
and (73, with q^ omitted. 



half-plane direction inequality 

vhi-2 7t — arctani 2x > y 

' 3 2 

vhi-3 arctan| — n —fx<y 

Vhl-i 7T x > 

vh.2-,3 — arctan| \x <y 

vh 2 -i -f x < y 

vh 3 . A f y < 

Table 4.1: The half- planes, the directions of the lines bounding the half- planes, and 
the inequalities of the half-planes defined by qi(t), q2(t), 53 (t) and 54 (t) at t = 0. 




Figure 4.2: The Voronoi diagram of the points qi(t), q2(t), qa(t) and 94(f) at t = 0. 



Remark 4.25 

Given the particular sites in Example 4.24, it is easy to check that the picture of 



the Voronoi diagram for small negative values of t is close to the image under the 
point reflection in (0, 0) of the picture in Figure 4.2. 
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Example 4.26 

Let S(t) be the set of polynomial sites of Example 4.6 and Example [T^, that is, 
u(t) = (0, —t), v(t) = (— t, 0), and w(t) = (t — t 2 ,—2t). Then vh uv is given by y < x, 
the half-plane vh UtW by y > x, and vh v ^ w by j/ > x as well. Therefore, V(u) is just 
the line y — x, while V(u) is given by y > x, and V(io) by y < x, all at t = 0. 



4.4.3 Polynomial bisector. 

Let u(t) and v(f) be two distinct polynomial sites. We have introduced the bisector 
of u and v at t = as the line passing through the bisection point b(u, v) and 
perpendicular to the direction <f> uv . We can derive directly an equation B t (u, v) for 
the bisector depending on t as follows. We construct two polynomial sites b(t) and 
a(t) that lie on the ordinary bisector for every t. The point b(t) is just the bisection 
point \{u(t) +v(t)). The other point a(t) is the image of u(t) under the | rotation 
around b(t), that is, 

a(t) = ( I -^.(u-bf + b. 

Finally, put B t (u,v) := lba(~t)- We call B t (u,v) the polynomial bisector of u(t) 
and v(t). 

Lemma 4.27 

Let u(t) and v(t) be two polynomial sites. Let B(u,v) be the perpendicular bisector 
of u(t) and v(t) at t = 0. Then 

B(u,v) = liraB t (u,v). 

Proof. By construction, the direction of B t (u,v) at t — is perpendicular to cf> uv . 
As B(u,v) and lim t |o B t (u, v) both pass through the bisection point b(0), the two 
lines coincide. H 

Example 4.28 

Let S(t) — {u(t),v(t),w(t)} be as in Example 4.26. The polynomial bisector 
B t (u, v) is given by 

B t (u,v) = {(x,y) el 2 \y = x}, 

while 

B t (u, w) = {(x, y) € R 2 | y = -2t + t 2 - -t 3 + x - fx}. 



4.4.4 Half-planes and Voronoi circles. 

Let S (t) be a set of polynomial sites that is in general position at t = 0. In Section 
4.4. l| we have introduced the Voronoi diagram V(S(t)) of S(t) by means of half- 
planes, while in Section [T^ we have defined the type of S(t) at t = in terms 
of Voronoi circles. The following lemma gives a connection between polynomial 
bisectors and circle centers. 

Lemma 4.29 

Suppose thatpi(t), Pj{t), and Pk(t) are distinct non-collinear polynomial sites. De- 
note by c the circle center ofpi(t), Pj(t), and pk(t) at t = 0. Let Bt(pi,pj) be the 
polynomial bisector of pi(t) andpj(t). Then 



c = \ua(B t (pi,Pj) n B t (pi,p k )). 



4.5. The positive area cells for one cluster. 
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Proof. For t positive and small enough, pi (t) , pj (t) , and pk (t) are distinct points 
in general position. For such points Lemma |2.4| states that the intersection of the 
bisectors of the points is the center of the circle defined by the points. But then it 
holds for t = as well. ■ 

Now we are able to define a Voronoi vertex so that a vertex is exactly one point in 



the plane, compare with Remark 4.23 . Let Pi(t), Pj(t), and Pk(t) be distinct non- 
collinear polynomial sites. Let c = Cyk(Q) be the clockwise oriented polynomial 
circle through Pi(t), Pj(t), and pk(t) at t = 0. Let x be the circle center of Pi(t), 
Pj(t), and Pk(t) at t = 0. Then a; is a Voronoi vertex v(pi,Pj,pk) at t = if c is 
a Voronoi circle at t = 0. 

Let u(t), v(t), and w(t) be distinct polynomial sites that coincide in (0, 0) at t = 0. 
Then the circle center of u(t), v(t), and w(t) is not equal to (0,0) only for special 
configurations. 

Lemma 4.30 

Let u(t), v{t), and w{t) be distinct polynomial sites that coincide in (0,0) at t = 0. 
If ' <Puv ^ <Puw mod 7r or <p uv ^ <p vw mod n, then the circle center c ofu(t), v{t), 
and w(t) at t = equals (0,0). 

Proof. We may assume without loss of generality that 4> uv ^ 4> uw mod it. Lemma 
L2S| states that c = lim t ^o(Bt(u, v) n B t (u,w)). From Lemma 4.27 it follows that 



c € B(u, v) n B(u, w) at t — 0. But then c is the unique intersection point (0, 0) of 
two non-parallel lines B(u,v) and B(u,w). ■ 

Example 4.31 



Let u(t) = (0,-t), v(t) = (-t,0), and w(t) = (t - t 2 ,-2t). In Example 4.8 we 



have shown that u, v, and w define a circle with center c = (—2,-2) at t = 0. 



We did so by evaluating the determinants a, d, and e introduced in Section 4.2.3. 
Alternatively, we can compute c(t) as an intersection of two bisectors, say B t {u, v) 
and B t (u,w). This gives: 

c(t) = B t {u,v)r\B t (u,w) 

= (x,x) n {x, -2t + 1 2 - i* 3 + x - tx) 

= i(-4 + 2t-£ 2 ,-4 + 2£-< 2 ). 
So c = c(0) = (—2, —2) as ex pecte d. We have already computed the equations for 



the two bisectors in Example 4.28 



4.5 The positive area cells for one cluster. 

In this section we assume we have a set S of polynomial sites that all coincide in 
the origin at t = 0. We show that in general some sites can be omitted without 
changing the shape of the Voronoi diagram of S. We characterize those sites that 
determine this shape and present an efficient algorithm to compute the shape of the 
Voronoi diagram. 

4.5.1 The direction hull and positive area cells. 

A zero cluster at t = is a set S(t) of polynomial sites in general position such 
that Pi(0) — (0,0) for every site Pi(t) in S(t). We want to know the shape of the 
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Voronoi diagram of a zero cluster. For this purpose it is enough to determine the 
boundary of the Voronoi cells of positive area at t = 0. 

Assume that S(t) is a zero cluster at t = 0. Fix some site v(t) in the combinatorial 
convex hull cCH(S(t)) at t = 0. There is one incoming edge and one outgoing edge 
in cCH(S(t)) at v(t) at t = 0. Suppose the incoming edge comes from u(t), while 
the outgoing edges goes to w(t). Site v(t) is called a corner site of S(t) at t = 
if the direction of the incoming edge is distinct from the direction of the outgoing 
edge, that is, if 4> uv ^ 4> vw . The corner sites of S(t) at t — are by definition the 
vertices of the direction hull DH(S(t)) at t = 0. Two vertices u(t) and are 
connected by a directed edge if there are no corner sites on the path from u(t) 
to v(t) in cCH(S(t)). 

Example 4.32 

Let S(t) = {q±(t), 54(i)} be as introduced in Example 4.13 



4.18 we know that cCH{S{t)) equals 132. Compare Table 



4.1 



From Example 
The directions 

<pi3 = arctan|, (^32 = 7r — arctan| and 4>2i = — arctang are all distinct. It follows 
that DH(S(t)) = cCH(S(t)) = 132 at t = 0. 

Example 4.33 

Let S(t) = {gi(t),? 2 (i),g 3 (*)} with 9l (t) = (-i,0), q 2 (t) = (0,t 2 ), and gj,(t) = 
(t, -t 2 ). Then cCH{S{t)) = 123, while DH(S(t)) = 13, as 12 = ^23 = 0. 

Lemma 4.34 

Let S(t) — {pi(t), . . . ,p n (t)} be a zero cluster at t = 0. Then any bisector B(pi,pj) 
passes through (0, 0) at t = 0. 

Proof. The bisector B{pi,pj) passes through the bisection point b(pi,pj) = (0,0) 
at t = 0. ■ 

Lemma 4.35 

Let S(t) = {pi(t), . . . ,p n {t)} be a zero cluster at t = 0. Then V(pi(0)) has positive 
area if and only if Pi G DH(S(0)). 

Proof. Assume that p { e DH(S(0)). As pt S cCH(S{0)), it follows that p, G 
CH(S(t)) for t > 0, small enough. That means that there exists to such that for all 
< t < to the Voronoi cell V(pi(t)) is unbounded. Moreover, pi is incident with two 
edges e(pi-i,pi) and e(pi,pi+i), say, that are unbounded for all < t < to. But this 
implies that these edges are unbounded at t = as well. The sites Pi-i and Pi+\ arc 
the direct predecessor and direct successor of Pi on cCH(S(0)). As pi G DH(S(0 j), 
the directions 4>i-i-% and <f>i-i+\ are distinct at t — 0. But this implies that V(pi(0)) 
has two unbounded edges of distinct direction on its boundary and by the convexity 
of V(pi(0)) this implies that V(pi(0)) has positive area. 

For the other direction, assume that V(j>i(0)) has positive area at t = 0. It follows 



from Lemma 4.34 that V(pi(0)) is unbounded, and therefore has two unbounded 
edges on its boundary of distinct direction. This implies that pi G DH(S(t)) at 

t = o. m 

The Voronoi diagram V(S(0)) of a zero cluster at t — looks as follows: 



1. If pi(t) & DH(S(t)) at t = 0, then area(V(S(£))) = by Lemma $M . 

2. If pi(t) G DH{S(t)) and if p,_i(t) and are its direct predecessor and 
its direct successor at t = in DH(S(t)), then 



V(pi(Q)) = vh (pi,Pi-i) r\vh (pi,p l+1 ). 



4.5. The positive area cells for one cluster. 
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Corollary 4.36 

The shape of the Voronoi diagram SV(S(0)) of a zero cluster S(t) equals the shape 
of the Voronoi diagram of DH(S(0)). 

We conclude that for computing the shape of the Voronoi diagram 6V(S(Q)) of a 
zero cluster, it is enough to consider sites on the direction hull DH(S(0)) C S(t) 
only. 



4.5.2 Determining the direction hull at t = 0. 

In Section |4.2| we have defined the combinatorial convex hull cCH(S(t)) of a set 
of polynomial sites by imitating the characterization of the ordinary convex hull as 
an intersection of half-planes. Exploiting this similarity we show how to compute 
cCH(S (t)) by an adapted version of the convex hull algorithm as presented in 
BKOSfl , page 6. 



The original algorithm is an incremental algorithm: first it sorts the n points in the 
input and next it handles the points one by one. The original algorithm sorts the 
input points lexicographically. This results in a sequence of points that is ordered 
from left to right and then from bottom to top. The algorithm determines the 
convex hull in two steps. In the first step the upper hull is determined. The upper 
hull is basically that part of the convex hull that is running from the leftmost point 
to the rightmost point of the sorted points. It constructs the upper hull by adding 
one point at a time and checking if the sequence of points that will be the upper 
hull in the end keeps going right. The lower hull is defined and constructed in a 
similar way. 

Our ordering < mimics the lexicographic ordering. Just think of substituting a 
very small positive value of t in the set of polynomial sites a nd s orting the resulting 



points lexicographically. An example is displayed in Figure AA where the position 
of some polynomial sites at t — 0.59 is indicated by the labels of the sites. In a 
similar fashion we define right turns for polynomial sites: we say that three sites 
u(t), v(t), and w(t) ordered by 0< make a right turn at t = 0, if w(t) is on the 
right of l uv (t) at t = 0. 

Algorithm 4.37 

Combinatorial convex hull for polynomial sites. 

Input: set S(t) of n polynomial sites in general position. 

Output: list L containing the vertices of cCH(S(0)) in clockwise order. 

1: Sort the sites by < , resulting in a sequence pi(t), . . . ,p n {t). 

2: Create list L up p Gr = {pi 1 p 2 \. 

3: for every i S {3, . . . , n} do 

4: Append p t to L U p pcr . 

5: while L U ppcr contains more than two sites and the last three sites in L U pper 

do not make a right turn, do 
6: Delete the middle of the last three sites from L U pper- 

7: Construct Li ovrel in a similar way and append the result to i upP er- Call the 

resulting list L;. 
8: return L. 

The running time of the algorithm is 0(n log n), cf. p3KOS| . Let S(t) be a zero 
cluster again. Given its combinatorial convex hull cC'H(S(t)) at t = 0, we compute 
in 0(n) time the direction hull DH(S(t)) at t = 0. 
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Algorithm 4.38 

Direction hull. 

Input: combinatorial convex hull cCH(S(0)) of a zero cluster S(t). 
Output: direction hull DH(S(t)) at t = 0. 

l: Compute the direction <p uv for every edge uv in cCH(S(t)) at t = 0. 

2: if two consecutive edges uv and vw have the same direction, that is, (j> uv = tf>. 
then 

3: Delete w and its incident edges. 
4: Create a new edge ttw;. 



4.6 Degenerate Voronoi cells and edges. 

Let S(t) be a zero cluster at t = 0. We have seen before that only sites in the 
direction hull DH(S(0)) generate a Voronoi cell of positive area. It is possible 
however that for a site pi € S(t) \ DH(S(0)) there exists an edge e(pi,pj) that does 
not collapse on (0,0) at t = 0. In this section we consider sites Piit) £ S(t) such 
that there exists, at t = 0, an edge e(j?i,Pj) with x £ e(pi,pj), x ^ (0,0). We 
first give an explicit example of an edge e(pi,pj) where even (0, 0) £ e(Pi,Pj), while 
Pi ^ DH(S(0)). Next we present an algorithm that determines all edges e(pi,pj) 
that contain some x £ R 2 with a; ^ (0, 0) for a zero cluster S(t) at t = 0. 




Figure 4.3: Edge of zero length outside (0, 0). 



Example 4.39 

Let S(t) = (p(t),q(t),r(t),s(t)), where p(t) = (-2t,0), q{t) = {-t,-{t 2 ), r(t) = 
(0, 0), and s(t) = (2t, 2t 2 ). The direction hull at t = consists of the two sites p{t) 
and q(t), while pg = 0. Therefore the shape of the Voronoi diagram JV^S^O)) is 
just a vertical line through the origin. There are two Voronoi cells of positive area 
at t = 0: V(p), whose cell is the left half-plane, and V(s), whose cell is the right 
half-plane. 

The type of V (S(t )) however is prq,psr. A schematic picture of the situation is 
given in Figure |4~3| . From the type it follows that combinatorially, V(S(t)) consists 
of two vertices or circle centers, prq and psr that are connected by an edge e(p, r). 
Besides, there are unbounded edges e(p,q) an e(q,r), incident to the circle center 
prq and edges e(r, s) and e(p, s), incident to psr. Both circle centers are situated at 
the point (0,2). This implies that the length of e(p,r) at t = equals zero, while 
the four other edges involved, have infinite length. 

Let S(t) be a zero cluster at t = 0. An edge e(pi,pj) with Pi,Pj £ S(t) is an outside 
edge at t = if there exists x £ e(Pi,Pj) such that x 7^ (0, 0). Here we allow x to 



4. 7. Example: 20 polynomial sites. 
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be any point of the form x — (a, b), x — (±oo,6), x = (a,±oo) or x — (±oo,±oo), 
where a, b <E M. Allowing these possibilities is motivated by Example |4.S| , which 
demonstrates that circle centers defined by polynomial sites can be at infinity. 



Algorithm 4.40 determines all outside edges of a set of polynomial sites in general 



position at t = 0. 
Algorithm 4.40 

Determining outside edges and their positions. 
Input: a zero cluster S(t). 

Output: all outside edges with their endpoints at t = 0. 
1: Determine the type T of S(t) at t = 0. 

2: Any edge uv of multiplicity 1 in the abstract Delaunay graph aD(S(0)) cor- 
responds to an unbounded edge e(u,v) in the Voronoi diagram V(S(t)). This 
edge e(u, v) starts at the circle center of the unique circle listed in T that has 
both u and v on its boundary. The direction of e(u, v) is perpendicular to <j> uv . 
Edge e(u,v) is oriented so that it makes an angle of -| with <fi uv ; 
3: Determine those circles in T that have positive radius; 
4: for a and b two labels occurring in some circle of positive radius do 
5: (We show how to 'draw' the edge ab). 

Check if we have not drawn ab before in step 2 as an unbounded edge, 
if the combination of the labels a and b occurs twice in the list of the circles 
of positive radius then 
8: there are two circle centers distinct from (0,0). We connect the two circle 

centers by an edge ab; 
9: else if the combination of the labels occurs only once then 
10: we connect the corresponding circle center c a & to the origin. 

Determining the shape of the Voronoi diagram is relatively easy, compare Section 



4.5. We have seen in this section, however, that the combinatorial structure on the 



boundary itself can get rather complicated. 



4.7 Example: 20 polynomial sites. 

In this section we analyze a slightly bigger example, illustrating the techniques and 
theory developed in the former sections. We consider a set S(t) of 20 polynomial 
sites. We determine the shape of the Voronoi diagram V(S(t)) at t = Q. Next 
we focus on the outside edges of V(S(t)) at t — 0. We visualize the polynomial 
sites and compare the outcome of the computation to the visualization. First we 
introduce the polynomial sites. S(t) is a set of 20 distinct sites (pi(t), . . . ,P2o(t)) 
such that pi(t) = (0,0) for i = 1, . . . , 20: 



Pi 


(0,i 5 ), 


Pn 


Pi 


(0,2t + £ 4 + 2i 5 ), 


Pu 


P3 


(2t,2t 3 + t 4 ), 
(t 2 ,2t + 2t 2 -t 3 ), 


Pl3 


Pa 


Pli 


P5 


(-2t\-t*), 


Pl5 


Pe 


(3t 5 ,t 2 + 2t 3 + 3t 5 ) 7 


Pl6 


P7 


(2i + 3t 2 -2t 3 ,2t-2t 2 -t 3 ), 


Pn 


P8 


(-2t-t 3 ,2t), 


Pis 


P9 


(-2t + t 3 ,-2t-t 3 ), 


Pl9 


Pw 


{t + 2t 2 + 2t 3 ,-t + 3t 3 ), 


P20 



-2t 2 -2t 4 ,-3t + 2t 2 -t 3 + t 4 ), 

-2t-t 4 ,"it 2 + 3t 3 -t 5 ), 

-2t 2 -t 4 ,-2t 2 + M 3 ), 

-t + it 2 - 2t 3 + 3fS t + 'St 2 - 3t 5 ), 

-t-2t 5 ,-3t 3 + 2t 5 ), 

t 2 - 2£ 3 + t 5 , -t 4 ), 

M 3 + 2^,-f 4 ), 

3t 3 + 3t 4 + 2i 5 ,2t-2^ 3 -t 4 ), 

t 2 +t 3 + 3i 5 ,0), 

M 2 + 2t 3 + 3i 5 , t 2 + 3i 3 + 2t A ). 
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We visualize S(t) by considering any Pi(t) £ S(t) as a plane curve, represented by a 
parametrization (pi, x {t),Pi,y(t))- I n Figure 4.4 these twenty curves are plotted for 
t £ [0, 1]. The position of a curve at f = 0.59 is marked by its label. 




\ 

Figure 4.4: The polynomial sites for t £ [0, 1]. 



4.7.1 The positive area diagram. 

Recall that the type of S(t) is defined a the set of all Voronoi circles at t = 0. 
We have computed it using brute force. That is, for every (™) unordered triples of 
polynomial sites pi (t) , pj (t) and pk (t) we first determine the orientation at t = of 
the circle C passing though Pi(t), Pj(t) and Pfc(t). Next we check if any of the n — 3 
remaining sites of S(t) is contained in C. If not, we conclude that C is a Voronoi 
circle and add the labels i, j and k ordered with respect to the orientation of C to 
the type in the making. It turns out that the type of S(t) is given by: 

((1, 5, 6), (1, 17, 5), (1, 6, 17), (2, 8, 4), (2, 4, 18), (2, 14, 8), (2, 18, 
14), (3, 20, 7), (3, 11, 10), (3, 10, 20), (4, 7, 18), (5, 13, 6), (5, 17, 13), 
(6, 13, 15), (6, 15, 14), (6, 14, 20), (6, 16, 17), (6, 19, 16), (6, 20, 19), 
(7, 20, 18), (8, 12, 9), (8, 14, 12), (9, 13, 11), (9, 12, 15), (9, 15, 13), 
(10, 11, 13), (10, 13, 20), (12, 14, 15), (13, 17, 16), (13, 16, 19), (13, 19, 
20), (14, 18, 20)). 

Counting multiplicities of the edges shows that all edges have multiplicity 2 except 
for the edges 

(4, 8), (3, 7), (3, 11), (4, 7), (8, 9), (9, 11). 

As site P7 is on the right of i^g at i = the combinatorial convex hull cCH(S(0)) 
is given by (8, 4, 7, 3, 11, 9). The directions of the edges of cCH(S(0)) are </>s ; 4 = 0, 
</> 4 ;7 = 0, </> 7 ;3 = 03;ii = -arctan(|) - |, <f> n . 9 = f + arctan 2, and 9;8 = |. 
The directions </>s ; 4 and 04 ; 7 are equal, so site p4 is not a corner site. We conclude 
that the direction hull DH(S(0)) is given by (8, 7, 3, 11, 9). This outcome can almost 



be 'guessed' from Figure 4.4 by considering the convex hull of the points pi(t) to 
pw(t) for small t. Knowing the direction hull and the directions of the edges of the 
direction hull, we know in fact the shape of the Voronoi diagram 6V(S(0)). It is 



depicted on the left in Figure 4.5 



4.8. The shape of a Voronoi diagram at t = 0. 
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4.7.2 Combinatorics of the edges outside (0,0). 

The only Voronoi circle that has a circle center outside of (0, 0) at t = is the 
clockwise oriented circle (8,12,9). Its center is at infinity. This means that the 
unbounded edge e(8, 9) starts at the Voronoi vertex at infinity, while two edges 
e(8, 12) and e(9, 12) run between vertices situated in (0, 0) and this vertex at infinity. 
We have seen before that p^ G cCH(S(0)) with p s its direct predecessor and p-j its 
direct successor in cCH(S(0)). As there are no further circle centers situated outside 
the origin, there are two unbounded edges e(4, 8) and e(4, 7) that, starting in (0, 0), 
both run upwards, but never meet. A schematic picture of the combinatorics of the 
edges outside (0, 0) is given on the right in Figure |4.5| . 

8 7 8 7 



8 
9 



11 




Figure 4.5: The positive area diagram at t = on the left. The combinatorics 
outside (0,0) included on the right. 



4.8 The shape of a Voronoi diagram at t = 0. 



In Section 4.5 we have given a method that determines the shape of the Voronoi 
diagram of a zero cluster at t = 0. That is, we assumed that all sites coincide in 
the origin at t = 0. In this section we apply this method in order to determine 
the shape of a Voronoi diagram at t = for an arbitrary set of polynomial sites in 
general position. 

Let S(t) be a set of polynomial sites of size n that is in general position at t = 0. 
The cluster locations 1(0) = {Zi(0), . . . , l m (0)} of S(t) at t — are the distinct 
positions of the sites in S(t) at t = 0. An lj-cluster, for lj G I 2 is a collection 
of polynomial sites that coincides at lj at t = 0. Write S(t) = UJL^j^), where 
Sj(t) C S(t) is the set of sites in S(t) such that Pj t (0) = lj(0) for all p ji G Sj(t). 
That is, we divide S(t) with respect to the ^-clusters at t = 0. 

Lemma 4.41 

The shape of the Voronoi diagram SV(S(t)) at t = of S(t) is given by 
SV(S(Q)) = |J (F(^(0))n5F(S,(0))). 



Proof. First of all, suppose that |<Sj(0)| = 1, for all j = 1, ...,m. In this case 
the lemma just states that the Voronoi diagram of distinct points coincides with 
the subdivision of the plane into Voronoi cells, see Section [2.2[ Next consider an 
arbitrary point x in the interior of, say, V(lj)(0). Then x is, by definition, closer to 



lj than to any other location Ik for fc ^ j. The lemma follows from Corollary 4.36 
adapted for an /j-clustcr Sj(t). 
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Example 4.42 

Suppose we are given the following sites: 

S(t) = ((-6,4),(-2,-6),(-l 1 -2-3t),(0 ) -3t),(3,5),(6,-3), 

(-1 - 3t, -2 + i), {-2t, -2t), (-2t, 2t), (2t, 0), (2t, 2t), (-1 + 2t, -2)). 

The cluster locations at t = are given by 

1(0) = ((-6, 4), (-2, -6), (-1,-2), (0,0), (3, 5), (6, -3)). 

The V oron oi diagram V(l(Q)) of the cluster locations is presented on the left in 
Figure 4.6. There are two clusters consisting of more than one point. The shape of 
the Voronoi diagram of the zero cluster, given by 



Si(t) 



((0, -3t), (-2t, -2t), (-2t, 2t), (2t, 0), (2t, 2t)), 



is presented in the top middle of Figure |4.6| , while the shape of the Voronoi diagram 
of the (— 1, — 2)-cluster 



s 2 (t) 



((-1, -2 - 3t), (-1 - 3t, -2 + t), (-1 + 2t, -2)), 



is shown in the bottom middle. According to Lemma 4.41, we can plug these two 
cluster diagrams in the cluster locations diagram in order to get the shape of the 
Voronoi diagram of S(t) at t = 0. This final diagram is depicted on the right in 
Figure 4.£. 



4.9 Generalizations and conclusion. 
4.9.1 Dropping general position. 

Throughout this chapter we have assumed that sets of polynomial sites are in gen- 
eral position. General position is not needed for the methods presented to work. 
Dropping general position just means that a lot of extra cases have to be checked, 
which distracts from the main line. Note that the notion of type as introduced in 



Section 2.3 is defined for point sets that are not in general position as well. As we 



can also compute convex hulls for point sets that are not in general position, there 



4.9. Generalizations and conclusion. 
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are no big obstructions for extending the methods to sets of polynomial sites in 
arbitrary position. 

Recall from Section 4.2 that general position for polynomial sites is defined in terms 
of cocircularity and collinearity polynomials: if some set S(t) is not in general posi- 
tion at t = 0, then there is at least one cocircularity polynomial I(t) or collinearity 
polynomial D(t) that is equal to the null polynomial. This implies that S(t) is not 
in general position for any value of t. 



4.9.2 Generalization to higher order Voronoi diagrams. 

In the /c-th order Voronoi diagram, the plane is partitioned according to the k closest 



sites, see Chapter H. Algorithm 3.5 computes the fe-th order Voronoi diagram of a 



set S of distinct points in general position, using 

• circles through three points of S, 

• points from S inside these circles, 

• directions of lines between points of S. 

We have shown in this chapter that these concepts can be generalized to sets S(t) 



of polynomial sites. Therefore, Algorithm 3.5 can also be generalized to sets of 
polynomial sites in general position. This will give us a combinatorial k-th order 
Voronoi diagram. 



Figure 4.7: No plugging of second order diagrams. 



It is also clear that we cannot directly generalize Lemma 4.41 for obtaining the 
shape of a fc-th order Voronoi diagrams. That is because the two steps approach, 
dealing with cluster locations and clusters apart, is not allowed anymore for k > 1. 

Example 4.43 



Let S(t) be as in Example 4.7. In Figure 4.7 we show on the left the second order 
Voronoi diagrams of the cluster locations Z(0) that we have found in Example 4.42. 
On the right, the second order Voronoi diagram for S(t), with t = 0.001 is shown. 



4.9.3 Polynomial arithmetic. 

So far, we did not pay any attention to how to compute efficiently with polynomial 
sites: higher order terms can often be skipped without changing the type or the 
position of the Voronoi vertices and edges. But at this point we have no rule 
available that tells us on forehand when it is allowed to skip higher order terms. 
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Example 4.44 

Let S(t) be as introduced in Example 4.13. Suppose that we change in Examples 



4.13 and 4.24 from q t = (t 2 - 2t 6 + t b , -t 4 ) into q 4 = (t 2 , -t 4 ). Then the type of 



S(t) at t = and the shape of the Voronoi diagram of S(t) at t — remains the 
same. 



This is one reason for considering certain properties of sites or point sets, rather 
than the sites themselves. We will do so in the following chapters. 



4.9.4 Conclusion 

The polynomial sites model enables us to introduce limit pictures of Voronoi dia- 
grams by means of the shape of a Voronoi diagram. Moreover, we have seen that this 
model can be used to extend the notion of 'type' to sets of points that sometimes 
coincide. The type gives us the complete combinatorics of a limit situation. 

The methods presented in this chapter however do only visualize the shapes of those 
cells that have positive area at t = 0. This demonstrates that not all cells in the 
limit Voronoi diagram at t — are treated in an equal way. The ingredient that is 
missing is the relation between shape and scale. Informally, think of dividing the 
polynomials defining the polynomial sites by powers of t until Voronoi cells that 
had zero area at t = before dividing get positive area. That is, we want to zoom 
in at a cluster in order to find out the shape of a Voronoi cell of some arbitrary 
site pi (t) . In order to solve this problem we will define clickable Voronoi diagrams 
in a following chapter. We do so by exploiting properties of point sets, like angles 
between two points or hooks between three points. 



Chapter 5 

Voronoi diagrams and angle 
compactifications. 



Given a configuration c of n distinct points, one can determine the set a n of angles 
of lines through any two points in c. We distinguish a n , the set of angles up to 
2ir, and a^, the angles up to ir. We analyze in both cases what sets of angles are 
possible. We show in which cases it is possible to reconstruct the Voronoi diagram 
V(c), knowing the angles only. We compactify the configuration space of distinct 
points by taking the closure of the graph of the map that associates the angles to a 
configuration. We present a variety T n as an algebraic alternative. We analyze the 
connection between boundary points of the compactification and singularities of T n 
for small n and give geometric interpretations. 



5.1 CDA n : space of angles on n points. 



In this section we introduce several spaces that will be important to us later on. 
We recall the notion of configuration space of n distinct labeled points in the plane. 
Both an intr oduc tion on and applications (in robotics!) of configuration spaces can 
be found in [AG|. Next we define two spaces by considering, for n distinct points 
in the plane, all angles between pairs of points. Here we distinguish angles mod 7r 
and angles mod 2it. 



5.1.1 The angle of two points. 

Given two distinct points Pi and pj , we determine the angle that the line that passes 
through pi and Pj makes with the positive s-axis. We distinguish the directed and 
undirected line. 



Definition 5.1 

See also Figure 5.1. 

(i) For any two distinct points pi and pj in the plane define the directed angle 
ctij S R/27rZ as the argument of the point 

(oosou.sinatf) = 

where d = (d x , d y ) ■= Pj — Pi- Note that ctij is determined up to multiples of 2-7T . 
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(ii) The undirected angle otij is defined as = etji = ciij mod tt. We often 
choose aij G (—§>§]• 

(Hi) If{p 3 -Pi)x ^ 0, then the slope a i3 is defined as a io = ^Ip*s" - If(Pj~Pi)x = 
0, t/iera ay =oo. -ft holds that S (— oo, oo]. 




Example 5.2 

Let pi = (1, 1) and pj = (—1, 0). Then d = pj — p,; = (—2, —1), so ^ = |, but the 
angle ay = arctan(i) + 7r as d x < 0. 

Remark 5.3 

If (pj — ^ 0, then ay = tana-TJ and 51J = arctana^. 



5.1.2 Configuration spaces of n distinct points. 
Definition 5.4 

The configuration space CONF n (M. 2 ) of n points in M. 2 is the set 

CONF n (R 2 ) {(p x , . . .,p n ) e (R 2 ) n \ Pl ? Pj if i £ j}. 

Two elements p,q € COJVF n (R 2 ) are ^-equivalent if they only differ by a scal- 
ing, s, combined with a translation, t. The reduced configuration space is the 
quotient space conf„ = CONF n / ^{s.t}- 

conf„ is a smooth manifold. We determine its dimension. A standard represen- 
tative for a class in conf„ is constructed as follows. Translate the configuration 
in such a way that p\ = (0,0). Scale the configuration in such a way that p2 is at 
distance one from p\. This shows that we need three parameters less than In to 
describe an element in conf„. 



Property 5.5 

dimension(conf n ) — 2n — 3. 

Remark 5.6 

Note that conf„ is not compact: take n 
but not bounded. 



= 3, then the ratio , Pl is continuous 

IP2-P1 



Remark 5.7 

Let c = {pi, . . . ,p n } be the standard representative of a class [c] in conf„ and let 
c = {pi, . . . ,p n } be an arbitrary element of [c). Then Pi —\ P2 — Pi \ Pi + Pi- If 
we write T CQn f for the set of standard representatives of classes in conf„, then 
CONF n = {cK> + M 2 | c € T conf J. 



5.2. From angles back to point configurations. 
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5.1.3 Compactification of the graph of the angle map. 

For an element c £ CONF n , write down for every pair of points Pi,Pj £ c, with 
i 7^ j, the angle ctfjj £ R/2ttZ or cEy 6 M./nZ. This gives (2) angles, one for every 
unordered pair of labels i and j. 

Definition 5.8 

DA n := (K/2ttZ)(2) is the space of directed angles; UA n :— the space of 

undirected angles. The directed angle map ipDA n is the map 

ipDAn : CONF n -> DA n , 

(pi,...,p„) l-> (atij) l<i<j<n- 

The undirected angle map ipUA n is defined in a similar way. 
Remark 5.9 

As a direct product of circles DA n and UA n are smooth. 
Remark 5.10 

ipDA n and "0ua„ are well-defined on conf n : both mapping are constant on classes 
of conf„ . 

Definition 5.11 

The graph of ip DAn is the set {(c,ip DAn (c)) £ (R 2 )" x DA n \ c e CONF n }. 
The compactification CDA n of the graph of the angle map, is the closure of 
graph(i/jrjA n ) in (K 2 )™ X DA n . The compactification CUA n is defined in a similar 
way. 

Remark 5.12 

CDA n is in fact not a compactification as it is not compact: CONF n is not bounded, 
compare Remark \j.6\ The projection map from CDA n to (R 2 )™ is proper however. 
Recall that a map is proper if it is continuous and if the preimage of every compact 
set is compact again. 

5.2 From angles back to point configurations. 

In this section we consider the mappings ^DA™ : conf„ — > DA n and V'UAn : conf n — > 
UA n . We describe the fibers of both mappings by showing which configurations c 
can be reconstructed from knowing the angles ipDA n (c) or V'UAn (c) only. 

5.2.1 Distinct points and angles in M./2nZ. 

A configuration c € CONF n is called collinear iff all points in c are collinear. A 
class of configurations [c] £ conf„ is collinear if the class elements are collinear. 
Define CLCONF n := {c 6 CONF n \ c collinear}, and similarly clconf„. 

Lemma 5.13 

The map ipDA„ ■ conf„\clconf„ — > DA n is injective. 

Proof. We have to show that for [c], [d] € conf„\clconf„ with [c] ^ [d] it holds 
that ipDA„([c\) ipV}A n ([d])- Let c and d be the standard representative of [c], resp. 
[d] . If [c] 7^ [d] then also c 7^ d. Let pi (c) be the i-th point of the configuration c and 
let ctij(c) £ R/27rZ be the angle between Pi{c) and Pj{c). Note that for standard 
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representatives p\(c) = pi(d) = (0,0). We show that in any case there exists some 
labels u and v such that a uv (c) 7^ a uv (d). 

Let i G {2, . . . , n} be minimal such that Pi{c) 7^ Pi(d). If i = 2, then pz(c) ^ P2(d), 
therefore 0.12(c) 7^ 0.12(d). Assume that i > 2. In this case P2(c) — P2(d) = 
(cosai2,sinai2) for some an G (— 7r,7r]. Let m be the line through p\ and Pi(c). 
If Pi(d) g 1 m, then a>n(c) 7^ au(d). Assume that Pi(d) € m. If P2 ^ to, then 

U2i{c) ^ Qt 2 i(d). 

We are left with the case that pi,p2,Pi(c) and Pi(d) arc all collincar on line m. As 
c,d£ CLCONF n , there exists j such that pj(c) ^ rn. If Pj(d) = Pj(c) then Oji(c) 7^ 
ctji(d). Let Z be the line through pi and Pj(c). If Pj(<2) € i, then 0(27(0) 7^ ot2j(d). 
Finally, iipj(d) £ I, then ay(c) 7^ aij(d). ■ 

Remark 5.14 

Suppose we are given an image point a = V>da„([c]) for some [c] G conf n \clconf n . 
We construct the standard representative from the angles in a as follows. Put 
Pi = (0,0) and P2 — (cos a±2, sin 012). A point is on the line I12 through p\ and 
P2 if and only if 

a.\2 = olu = a2i mod it. (5-1) 
As [c] G conf„\clconf„, we know that there exists pi that is not on the line Z12 . We 



find it by checking Equation 5.1. We construct pi as the intersection of the lines 
In and ?2i- Here In is the line that passes through p\ and has direction an and 
I21 is the line that passes through P2 with direction 021- Any other point pj of the 
standard representative is constructed in a similar way: there is always a pair of 
vertices from the non-degenerated triangle p\P2Pi such that the line through these 
two vertices does not contain pj . 



Let c G CLCONF n . Then all points pi, . . . ,p n from c are on a common line l c 
that makes some undirected angle G (—§,§] with the positive x-axis. Order 
the labels of the points according to the order of the points on l c as encountered 
from left to right, or from bottom to top, in case = ^. This defines an ordered 
n-tuple a(c). The n-tuple cr([c]) is defined as <r(c), for some representative c G [c]. 
Define an equivalence class ^ a ,a on clconf„ as follows: [c] ^ Qj(T [0?] if and only if 
both ajj = opj", and cr([c]) = <r([d]). 

Lemma 5.15 

Consider ij) DAn ■ CLCONF n -> DA„. 

ipDA n is constant on classes o/~ Qj(T . 
fiij TTie map ipDA n '■ clconf„/ ^ Qi(T — ► -DA„ is injective 

Proof. We prove the two claims. 

(i) Let [c],[d] G clconf n such that [c] ~ QiCr [d]. Let Pi(c),pi(d) and pj(c),pj(d) 
denote the z-th and j-th points of some representatives c and d of the classes [c] 
and [d]. Then ctij(c) — aij(d), as all four points are collinear. Moreover, Pj(c) 
is on the right of Pi(c) whenever Pj(d) on the right of Pi(d), as a(c) = a(d). So, 
aij(c) = Oij(d). 

(ii) Suppose that [c] 7^,0- [d]- Let pi(c),pt(d),pj(c) and J3j(d) be as above. If 
57 7^ cEd, then clearly ipT>A n ([c\) 7^ V'DA„([rf])- Assume ct^ = cxd. Then cr(c) 7^ <r(d). 
Let z be minimal such that k :— a(c)i 7^ cr(d)i —: I Then aki(c) = ctki(d) + 7r. 
Therefore Vda„ ( [c]) + ^da„ ( [d] ) . ■ 



We use Lemma 5.13 and Lemma 5.15 to answer the following question. Is it possible, 
given some a = i/^da^c), with c G CONF n , to determine the Voronoi diagram V(c)7 



5.2. From angles back to point configurations. 
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Let the reduced Voronoi diagram, notation V(c), of a point c £ CONF n be the 
Voronoi diagram V(c) up to scaling and translation. A reduced configuration of 
points [c] or a reduced Voronoi diagram V(c) is reconstructible from a list of 
angles a iff [c] resp. V(c) is uniquely determined by a. 

Property 5.16 

V(c) is independent of the choice of c € [c] . 

Notation: By V([c]), for [c] S conf„, the reduced Voronoi diagram of some repre- 
sentative c G [c] is indicated. Note that if one can reconstruct the class [c] € conf„ 
from ^DA„(c), for c e CONF n , then one can also construct the reduced Voronoi 
diagram V([c]). 

Corollary 5.17 

Lei a = V'.DAnQc])) lOT'tt [c] (E conf„. 

(i) If [c] £ conf„ \ clconf„, then [c] is reconstructible. 
(ii) If [c] £ clconf„, then [c] is not reconstructible, except for n = 2. 



Proof. The first claim follows from Lemma 5.13 and Remark 5.14. The preimage 
-0nA „ ([c]) consists of distinct classes of conf„, except for n = 2, consult Lemma 



5.15|. This proves the second claim. 



Corollary 5.18 

V([c]) is reconstructible if and only if either [c] G conf„ \ clconf„, or [c] € clconf„ 
and n — 2 or n = 3. 



Proof. From Corollary 5.17 it follows that [c] is reconstructible in the cases men- 
tioned except for [c] G clconf„ in case n = 3. But in this case, V([c]) consists just of 
two parallel bisectors perpendicular to o^. From 03 = -0da„([c]) we can determine 
the order of the three points on the line, which gives us the labels of the bisector. 
Changing the relative positions of the three points in [c] has no influence on V(c) 
as long as the order of the points is maintained. Suppose, on the other hand that 
[c] € clconf„ for n > 3. It is easy to show that in this case mutually distinct classes 
of configurations in the fiber V'da (M) correspond with mutually distinct reduced 
Voronoi diagrams. H 



Figure 5.2: The Voronoi diagram of four collinear points. 



Example 5.19 

Consider the element d € DA4 defined by olyi = ctn = an — 0123 = ct24 = 0. 
Suppose that we try to construct a configuration S of labeled points having those 
angles. It is clear that any such configuration consists of four points on a common 
horizontal line. An example is given in Figure 5.2. Moreover, it is possible to 



reconstruct the order of the points on the line: from left to right we encounter the 
points pi, P2, P3 and p^ in that particular order. This shows however that from the 
information present in d we cannot determine the ratio 1/IP1P3 1 of the lengths 
of the line segments P1P2 and pip^: whatever this ratio is, the angles cti2, aia and 
Cf23 do not change. The Voronoi diagram of p\ to p± consists of the three vertical 
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bisectors B(pi,p 2 ), B(pi,p 3 ) and B{p 2l p 3 ). Again it is impossible to determine 
the ratio \B(pi,p 2 ) — B(p 2 ,p 3 )\/\B(p 2 ,p 3 ) — B{p 3l p±)\ of the distances between the 
bisectors. This shows that it is impossible to associate a unique reduced Voronoi 
diagram to d £ DA 4 . 



5.2.2 Distinct points and angles in M/7rZ. 

Suppose we start with some c £ CONF n . Assume we have determined a = ip\jA n ( c ) > 
that is, the set of all angles oTJ mod ir between pairs of points (pi,Pj) £ c. We 
introduce another equivalence class on conf n . Let [c], [d] in conf n . The classes [c] 
and [d] are reflection equivalent, notation [c] ~ R [rf] if and only if the standard 
representative c of [c] equals the standard representative d of [d] up to a reflection 
mp 1 (c)=p 1 (d) = (0,0). 

Lemma 5.20 

The map ipuA n ■ (conf„ \ clconf„)/ ~r — > UA n is injective 



Proof. Proceed as in the proof of Lemma |5.13 , but take those standard represen- 
tatives that have a(c), a(d) £ (— f , § ]■ It is possible to put mod tt bars on all o^'s 
that occur in the proof. ■ 



If we take angles between points in R/7rZ, all collinear configurations that have their 
points on a line with angle a £ (—?,§] are mapped to the same a = (a, . . . , a) in 
UA n . So it is impossible to reconstruct the original order of the points on the line. 

Corollary 5.21 

Let a = i/JuAndc]), for [c] £ conf„. 

(i) If [c] £ conf„ \clconf„, then [c] is reconstructible up to a point reflection, 
(ii) If [ [c] £ clconf„ and n ~ 2, then [c] is reconstructible up to a point reflection. 
(Hi) If [ [c] £ clconf„ and n — 3, then [c] is not reconstructible but the reduced 
Voronoi diagram V([c\) is reconstructible 

(iv) If[ [c] £ clconf n and n > 3, then V([c\) is not reconstructible. 




Figure 5.3: V(c) reconstructed from angles mod tt. 



Example 5.22 

Let c = {pi,p 2l pz\, with p x = (2,0), p 2 = (0,0) and p 3 = (1,1). The Voronoi 
diagram V(c) is shown on the left in Figure 5.2. Then a 12 = tt, a 13 = ^ and 
a 23 — j- So Txyi = 0, while aTjJ = «i3 and 023 = 0^23- We get a 'reconstructed' 
configuration c' consisting of p[ = (0,0), p' 2 = (1,0) and p' 3 = (— 5,— |)- This 
configuration c' together with its Voronoi diagram V(c') is shown on the right in 
Figure 5.2. 



5.3. Angle models for small n. 
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5.3 Angle models for small n. 

In this section, we analyze the compactifications CDA n and CUA n for n = 2. 
Moreover, we give a complete description of the image ipriA 3 (CONF3). For points 
on the boundary of ipT>A 3 (CO JVF3), we give a geometric interpretation in terms of 
coinciding points. 

Remark 5.23 

D. G. Kendal has introduced the 'The theory of shape' in a statistical context. As a 
specific example, the shapes of triangles arc analyzed. To compare with this result, 



consult [Kcl| and |Ke2] 



5.3.1 n=2. 

In the directed case, we have the following diagram: 

CONF2 C l 2 xM 2 

n n 

graph(?/W 2 ) C CDA 2 C I 2 x I 2 x R/2ttZ. 

il>DA 2 maps a configuration c = (pi,P2) G CONF n to the directed angle 01.12. For 
simplicity we assume that p\ = (0,0). As c € CONF2, this implies that P2 € 
R 2 \ {(0, 0)}. This space R 2 \ {(0, 0)} is homeomorphic to a doubly open cylinder: 
imagine the omitted point (0, 0) as one side of the cylinder and infinity as the other 
side. 

A classic construction is the following: The blow-up of R 2 at (0, 0) is by definition 
the closed subset X C R 2 x P 1 of all points (x\, x 2 , V\, ya) G R 2 x P 1 such that 
x iV2 = X2lJi- The projection it : X — » R 2 onto the first factor has the following 



properties, compare [Ha 



1. 7r 1 (p2) consists of a single point, if P2 € R 2 \ {(0, 0)}. 

2. ? r- 1 ((0,0))=P 1 . 

3. The points of 7r _1 ((0, 0)) are in 1-1 correspondence to the set of (undirected) 
lines through (0,0). 

We consider CDA 2 , the closure of graph(^DA 2 )- 
Lemma 5.24 

CDA2 is homeomorphic to a half-open cylinder times a plane. 

Proof, (sketch). The plane comes from varying pi, so assume p% — (0,0). Then 
P2 € R 2 \ {(0,0)}, provided that (pi,p 2 ) G CONF 2 . Write p 2 = r(cosa,sina), with 
r £ R>o, and a € R/27rZ. Any configuration ((0,0), a), is added exactly once in 
order to obtain the closure of graph(V>DA 2 )j as lim r ^o ^(cos a, sin a). The punctured 
plane R 2 \ {(0,0)}, is homeomorphic to a doubly open cylinder. Adding all points 
of the form ((0, 0), 0:12), with ayi € R/27rZ means that we attach one full circle to 
that end of the open cylinder that corresponds to P2 = (0, 0). ■ 

Remark 5.25 



Lemma 5.24 demonstrates that CDA n can have a 'boundary'. 



In the undirected case, we consider CUA 2 , the closure of graph(V>uA 2 ) ■ 



Lemma 5.26 

CUA2 is homeomorphic to a Mobius strip times a plane. 
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Figure 5.4: Antipodal identification. 



Proof, (sketch). Start as in the proof of Lemma 5.24 . Write P2 in polar coor- 
dinates, that is, P2 — r(cos a, sin a), with r € M>o , and a 6 M/27rZ. Any undi- 
rected angle «i2 is added twice, once for p2 = lim r _>o r(cos a, sin a) and once for 
P2 = Um r ^o ?*(cos a + tt, sin a + tt). This gives the antipodal identification as indi- 
cated on the left in Figure |5.4| by the little solid arrows. It is allowed to cut the 
space as long as it is eventually pasted back together in the same way as it was 
cut. We make two cuts, c\ and C2, through the side of the cylinder. This enables 
us to perform the antipodal paste, resulting in the rectangle on the right of Figure 



5.4. Pasting back the cuts we have made, results in the Mobius strip, compare e.g. 



[Mn] 



Remark 5.27 

Note that the construction in the proof of Lemma |5.26 is exactly the blow-up of 1 
at (0, 0) that we have discussed above. 



5.3.2 n=3. 

We consider the possible angles mod 2tt between three distinct points p\ , p2 and p 3 
in the plane. That is, we determine the image iprjAsiCONF^) in DA3. Let 1^ denote 
the directed line that passes first through pi and last through pj . The triangle with 
vertices pi, P2 and p 3 is oriented clockwise if p$ is on the right of l±2- The triangle 
is oriented anti-clockwise if p 3 is on the left of Z12. 

Lemma 5.28 

Let c — (pi,P2,P3) G CONF3. Let A denote the triangle with vertices p%, P2 and 
p 3 . For (ai2, ai3, CK23) = ipDA 3 (c) the following holds. 

(i) Lf A oriented clockwise, then < ai3 — «23 < «13 — ot%\ < tt mod 27r. 

(ii) Lf A oriented anti- clockwise, then < a%\ ~ 023 < «2i — cxi3 < tt mod 2-7T. 

(Hi) Lfpi, P2 and ps collinear then 

ai2 = «13 = "23 V a 12 = ai3 = a 32 V a 2 \ = a 13 = a 2 3- 

Moreover, for any image point («i2, "13, "23) € r 4 ) DA 3 (CONF 3 ) either (i), (ii), or 
(Hi) holds. 



Proof. Note that the differences of angles as they occur in the statements in the 
lemma does not change if A is rotated. 

(i) Assume A is oriented clockwise. Rotate A such that ai3 = tt, compare 
the triangle on the left in Figure 5.5. As p\ is on the right of I23 it follows that 
"21 < "23 < tt- 



5.3. Angle models for small n. 
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2 3 
Figure 5.5: The clockwise and anti-clockwise case. 



(ii) Assume A is oriented anti-clockwise. Rotate A such that an — tt, compare 
the triangle on the right in Figure [Tf^. As p\ is on the right of I32 it follows that 
031 < 0/32 < tt. It follows that ai2 — 032 < ctu ~ O31 < t mod 27r, or equivalcntly, 
that 021 — 0:23 < «2i ~ a i3 < 7f niod 27T. 

(iii) Proof by inspection. 

As any triangle on the vertices pi,P2 and p$ is either oriented anti-clockwise, or 
oriented clockwise or degenerate, the lemma follows. ■ 

Boundary points of tpr>A 3 (CONF3) are characterized as follows. 
Corollary 5.29 

Let 03 = (0:12,0:13,0:23) be a boundary point of ipDA-siCONF^). 
( i) If a 3 = ip D A 3 (c) for c G CONF 3 , then a^i = oJi = a^- 
(ii) If 0,3 ipDA 3 {CONF3), then 012 = 013 V 013 = 023 V 023 = 012 + tt. 
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Figure 5.6: Possible values for 0:13, 013, «23 and 0:23, in case 012 = ? and 012 = ^ 
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j, on the left, and, 



Example 5.30 

Figure 5.6 shows the possible values of 0:13 and 023, for 012 
for 012 = ^r, on the right. Note that 

(i) Configurations of points that are not collinear are mapped by i/>da 3 to the 
interior of the triangles and rectangles in the figure. 

(ii) Configurations that are collinear are mapped to the vertices of the big central 
triangle on the left and on the right. 

(iii) Other points on the boundary of the triangles and rectangles do not corres- 
pond to configurations in conf^. 



Example 5.31 

The picture in the middle of Figure 5.6 shows the possible values of 013 and 
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ct23 for a±2 = \- The only configurations in UA3 that are not image points of 
V'UAn (CONFz) are the points where two angles mod tt coincide but not all three. 



Recall that conf3 \ clconf3 consists of the configurations of three distinct, non- 
collincar points up to scalings and transformations. 



Corollary 5.32 

Let a 3 = (x, y, z) 
if 



Then 0,3 — ipDA 3 ([c]) for some [c] G conf3 \ clconf3 if and only 



V 
V 
V 
V 
V 
V 
V 



< x < tt; 

< x < tt; 

< x < tt; 

< x < tt; 

TT < X <2tT 

TT < X < 2lT 

TT < X < 2TT 

TT < X < 2tT 



y < x; 

y > x; 
V>0; 
z > tt - 
y >0; 
y > x; 
y > z; 
y > x; 



x; 



y > z; 
z < x + tt: 
y < x; 
y > z; 
y < z; 
z>0; 
z > x — tt; 
y < z; 



z > 0, 

y < z, 

z > x + tt; 

y < 2tt, 

Z < X — TT, 

y < 2tt; 
y <x, 
z < 2tt. 



z < 2tt, 



y < x 



The boundary of ipDA 3 ([c\) with [c] G conf3 \ clconf3 is depicted in Figure 5.1. 




Remark 5. 33 

In Corollary 5.29| we have classified those points 03 G DA3 that are on the boundary 
of ipuA 3 (CONF 3 ) but not in tpnA 3 (CONF 3 ) itself. These triples of angles can be 
interpreted as the triples of angles that correspond to configurations of three points, 
such that exactly two points coincide. Consider the angles between the points in 
the three configurations presented in Figure 5.S. In the leftmost configuration, 
ai2 = ai3, in the configuration in the middle, = Q?23, while in the rightmost 
configuration, = 032- This interpretation holds, as long as the three points are 



not collinear, compare the first statement of Corollary 5.29 



Example 5.34 

Fix pi = (0,0) and P2 = (1,1)- This assures that 0:12 = \. Suppose that the 
position of j<3 is given by p^it) — P2 + 1. 1 (cost, sin i), for t G [0, 2tt). So «23 = t. 



5.3. Angle models for small n. 
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2 = 3 



1 = 3 



1 = 2 



Figure 5.8: Configurations for two coinciding points. 



Then (ai 3 (i),f) is a curve, that is defined for any t. It is depicted in Figure 5.9, on 
the left. 

Example 5.35 

Fix p\ — (0,0) and pi = (1,1) again. Let P3(t) — v2(cos t, sin t) . I n th is case, 
the curve (oti3(t) , t) is built up out of straight line segments, see Figure 5.9, on the 



right. At t = 2?-, the point P3 passes through p%. At this t, the curve is not defined. 
This t coincides with a swap ot\3 = ot\3 + w. 



5tt 



571 



ai3 



"13 



5tt 



4 4 4 

Figure 5.9: Configurations for moving third point. 



There are some obvious geometric transformations that relate points in the image 
tpr>A 3 (CONF 3 ). We use these relations for analyzing the structure of ipBA 3 ( CONF 3 ) . 
Here we consider DA 3 as a quotient of R 3 . That is, identify points a = (a x , a y ,a z ) 
and b = (b x , b y ,b z ) in R 3 iff a x = b x mod 2ir etcetera. 

Lemma 5.36 

Let (x,y, z) = ipoA 3 (c), with c = (pi,P2,P3) G CONF 3 . Let a be an element of the 
symmetric group S3. The consequence of a permutation a of the labels of the points 
in c is as follows. 



aeS 3 o-(c) e ^da 3 (CONF 3 ) 

y, z) 

(12) {x + tt, z, y) 

(13) (z + ir, y + 7r, x + 7r) 
(23) (y, x, z + ti) 

(123) (z, x + tt, y + n) 

(132) (y + 7T, Z + IT, x) 



geometric action in DA3 

^(7r,0,O) S ai3 — a23 
^(0,0,7r) Sa 12 =a 13 



i?2, oT, 



(t,w,0) 



A T ' indicates a translation over a given vector, an l S ' a reflection in a given 
plane, while 'R' indicates a rotation with respect to the axis through (0,0,0) and 
(1,1,1) over a given angle. 



Proof. A sketch of the six permutations of the labels of the components of p suffices 
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to determine o~(p). The geometric action that corresponds to applying the three 
order two elements on the labels is clear. The following formula, see |Gc| ], gives the 
image r' of the rotation of a vector r through an angle <p about an axis n. 

r' = r cos (j) + h(n ■ r)(l — cos (p) + (r x n sin0). 

In our case, n is the normalized vector 1, 1), while r = {x, y, z}. Applying 

this formula for (j) — — =£- gives r' = {z,x, y}, while 4> = leads to 7-' = {y, z,x}. 
This explains the action corresponding to the order three elements of S3 . ■ 

Lemma 5.37 

If (x, y, z) e ij^DA 3 (CONF 3 ), then also (x + p,y + p, z + p) e ip D a 3 {CONF 3 ) for any 
P€[0,2tt). 



Proof. A rotation of p 6 [0, 2ir) of the baseline, that is the axis with respect to 
which angles are measured, corresponds to adding p to all angles (X12, otis and 

«23- ■ 



Lemmata 5.36 and |5.37 suggest to consider the set of all possible angles between 



three distinct points up to the permutations and baseline rotations introduced in 
the lemmata. Define $^ r :— tpuA 3 ( CONF3 ) / ^{p.r}, the fundamental image of 
4>da 3 (CONF3), as the image tpuA 3 ( CONF3 ) up to permutations and rotations. 

Corollary 5.38 

■0-0^3 ( CO JVF3) up to baseline rotations is obtained by projecting ipr)A 3 {CONFy,) on 
the plane cx.12 + a\ 3 + «23 = 0. 



Proof. This follows directly from Lemma 5.37, as any vector (p, p, p) is orthogonal 
to the plane an + a±3 + 0^3 = 0. ■ 



We can actually construct the projection ir r on the plane P : ct\i + 013 + 0:23 = 
as follows, iti = (1, —1, 0) and u-i = (h, i, —1) form an orthogonal basis of P. Let 
ei and e2 be the corresponding unit length vectors. Then n r is given by 

7iy : JR 3 -> M 2 , 

c ^ (< c, e 1 >, < c, e 2 >). 



^ r (ipDA 3 (CONF 3 )) is shown on the left in Figure 5.1C. 



Remark 5.39 

7ty maps (-7T, 0,0) to p = (ttA^ttB), maps (0, tt, 0) to q = (—t:A,ttB) and maps 
(0,0,71-) to (0,2ttB), where A = ±V2 and B = As a consequence, any 0,3 £ 

ipT>A 3 (CONF3) has a r epres entative mod 2ir that is mapped by tt t in the area shown 



on the right in Figure 5.1C 



Next lemma is the 'projected' version of Lemma 5.36 



Lemma 5.40 

Let (x,y,z) = iJjda 3 ((.Pi j P2,Ps)), with c = (pi,p 2 ,P3) € CONF 3 . The orbit of 
ir r (x,y, z)) under the action of the symmetric group S3 on the labels of p\ 7 pi and 



5.3. Angle models for small n. 



09 



-to 



«13 = «23 



^Cti2 = «23 



«12 = «13 



Figure 5.10: w r (ipBA 3 (CONF 3 )) and six times $>p r . 



P3 is given in the following table. 



ieS 3 


(12) 

(13) 

(23) 

(123) 

(132) 



geometric action, coordinates 



Tftp + S ai3 — a23 

^Ql2— CK23 



B(x + y-2z)) 
(A(tt + x- z), B(tt + x - 2y + z)) 
(A(z-y), B(z + y-2x)) 

(A(y-x), B{x + y-2z-2ir)) 

(A(z-x-tt), B(x + z - 2y - w)) 
{A(y~z), B{2tt + y + z - 2x)) 



Here p, q, r and A, B are as in Remark 5. 3!\ and m is given by m = (^^, —ttB). 
Moreover, 5 = 1, if the second coordinate of ir r {x, y, x) < and S = —1, if the 
second coordinate of TT r (x, y, x) > 0. 

Proof. We can easily compute the coordinates of the points in the orbit by pro- 
jecting the orbit in tpY)A 3 {CONF3) that we found in Lemma 5.36. Concerning the 
geometric action: the reflection matrix S a for a reflection of a point in the line 
y = ax is given by 



1 



l + a A 



I - a 2 2a 
2a a 2 - I 



It is easy to check that the action of (12) boils down to a translation over p followed 
by a reflection in the line ai3 = 0:23- Similarly for the action of the other two 
involutions. The rotation matrix R a for a rotation of a point around the origin is 
given by 



Ra 



cos a — sm a 
sin a cos a 



A rotation of a vector v around an arbitrary center m is given by R a {v — m) + to. 
This allows us to check the action of (123) and (132). ■ 

Corollary 5.41 

<frp r is given by the triangle in the picture on the right of Figure 5.1C with vertices 
(0, 0), —m and p. 
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Proof. Every element of the orbit of S3 acting on a % r (x, y, z) for some (x, y, z) = 
■0da 3 (c), with c G CONF3 lives 'at the same place' in its own triangle in the picture 



on the right of Figure 5.10: the three order 2 elements map ir r (x, y, z) on one of the 
three gray triangles on the other site of the dotted line a 12 — c*23- The two order 
three elements rotate w r (x, y, z) clockwise or anti-clockwise into the two adjacent 
gray triangles. ■ 

Corollary 5.42 

Let P be the plane orthogonal to the vector (1, 1, 1). Then the 2 -dimensional wall- 
paper group p6m is acting on P D ip da 3 {CONF^) . 

Proof. There are order 6, order 3 and order 2 rotations, and 6 axes of reflection. 



This characterizes the 17-th wallpaper group, see for example [CM| 



5.4 The triangle variety T n . 

We want to answer the question: 

are CUA n and CDA n smooth manifolds? 

Moreover, we are interested in their algebraic counterparts. In this section we 
describe an algebraic variety that is very similar to CUA n . We know that CONF n 
is contained in CUA n : by definition, CUA n equals the closure of the graph of the 
undirected angle map: 

CONF n c CUA n c (R 2 )™ x (R/ttZ)(=). 

Therefore, CUA n contains a 'CONF n part' that is smooth. The remaining points of 
CUA n lie above the diagonal A C (R 2 )" consisting of configurations with at least 
two coinciding points pi and pj. We make an algebraic description for R/irZ = P 1 
by taking coordinates (<2y : 1) and (1 : bij), where 

dij — tancfij ■; bij - 



tan Uij 



For simplicity, we consider only the case where fey 7^ on each P 1 , so we work on 
the (aij : l)-chart. We have transformed ip\jA n in a rational map Aslope given by 

Aslope {{xo,yo),---,(Xn-l,y n -l)) = {( { _ % )}o<i<j<n-l, 

X j X/% 

where ((xq, yo), ■ ■ ■ , (x n -i 7 y n -i)) € CONF n . Without loss of generality we assume 
throughout this section that xq — yo — 0. That is, we consider configurations up to 
translation. The dimension of CONF n (M. 2 ) up to translations equals 2n — 2. The 
slope aoj, for i G {1, . . . ,n— 1} is denoted short as c^. The triangle is the triangle 
with vertices (xo, 2/0)1 (xi,j/i), and (xj,yj). The following lemma shows that there 
exists a relation between the x-coordinates of the vertices of and the slopes of 
the lines bounding Ty. Let 

t ij — Cbf^X 1 O, j X j CL^^ X^ I *^ j ' 

Lemma 5.43 

tij — on the (a.y : l)-chart of CONF n . 

Proof. It holds that j/j = aiXi, j/j = cijXj, and j/, — yj = dij(xi — Xj). Substitute 
the former two equations in the last equation. ■ 
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Corollary 5.44 

On the (a,ij : l)-chart we have that CUA n C {iy = 0} for l<i<j<n — 1. 

A question is if equality holds. That is, if the closed algebraic set {tij = 0}i<i<j<„_i 
is contained in CUA n . The answer is no. We prove this later on by means of the 
six-slopes formula: 




Figure 5.11: For four distinct points the six-slopes formula holds. 



Lemma 5.45 (six-slopes formula) 

Let po, pi, P2 and p 3 be distinct points in the plane. Then A = A0123 = 0, where 
A is given by 

A = (<zi - ai 2 )(a 2 - a 23 )(a 3 - 013) - {a\ - a 13 )(a 2 - ai 2 )(a 3 - a 23 ). (5.2) 

Proof. Assume that p = (0, 0) and pi — (1, ai). We compute coordinates for the 
points P2 and p 3 . Let Uj be the line through the points Pi and pj. The lines Iq2 & n d 
l 03 are given by 

I02 ■ 2/2 = a 2 x, l 3 ■ 2/3 = a 3 x, 

and l\2 and I13 by 

li2-y-ai = ai 2 (x-l), h 3 :y~ai = a i3 (x - 1). 

We compute the intersections (x 2 , y 2 ) = I02 H Z12, and (x 3 , j/ 3 ) = /o 3 n Zi 3 : 

= 01-012 x = 01-013 

" 02—012 J 03— ai3 

2/2 = 02X2 2/3 = 0-3X3 

After some formula manipulation the expression follows as a 23 = ■ B 

Remark 5.46 



Some remarks on Equation 5.2 



(i) Interchanging indices 1 «-> 2, etcetera, changes the appearance of the expres- 
sion for Aoi 23 , but does not change the expression itself. 

(ii) By Aijki we denote Aoi 23 with 0, 1, 2 and 3 replaced by i, j, k and I. 

(iii) If no three points of po, pi, P2 and p 3 are collinear, then we can express a 23 
as follows in terms of di, a 2 , a 3 , ai 2 and ai 3 : 

a 3 (ai - oi 3 )(a 2 - a i2 ) - a 2 (ai - ai 2 )(a 3 - a i3 ) 

° 23 ~~ — 7 \7 ^ 7 \7 ^ — ■ 

(ai - ai 3 )(a 2 - ai 2 ) - (ai - ai 2 )(a 3 - ai 3 ) 

Corollary 5.47 

Aijkl = on CUA n . 



Proof. This follows from Lemma 5.45 and the definition of CUA„ 
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Instead of just looking at the zeros of fy = 0, we add the condition that all A^-fe; 
equal zero as well. This leads to the following definition. 

Definition 5.48 

The triangle variety T n is the set of common zeroes of the polynomials tij for 
1 < i < j < (n — I ) and for 0<i<j<k<l< (n — I). Any variable, of the 
form aij or Xi, takes value in R. 

Example 5.49 

The zero set of a collection of polynomials fx. ... , /„ is indicated by V(fi, . . . , /„). 
SoT 4 = F(ti 2 ,ti3,i23,A 123 ). 

Note that the ideal of the triangle variety T n contains: one polynomial iy for every 
triangle with vertices po, Pi and pj; one polynomial A«jy for every quadrilateral 
with vertices pi, pj, p k and pi. We could consider explicitly triangles with vertices 
Pt, Pj and p k by including: 

tijk — ^ik{Xi ^k) &ij(xi %j) ^jk{Xj 37fe)' 

The following lemma shows that this is not necessary however. 
Lemma 5.50 

tijk tij tik ~l~ tjk- 

Proof. This follows directly from the following equations: 

^ik Qi^C'i Qk-^k Q'ik'Ei ^ik^k^ 

t jk = ajxj - a k x k - a jk Xj + a ]k x k . ■ 

5.4.1 Singularities of T n for small n. 

In this section we determine singularities of the triangle variety T n for n = 3 and 
n = 4. For this purpose recall the definition of singularity. 

Definition 5.51 

See / jffcj /- Let Y C A n be an affine variety, and let fi, .... f t £ A — k[x\, . . . , x n ] be 
a set of generators for the ideal ofY. Then Y is nonsingular at a point P S Y 
if the rank of the Jacobian matrix ((dfi/dxj)(P)) is n — r, where r is the dimension 
ofY. IfY is nonsingular at every point, then Y is nonsingular. 

5.4.2 n=3. 




Figure 5.12: Three coinciding points with coinciding directions. 
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Lemma 5.52 

An element c — (xq, X\, X2, ai, a 2 , 012) € T 3 is singular iff all points and all slopes 
coincide. That is: 



xa = xi = x 2 , 
The type of this singularity is A^ . 



a\ = a 2 = ai 2 . 



Proof. T 3 = V(ti2), with ti% = a\X\ — 02X2 — 012^1 + 012X2- T 3 is singular if the 
rank of the Jacobian of £3 is smaller than 5 — 4 = 1, that is if all partial derivatives 
of ti2 vanish. 

variable x\ X2 di 02 012, 

partial derivative a\ — a\2 012—02 X\ —X2 x 2 — x±. 

Solving this system gives the solutions x± = X2 — A ai = 02 — &i2- For the type 
of the singularity: note that t\2 — (a\ — avi)x\ + (012 — 02)X2- From this it follows 
that the singular set is given by Oj = Oj = a,j, and that the singularity is of type 
■ U. cf. @. ■ 

Remark 5.53 

Geometrically, the singularities of T3 correspond to degenerated configurations 
where bot h the three points coincide and the directions between the points coincide, 
see Figure 5.12 . 



5.4.3 n=4. 

Lemma 5.54 

Define 



Ax 


= (a 2 


- ai2)(d3 


- 023) 


- (02 


- a23)(d3 


- 013), 


A 2 


= (ai 


- ai2)(a3 


- 013) 


- (ai 


- 013) (a 3 


- a 23 ), 


A 3 


= («2 


- 023) (ai 


- 012) 


- (02 


- Ol2)(oi 


- ais), 


A 12 


= (ai 


- ai3)(a3 


- 023) 


- (a 2 


- 023) (03 


- 013), 


Al 3 


= (02 


- 012) (a 3 


- a 23 ) 


- (ai 


- Ol2)(02 


- a 23 ), 


A23 


= («i 


- ai 3 )(a 2 


- 012) 


- (ai 


- ai 2 )(a 3 


- a 13 ). 



(1) If A x ^ 0, then 

x 2 = Xx 
If A 2 ^ 0, then 

Xx = x 2 
If A 3 ^ 0, then 

x 2 = x 3 

(*) £t = A 



(012 


- a 13 )(a 23 


- 03) 




Ax 




(aia 


- a 23 )(a 13 


- 03) 




A 2 




(a23 


- a 13 )(ax2 


- ax) 


A 3 



X 3 = XX 



x 3 = x 2 



xx = x 3 



{ail 


- a l3 ){a2 3 


- a 2 ) 




Ax 




(ai2 


- a 23 ){a 13 


- ax) 




A 2 




(023 


- 013) (012 


- a 2 ) 


A 3 



13 • 
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Proof. As £23 — on T n , it follows that (023 — 03)2:3 = (023— 02)2:2. As £12— £13 = 0, 
we arrive at 

(ai 2 - 02)2:2 - (013 - 03)2:3 = 0122:1-0132:1. 

Add (023 — 03) times this last equation to /23 and substitute (023 — 02)2:2 for 
(023 — 03)2:3. This gives that ^12:2 = (012 — 013X023 — 03)2:1, resulting in the 
formula for x 2 . The formula for 23 is obtained by adding (023 — 02) (/12 — /13) to 
/23 and substituting (a 2 3 — 03)2:3 for (023 — 02)2:2. Substitute the expression for 22 
in /12, multiply by A\ and reorder: 

xioi^i + 2:1(012 - a 2 )(ai2 - ai 3 )(o23 - a 3 ) = Xiai 2 Ai, 

xiA = 0. 

The last equation holds in two cases. In the first case, x\ — 0. But from the 
expressions for 2:2 and 2:3 it follows that then 2:2 = and 2:3 = as well, while no 
conditions are set on the slopes. In the other case A = 0, while for 2:2 and 2:3 the two 
formulas in the Lemma hold. The equations in case of A2 ^ and A3 7^ follow by 
relabeling. This proves the first claim. The second claim follows by inspection. ■ 

Corollary 5.55 

CUAi ^ y(Wi 3 ,£ 23 ). 



Proof. This follows from the proof of Lemma 5.54: The variety V(£i2, £13, £23) 
contains = x\ = X2 — 2:3 as a component, without the condition that A0123 
I). ■ 

Remark 5.56 

Corollary 5.55 and its proof explain why we have added the Ay^'s in the definition 



of triangle variety T n , cf. Definition 5.48. 



Lemma 5.57 

A configuration c G T4 is singular iff, up to relabeling, both Xq = x\ — X2 and 

Ol = 2 = Oi 2 . 



Proof. Consider the Jacobian J4 of T4 and apply Definition |5.51| . The Jacobian of 
T4 where the nine variables are in the order xi, X2, 2:3, 01, 02, 03, 012, 013, 023 is 
given by: 



A'i 


- Oi2 


Ol2 - 2 









Xl 


-x 2 





x 2 - xi 





\ 


Ol 


- Oi 3 





O13 




a-3 


Xi 





-X3 





2:3 - Xl 










02 - 2 3 


023 




a-3 







-X3 








X3 - x 2 


V 















-A x 


-A 2 


~A 3 


±Ai2 


±^13 


±A 2 3 J 



The dimension of T4 equals 2-4 — 2 = 6. Therefore, a configuration is singular 
iff the rank of its Jacobian is equal to or smaller than 2. Denote by m(i,j,k) the 
submatrix of J4 consisting of columns i, j and k where the fourth row is deleted. We 
distinguish several cases, by considering the number of distinct clusters of coinciding 
points in c. 

(i) [> 3 clusters] We may assume that xi ^ 0, 12 / 0, 13 ^ and xi ^ x 2 - 
Consider m(5, 6, 7): 



i(5,6,7) = 




The determinant of m(5, 6, 7) equals 2:22:3 (2:2— x{) 7^ 0. This implies non-singularity. 
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(ii) [2 clusters: (2,2)] We may assume that x\ = and that 12 = 13^ 0. In 
this case the determinant of m(5,6, 7) equals (x 2 ) 3 , which implies non-singularity . 

(iii) [2 clusters: (3,1)] See also Figure [5.13| . Assume that for c 6 T4 it holds 
that Xq = x\ = x 2 ; that <X\ = a 2 — a\ 2 and that x 3 ^ Xq. As a consequence 
13 = 013 = 023- It follows that the first and last row J4 vanishes. Therefore, c 
is singular. If we only assume that xq — x\ — x 2 and that x 3 7^ xq, then the 
determinants det(m(l, 6, 8)) = (x 3 ) 2 (a\ — an) and det(m(2, 6, 8)) = (iE3) 2 (ai2 — a 2 ) 
show that the condition a\ = 02 = ai2 is necessary for a singularity. 

(iv) [1 cluster] It holds that xq — x\ — x 2 — x 3 = 0. Consider J4. First assume 
that the last row of J4 equals zero. As A = on T4 this means that we are looking 
for solutions of the system of equations 

A 1 = A 2 = A 3 = A 12 = A 13 = A 23 = A = 0. 

Applying Mathematica's Reduce function results in 10 reductions. Each of these 10 
reductions contains equations of the form a± = a 2 = an- We are left with the case 
that some Aij ^ 0. If rank( J4) < 2, then the determinants of all 2 x 2 submatrices 
of m(l,2,3) vanish. This gives a system of nine equations that can be reduced to 
the six equations of the form a\ = a 2 = a 3 = an = ^13, that is, five slopes are equal. 
Again, each of these equations contains a condition of the form a\ = a 2 = an- B 




3 



Figure 5.13: Typical singular configuration: three coinciding, collinear points with 
a fourth point. 



Remark 5.58 

The singularities of T4 are closely related to the singularities of T3. The singularity 
in the (3,1) case is of type ^loo just as for T3. But if all points coincide a more 
complicated singularity occurs as several A^ singularities 'meet': one can move 
any of the four points away in such a way that the three remaining points are as in 
the configuration of Lemma 5.52. 



It is still an open question whether CUA4 — T4. Maybe we need to add some 
relations or inequalities to T4 to obtain equality? 
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Chapter 6 

Continuity of the Voronoi 
map. 



In Chapter 4 we have introduced the compactification CDA n of the graph of the 
angle map, applied to n distinct point in the plane. By means of the extended 
definition of Voronoi diagram, as introduced in Chapter 3, we are able to consider 
the Voronoi diagra m VO jn) of a data set 7„ e CDA n . The main result of this 
chapter is Theorem |6.1l[ It states, up to a compactness condition, that two data 
sets 7„ and r\ n in CDA n that are Euclidean close, have Voronoi diagrams that are 
Hausdorff close. That is, the Voronoi diagram generated by a set of not necessarily 
distinct points and the angles between those points does not change dramatically if 
we perturb both the points and the angles only slightly. 



6.1 Introduction. 



In Chapters 3 and 4 we have extended the notion of Voronoi diagram to configu- 
rations of points together with the angles between the points. This allows us to 
consider Voronoi diagrams of point sets that include coinciding points. More pre- 
cisely, in Section 5.1.3, we have introduced a compactification CDA n , as the closure 
of the graph of the angle map 



V>da„ : CONF n -> DA 



The angle map V'da,, maps a set c of n distinct points in the plane to the angles in 
R/27rZ between the points in c. 

In this chapter we call the elements 7„ G CDA n , consisting of n points in K 2 and 
(™) angles in R/2ttZ, data sets. With any data set j n € CDA n , we associate a 
Voronoi diagram V(j n ) in Section 3.2. 

One can wonder whether the data sets are a robust way of storing Voronoi diagrams. 
That is, we consider the following question: if we perturb a data set 7„ 6 CDA n 
slightly, how does the corresponding Voronoi diagram V(7n) change? By perturbing 
a data set 7„ € CDA n we mean that both the points components and the angle 
components of 7„ are allowed to be perturbed slightly, as long as the perturbed data 
set 7„ is again in CDA n . A first result is given in Theorem 6.1, that is a kind of one 
point continuity theorem. It states that a point x that was on the one-skeleton of a 
Voronoi diagram V("f n ) before perturbing, cannot be too far from the one skeleton 
of the Voronoi diagram of the perturbed data set. 
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Before we are able to give a more general continuity result, we first need a suitable 
metric on the set of Voronoi diagrams. A metric that is often used to compare 
pictures, is the so-called Hausdorff metric. Two subsets A and B of R 2 are close 
to each other in the Hausdorff metric, iff the maximal Euclidean distance from any 
point b 6 B to the set A is small and vice versa. For a precise definition, see Section 



f5.4 



In order to prove a continuity theorem with respect to the Hausdorff metric, we have 
to add some restrictions to our underlying point configurations. We add four so- 
called camera points, points that are very far away around our configuration. Now 
we restrict the configurations we are interested in to a relatively small bounded 
subset of the plane. If we only change the configurations within this subset, the 
camera points guarantee us that the Voronoi diagrams of the configurations are not 
changing outside a second, larger bounded subset of the plane. But then we are 



ready for proving the main result of this chapter: Theorem 6.11 states that, up to 
the compactness condition, two data sets that are Euclidean close, have Voronoi 
diagrams that are Hausdorff close. Note that it is essential that the continuity 
theorem holds on the whole of the chosen compact subset of CDA n : that is, the 
points added to the graph of ipr>A„ by taking the closure are essential as there 



availability is used in the proof of Theorem 6.11 



6.2 Preliminaries. 

Throughout this chapter we compare data sets 7„,?7 Tl G CDA n . Write 

In = (pi(7„), . . . ,p„(7„),ai 2 (7„), . . . ,a ( „_ 1)n (7„)), 

and similarly r\ n . We say that n n is within distance 8 of 7„, notation d(^f n , rj n ) < 5, 
iff 

max ||pfe(7„) -Pk(ih.)\\ < & A max \\ a ij(ln) - <Xij(Vn) II < 5. 

l<k<n l<z<j<n 

Here, ||.|| denotes the ordinary Euclidean distance on M 2 and R/2irZ. 



From Section iA we recall the definition of Voronoi diagram in terms of not neces- 
sarily non-coinciding points in the plane and directed angles between those points. 
This enables us to define for any 7„ 6 CDA n a Voronoi diagram V(~f n ). In Sec- 
tion we have seen that the Voronoi cell of a generator pi can be written as the 
intersection of all Voronoi half-planes vh(j>i,pj). We use this characterization to 
introduce the Voronoi cell of a point Pi(7«). 

Fix two points pi — Pi{~f n ) and pj — Pj(jn)- The bisection point b(pi,pj) is 
the point \(pi + Pj)- If Pi 7^ Pj, the bisection point is just the middle of the line 
segment PiPj. If Pi = Pj then b{pi,pj) coincides with the double point Pi = pj. 
The perpendicular bisector B(pi,pj) is the line through b(pi,pj) perpendicular 
to the angle ctij = aij{^f n ). Let n be any non-zero vector, pointing in the direction 
otij. The Voronoi half-plane vh(pi,Pj) is the half-plane defined by 

n • (x - b x , y - b y ) < 0. 

The Voronoi cell V(pA is defined as 

V{pi) = f]vh(pi,pj). 

A point x is on the Voronoi edge e(pi,pj) iff it is on the intersection of the Voronoi 
cells V(jpi) and V(pj), that is 

x 6 e(pi,Pj) x G V(pi) n V(pj). 



6.3. One point continuity. 
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The Voronoi diagram is the family of subsets of R 2 consisting of the Voronoi cells 
V(pi) and all of their intersections. The one-skeleton of the Voronoi diagram is the 
union of the boundaries of the Voronoi cells. Note that the notion of one-skeleton 



equals the notion of shape of a Voronoi diagram, introduced in Section 4.4.1 



6.3 One point continuity. 

In this section we prove the following theorem. 
Theorem 6.1 

Let 7„ £ CDA n be a set of data representing n generators in R 2 and Q) angles 
in R/27rZ between these generators. Suppose that x is on the one-skeleton of the 
corresponding Voronoi diagram V(j n ). Then for every e > there exists S > 0, 
such that when we perturb 7„ by not more than 8, the one-skeleton of any perturbed 
diagram is within distance e of x. 

In order to prove this theorem, we perform the following steps: 

1. We decrease e a finite number of times: for every pair of generators (p,q) 
whose bisector B(p, q) does not pass through x we decrease e, if necessary, in 
such a way that 

d{x,B(p,q)) > 2c. 

2. For every pair of generators (p, q) we show how to choose S pq . If x £ B(p, q) 
before perturbing then after perturbing by S pg , 

B(p,q)nB e (x) ? 0. 

If x ^ B(p, q) before perturbing then after perturbing by S pq , 

B(p,q)nB e (x) = 0. 

Of course this last statement can only hold if we adjust e as indicated above. 

3. Set 5 := min Pjg S pq . 

These steps gives us enough control on the bisectors to complete the proof. We first 
show that the bisection point of two generators p and q can not get more perturbed 
than the generators themselves. 

Lemma 6.2 

Let p and q be two points in the plane. Letb(p,q) = (p + q)/2. 7/max( \\p— p\\, \\q — 
q\\ ) < S, then 

\\b(p,q)-b(p,q)\\ < 5. 
Proof. ||6(p,«)-6( P) g)|| - |||(p + q) - \(j> + q)\\ - |||(p - p) + ±(q - q)\\ < 

\\\P-p\\ + ill9-g|| < s - ■ 

We apply Lemma |6^2| in choosing 8 pq such that a bisector B(p, q) that passes through 
x before perturbation still passes through B t (x) after perturbation. 

Lemma 6.3 

Let p, q and a pq be data from j n £ CDA n such that x £ B{p 1 q). Write r = 
\\x-b{p,q)\\. Let 

S pq = min(— -— ,1). 
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//max( \\q — q\\, \\a pq — ap q \\) < S pq then the bisector B{p,q) passes through 

B e (x). 

Proof. We treat the perturbation of the bisector angle a pq and the bisection point 
b(p, q) separately and add up the maximal effect. Let l &pq be the image of a rotation 
over ±S pq around b(p,q) of the line B(p,q). Then x is at distance rsm.S pq of 
la pq - Moreover, from Lemma 6.2 it follows that \\b(p, q) — b(p,q)\\ < S pq . Let 
d = \\x — B(p,q)\\. Adding up the two effects gives that d < 5 pq + rshiSp q . We want 
d to be smaller than e. As S pq is a majorant for sm.S pq , we are safe if we ensure, as 
claimed, that 5 pq < e/(l + r). ■ 

The next lemma gives a value of S that ensures that a bisector that misses x before 
perturbation, stays away far enough from x after perturbation. 

Lemma 6.4 

Let p, q and a pq be data from j n € CDA n . Let x £ R 2 be such that 

\\x-B(p,q)\\ > 2e. 

Write r = \\x— b{p, q)\\ . The angle between the bisector B{p, q) and the line through 
b{p, q) and x is denoted by 7 € (0, Let 

7r — 2e 

m = TT2 - ' 

7/max( \\q~ q\\, \\a pq - a pq \\ ) < 6 pq then \\x - B(p,q)\\ > e. 



Hp, 1) 




B(p,q) 



Figure 6.1: d is the distance from the bisector B to x. 



Proof. First we show that indeed 5 pq > 0. As 2e < rsin7 < rj it follows that 
7r — 2e > 0. We treat the perturbation of the bisector angle a pq and the bisection 
point b(p, q) separately and add up the maximal effect. Let l & be the image of 

T he d istance d of x to 
see also Figure 3T. Moreover, from 
5 n 



a rotation over ±5 pq around b{p,q) of the line B(p,q) 



la pq s given by rsin/3. Here (3 := 
Lemma ^2 it follows that \\b(p,q) 



7 ~ o pq , 
-b(p,q)\\ < 



j pq This implies that the composite 
minimal distance, that is the distance of x to B(p, q) could become as small as 
d\ := rsin(7 — 5 pq ) — Spq. We want this distance d\ to stay bigger than e. As 
5(7 — S pq ) is a minorant of sin(7 — S pq ) we are safe if we ensure, as claimed, that 



S pq < 



jr — It 



The following corollary is obvious but useful. 



6.4. The distance between two Voronoi diagrams. 
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Corollary 6.5 

Let x,p,q,a pq and 5 pq be as in Lemma Suppose that (p,q,a pq ) is perturbed by 
a vector of length at most 5 pq . For every point y G B e (x) it holds that 

y G vh(j>, q) before perturbation y G vh{p, q) after perturbation, 

where vh{p 1 q) denotes the Voronoi half-plane of generators p and q containing p. 



We are ready to prove Theorem 6.1 



Proof of Theorem 6.1. Let x be on the 1-skeleton of V(j n ) before perturbation. 
Let e > 0. Adjust e a finite number of times: for every pair of generators (p, q) 
whose bisector B(p, q) does not pass through x we set 

e = min(e, -d(x,B(p,q))). 

This ensures that for every such pair (p, q) it holds that d(x, B(p, q)) > 2e. 

Next determine a S pq for every pair of generators in such a way that i) if x G 
B(p,q), then the condition in Lemma O is fulfilled; ii) if x ^ B(p,q), then the 



condition of Lemma 3.4 is fulfilled. Let 5 be the minimum of all these S pq , that is, 



P := {p £ 7„ | fl B t (x) ^ 0, before perturbation}. 

Consequently, P has at least two elements. Say P = {pi, . . . ,p m }, for some m > 2. 
Before perturbation it holds that 



B e (x) = \J(V( Pl )nB t (x)). 



Suppose that, after perturbation, B e (x) C V(s), for some generator s. We show 
that this leads to a contradiction. First note that s E P leads to a contradiction 



immediately: we can apply Lemma 3.3. This implies that after perturbation some 
bisector 6(s,p,) passes through B t (x). We are left with the case that s ^ P. Fix 
a point y S B e (x) such that y is not on any bisector before or after perturbation. 
Then 

y G V(pi) before perturbation, say, 
y G V(s) after perturbation. 

We concentrate on the bisector B(pi, s) now. If x £ B(j>i, s) before perturbation, 
we apply Corollary |6.5| : any point in B e (x) stays in the same half-plane vh(p%, s). As 
a consequence, y G vh(pi, s) before and after perturbation. This is in contradiction 
with our assumption that y G V{s) after perturbation. So x G B{pi,s) before 
perturbation. This implies that 

B(p u s)nB e (x) ^0, 

after perturbation: a contradiction with our assumption as well. ■ 



6.4 The distance between two Voronoi diagrams. 
6.4.1 Hausdorff distance. 

If we want to compare two Voronoi diagrams, a suitable notion of distance is the 
Hausdorff distance: two sets A and B are within Hausdorff distance r iff r is the 
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smallest number such that any point of A is within distance r from some point of B 
and vice versa. Let us give this definition more formally. Suppose we have a metric 
space (X, d). For A C X and r > we define the open neighborhood N r (A) as the 
set 

N r (A) := { y | d(x,y) < r, for some x G A}. 
Think of N r (A) as the territorial border of A. 
Definition 6.6 

Let (X, d) be a metric space. The Hausdorff distance h(A < B) between two 
subsets A,B C X is defined as h(A, B) := inf{r : A C N r (B) and B C N r (A)}. 



The Hausdorff distance h defines a metric on the set of nonempty compact subsets 
of (X,d), see [Ca]. Note that 

h(A,B)<r AcN r (B) A B c N r (A). 



6.4.2 A compactness condition. 

In this section we put a compactness condition on the configuration space of data 
sets. It consists of two ingredients. Suppose that {pi, ■ ■ ■ ,p n } denotes the un- 
derlying point set of some data set 7„ £ CDA n . First we restrict the domain of 
{pi, . . . ,p n } to a closed disk U . M oreover, we add four camera points c%, c%, C3 and 
C4 far away outside hi. Lemma |6.8| shows that we can choose hi and ci, C3 and C4 
in such a way that it is guaranteed that the Voronoi diagram of j n , extended with 
Ci,C2,C3 and C4, changes inside some compact subset of the plane only, provided 
that we perturb {pi, . . . within hi. 

For U we choose the unit disk, that is U := {x E M 2 | ||a;|| < 1}. 
Definition 6.7 

Let 7„ £ CDA n be a set of data representing n generators {p\, . . . ,p n } in Li and 
(2) angles in R/27rZ between those generators. Let furthermore 

Cl = (-JV,0), c 2 = (0,7V), 
c 3 = (N,0), c 4 = (0,-iV). 

for some N ^>> 1 be the so-called camera points. The Voronoi diagram with 

camera points Vjv(7n) is the Voronoi diagram of the data set 7 n , which consists 
of the generators {pi, . . . ,p n } U {c%, . . . , C4} and the (™J 4 ) angles between those 
generators. 

Lemma 6.8 

Let N > 2 + 2 v / 2 and let B% := {x e M 2 | d(x, 0) > N}. Then V N {j n ) n B% does 
not change, provided that j n is perturbed inside hi. 

Proof. Let z = i?(cos 4>, sin (f>) be an arbitrary point in B^. It is enough to show 
that z is always closer to at least one camera point than to 14. Using the symmetry 
we can assume, without loss of generality, that < <f> < ?. So it is enough to 
show that d(z, C3) < d(z, IX) or, equivalently, that d(z, IX) 2 — d(z, C3) 2 > 0. Now, 

d(z,c 3 ) 2 = R 2 sin 2 <t>+ (N - Rcoscj)) 2 , 
= R 2 + N 2 - 2N Rcoscj), 

is maximal when <p = ~ . From now on assume that 



d(z, c 3 ) 2 = R 2 + N 2 - V2NR. 



6.5. Continuity of the Voronoi map. 



83 



As d(z,U) = R — 1, we have to show that 

d(z,U) 2 -d(z,c 3 ) 2 = (R 2 -2R+l)-(R 2 + N 2 -V2NR) > 0. 

Substituting R = N + X , where X > 0, gives that 

d(z,c 3 ) 2 - d(z,U) 2 = N((V2-l)N~2) + (V2N-2)X + l. 

Now, (V2 - l)iV - 2 > 0, whenever AT > 2 + 2\/2, while (V2A - 2) > whenever 
TV > y/2. This proves the lemma. H 

6.5 Continuity of the Voronoi map. 
6.5.1 The Voronoi map. 

In our original configuration space, generators live in M 2 . But under the compact- 
ness condition above, generators all live on the unit disk U and we have four addi- 
tional camera points, c\, C2, c 3 and C4. That is, we have a restricted configuration 
space 

CDA" := {ci, . . . ,c 4 ;pi, . . . ,p„;ai2, . . . ,a(„ + 3)( Il+ 4)} ; Pi € W; ay G M/27rZ. 

Let V(n) C -P(R 2 ) be the set of one-skeletons of Voronoi diagrams, defined by data 
sets 7„ € CDA^. Denote by d be the Euclidean metric on CDA^ and by h the 
Hausdorff metric on P(M 2 ). We show that the map fv, 

f v : {CDA u nl d) -> (V(n),h), 

In l-> V(7„), 

that maps a data set 7„ to its Voronoi diagram, is continuous. By definition, this 
means that we have to show that 

V 7n G CDA%, Ve > 0, 38 > 0, Vr?„ G CDA% 

d{ ln ,n n )<5 => h(V(~/ n ),V(r) n ))<e. (6.1) 



In Theorem 6.1 we have proved that 



V(7„,a; G ^( 7 „)), Ve > 0, 3<5 > : d{ ln , Vn ) < 6 d(x,V(n n )) < e. 

In this formula, 5 really depends on both the particular position of x on V("f n ) and 
the particular diagram V(j n ). The uniform version of this claim is given by 

Ve>0, 35 >0, V(7 n)I e%)): 

dhn,Vn)<S =► dCa;, V(ffc)) < e. (6.2) 

Remark 6.9 



It is clear that (6.2) does not hold if we regard the one-skeleton of Voronoi diagrams 
that correspond with data sets that do not contain the camera points. Suppose for 
example that our data consist of two points p and q that are close to each other. If 
we fix p, move q slightly, thereby changing a pq slightly as well, some point that was 
far away on the bisector B(p, q) before moving q, is on a big distance from B(p, q) 
after disturbance. 
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6.5.2 {(7 n ,x) | x G V(7„))} is closed in CDA n x M 2 . 

In this section we show that the set {(7m x) \ x G V(^ n )} is closed in CDA n x R 2 . 
Or, equivalently, that the complement {( 7n , x) \ x £ V(^ n )} is open. But this means 
exactly that we have to prove the following: 

V 7 „, x<£V(i n ), 3e>0, 36 >0: d(7„,?? n ) < 6 d(x,V[r) n )) > e. 

We do so in the following lemma. 

Lemma 6.10 

Let 7„ be a data set representing n generators in the plane and Q) angles between 
those generators. Suppose that x is in the interior of Voronoi cell V(p), for some 
generator p € j n . Then there exists e > and 6 > such that the following holds: 
if we perturb 7„ by not more then 6, the one-skeleton of any perturbed diagram is 
at distance at least e of x. 

Proof. [Determine e.] Start with some e > 0. For every pair of generators (p, q) 
whose bisector B(p, q) does not pass through x we set 

e = min(e, \d(x, B(p, g))). 

This ensures that d(x, B(p, q)) > 2e for every such pair (p, q). 

[Determine 5(e).] As a consequence, we can apply Lemma [TJ for every pair of 
generators (p, q) whose bisector B(p, q) does not pass through x. The lemma gives 
us S pq such that any perturbation of the data by 6 pq implies that d(x, B(p,q)) > e 
after perturbation. Let 6 := min^ pq ^ 6 pq be the minimum of all these 5 pq . 

[Show that this 5 works.] Suppose that the claim of the lemma is not true. This 
means that after perturbation there are generators q\ and q 2 such that both 

B e (x) n V( gi ) + 0, and B e (x) n V(q 2 ) ± 0, 

hold. It follows that, after perturbation, B(qi, q 2 ) fl B e (x) ^ 0. Because of our 
choice of 5 and e this implies though that x G B (91,52) before perturbation. As 
x G V(p), before perturbation, we conclude that qi ^ p and q 2 7^ p. But then there 
exists y G B t (x) such that 

V e h(p,qx), 

before perturbation, and 

y e h(q 1 ,p), 



after perturbation. This in contradiction with Corollary 3.5 however 



6.5.3 Proof of the continuity. 

In this section we prove the following theorem: 
Theorem 6.11 

Let 7„ G CDA 1 ^. Then the Voronoi diagram V(r) n ) of any data set r\ n G CDA^ that 
is Euclidean-close to 7„, is Bausdorff -close to the Voronoi diagram V("f n ): 

Ve>0, 36 >0, : d(7n,r?„)<<5 => W,%))<e 

Proof. The proof is divided into a number of steps. 

Step 1. It is enough to prove the following assertion. 

Ve > 0, 36 > 0, V 7 „ G U, Vir G F( 7 „), Vi)„eW: 

d(jn,Vn)<6 => d(x,V(rj n ))<e. (6.3) 
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In (| 

fact 



the implication holds for all choices of x on V(j n ). This implies that in 



d(7n,f?n) < 5 



d(x, V(r) n )) < e Vx. 



This together with the fact that we can interchange 7n and r\ n defines the Hausdorff 
distance and therefore implies the claim of the theorem. 

Step 2. Construct a con vergent sequence. 

Suppose that Assertion 6.3 does not hold. That means that the negation must be 
true, where we replace V<5 by Vm: 



3e > : VmeN, 3( 7n , x G V(-y n )), Brj n : 

(d{jn,Vn)<— A d(x, V {r) n ))> e). 
m 

Fix such e and call it eg- We can find for every m G N a triple i T 



(6.4) 



(6.5) 



such that Assertion 6.4 is true. 



G 



d(cc,0) < AT}. Suppose that 



Step 3. x m G B N , Vm. 

Recall that denotes the disk Bm 
x m <£ B N for some x m G V{j!?). Then Lemma |6j tells us that x m € V(rj n ), for all 
7] n G CDA^. So Assertion 3.4 can never be true. 

Step 4. Remarks on compactness. 

7 n and r) n both live in the compact set 

W x {ci,c 2 ,C3,c 4 } x (R/2ttZ)("^ 4 ). 

We have proved in Lemma [6.10 that the set {(7„,a;) | a; G V(^f n )} is closed. This 
means that we get a sequence (t m ) of triples of the form ( |6.5| ) that live on a com- 
pact set. So, (t m ) has some convergent subsequence t m (k) that converges to, say, 
{ln-,x,fj n ). Note that x G V(%). 



Step 5. A contradiction by combining Assertion 6.4 and Theorem 6.1 



Because of convergence, V(j n ) = V(fj n ), so x G V r (f? n ). Apply Theorem 6J with 

• x G V(77 n ), 
l . 



The theorem gives us a So such that 

d(fj n ,0 n ) < S d{V(9 n ),x) < ie . 

Choose m so big that the following condition hold: 

d( {l™iX m ,r]™),{%,x,f) n ) ) < mm(S , |e ). 

We show that this leads to a contradiction on d(x m , V(t]™)). 

(i) The sequence t m is constructed in such a way from Assertion 6.4 that 

d(x m ,V( V ™)) > 6 . 



(ii) We have chosen m so big, (6.7), that 

d(x m ,x) < 



(6.6) 
(6.7) 

(6.8) 
(6.9) 
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(iii) But m is also big enough such that d(fj n ,r)™) < Sq. Therefore we conclude 



from mm, with 9 n := that 



(6.10) 



Now ( p.9| ) combined with ( p. 1C| ) gives a contradiction with (6.8). 
This proves the theorem. 



Chapter 7 

Clickable Voronoi diagrams 
and hook compactifications. 



We define a compactification space XAH n of the configuration space of n distinct 
points in the plane, by considering data elements of pairs of points and triples of 
points. For every pair of points we write down the angle mod tt between the two 
points, and for every ordered triple (pi,Pj,Pk) of points, we specify a hook . 
This hook expresses how to construct the point pk given pi and pj. Now XAH n 
is defined as the closure of the image space of all angles and hooks on n distinct 
points. We show that configurations that are added by taking the closure have a 
natural nested structure, easily revealed by analyzing the hooks. The main result 
of this chapter is an explicit construction establishing XAH n as the graph of a 
function. This construction shows that XAH n is a smooth manifold. If we replace 
the angles mod ir by angles mod 2ir we get a compactification space XEDAH n . This 
space is isomorphic to the manifold with corners FM2 (n) , introduced by Kontsevich- 
Soibelman. On XEDAH n we define clickable Voronoi diagrams. 



7.1 Fulton-MacPherson related models. 



In this chapter we describe a compactification XAH n of the configuration space of 
n distinct points in the plane. The basic idea of the compactification consists of 
considering the geometry of all pairs and all triples of distinct points. In section 
7.2 we give an overview of the contents of this chapter: it describes in an informal 
way the construction and some properties of this compactification XAH n . And it 



explains how clickable Voronoi diagrams are related to XAH n . Sections 7.3 to 7.9 



contain the mathematical constructions referred to in Section \l while Section 
7.1C| is devoted to concluding remarks. We start however in this section by recalling 
the famous compactification due to Fulton an MacPherson, and by presenting a 
related compactification as described by Kontsevich and Soibelman. These two 
compactifications have some important features that we mimic and extend in our 
compactification XAH n . 



7.1.1 The Fulton-MacPherson compactification. 

The following compactification is defined in terms of a nonsingular algebraic variety. 
The authors note however, cf. [FM], page 188/189, that the same constructions work 
for complex manifolds, as well as for real manifolds. So, let X be a nonsingular 
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algebraic variety, and let 

CONF n (X) = {(p X) . . . )Pn ) eX n | Vi ± Pj if i + j}, 



be the configuration space of n distinct labeled points in X, cf. Definition |5-4 For 
a subset S of the set of labels {1, . . . , n}, the small diagonal As is given by 

A s = {(pi,.. .,p n ) £ CONF n {X) | Pll = ■ ■■ =p im for any i 3 e S}. 

For any such S with at least two points, let B1a(A s ) be the blow-up of X s along 
As- For the definition of blow-up, consult [Hafl. There is a natural embedding 



CONF n (X) C X n x Y[ B1 A (A S ). 

|S|>2 

The Fulton-MacPherson compactification X[n] is defined as the closure of 
CONF n (X) in this product. We list some properties that are stated respectively as 



Theorems 1, 2 and 3 in FM|. The definition of divisor can be found in [Ha], 
f. X[n] is nonsingular. 

2. For n > 2, X[n] is the closure of CONF„{X) in the product of the Bl A (X s ) 
for S C {1, . . . , n} of cardinality 2 and 3. 

3. For each S C {1, . . . ,n} with at least two elements, there is a nonsingular 
divisor D(S) C X[n] such that 

(i) The union of these divisors is X[n]\CONF n {X). 

(ii) An intersection of divisors D(Si) D • • • H D(S r ) is nonempty iff the sets 
Sk are nested in the sense that each pair Si and Sj is either disjoint or one 
is contained in the other. 

Im porta nt for us is the 'first geometric description' of A[n], as presented in Section 1 



of FM |. A configuration in X[n] is called degenerate, if not all points in the 
configuration are distinct. Any such configuration c is described by a set of data 
given the locations of the points, and if two or more points coincide at say x, then a 
screen is specified for the set S C {I , . . . , n} of labels of the coinciding points. The 
data elements describing a screen for S at x consist of a labeled set of points x a in 
the tangent space T x of x, such that not all points in x a are equal. This process is 
repeated for x a until all points in the configuration c are separated in some screen. 
Think of repeatedly zooming in at clusters of points until all points are separated, 



cf. Figure 7.9. The divisor D(S) of A[n], mentioned in the property above consists 
of all configurations c in X[n] that have a screen that contains exactly the points 
x a for a in S. 



7.1.2 The Fulton-MacPherson operad. 



In [Ko and [KS| the Fulton-MacPherson operad FMd is discussed. A definition 
of operad can be found in |Kcfl . As we consider only point configurations in M 2 , we 
set d = 2 in the sequel. The operad FM 2 = {FM 2 (n)}> is defined as follows, cf. 



[KS|: 



1. FM 2 (O):=0. 

2. FM 2 (I) := point. 

3. FM 2 (2) := conf 2 = S 1 . 

4. For n > 3, the space FM 2 (n) is a manifold with corners, its interior is conf fl , 
and all boundary strata are certain products of copies of conf„/ for n' < n. 

A manifold with corners is a topological space that is locally homeomorphic with 
R" . The space FM 2 (n), n > 2 is defined explicitly as follows. 



7.2. Informal introduction. 
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For n > 2, the manifold with corners FM 2 (n) is the closure of the image o/conf„ 
in the compact manifold (S )\ 2 ) x [0, +oo under the map 

[(pl,...,P„)] H-> ((atij)l<i<j<n,Pij), 

where i, j, and k are pairwise distinct indices, ctij € R/27rZ, and Pjj — j^E^j ■ 

That is, Kontsevich and Soibelman define the space FM 2 in terms of data elements 
for pairs and triples of points in any reduced configuration c S conf„ : for any pair 
of points they write down the directed angle between the two points, while for every 
ordered triple {pi,Pj,Pk) of points they specify the ratio /?]| of the line segments 
PiPk and pipj. 



7.1.3 Combining the models. 

From the Fulton-MacPherson compactification we borrow the screen model, while 
from Kontsevich-Soibelman we borrow the angles and ratios. We extend, however, 
the ratios to hooks, by marking both the ratio [if* of the line segment piPk and piPj 
and the angle afj between the line segments. In this way every hook has two repre- 
sentatives. One with a positive ratio and one with a negative 'ratio' /3|| . In the 
latter case we just add tt to the hook angle afj . These two distinct representatives 
are identified by an equivalence relation This will be important for getting a 
smooth model. Another adjustment that we make for obtaining smoothness is that 
we write down the angle between two distinct points up to multiples of tt instead 
of up to 2tt. To summarize, we get a map 

ipAHn '■ conf„(R 2 ) -> AH n := (E/ttZ)(=) x (([-co, oo] x M/27rZ)/ ^ k f^ 3 \ 
and define XAH n as the closure of the image of conf„ in AH„ . The data elements 



and spaces introduced so far are discussed in detail in Section 7.3. In section 7.4 we 
show that we can associate with any x £ XAH n a nested set of subsets of {1, . . . , n}. 
This is done by analyzing the ratios. It turns out that for degenerate configurations 
certain ratios equal zero. If [3$ is close to zero, then the length of the line segment 
PiPk is very small compared to the length of the line segment PiPj. We then say 
that pi and pk coincide with respect to Pj . 

Combinations of compactifications of configuration spaces and combinatorics receive 



a lot of attention in recent years. Elegant examples can be found in [Del] and [Dc2] 



7.2 Informal introduction. 

In this section we introduce most of the terminology that we use in the rest of the 
chapter, in a more intuitive setting. 

7.2.1 Screens, clusters and nests. 

Suppose we are looking at a plane containing n > 2 sites. By sites we mean points, 
but in order to be able to use points in other contexts, we prefer the word sites. 
When we see all n sites at first glance we are happy, but when we seem to see fewer 
sites we tend to look more closely. That is, we zoom in at a site in order to see if it 
is really one site, or in fact consists of a number of sites. At the same time we will 
lose sight of other sites. We repeat this procedure until we have found all n sites. 
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Imagine that the plane is in fact a computer screen. If some sites coincide, that 
is, if they form a cluster, we click on the cluster. Another screen pops up that 
contains exactly the sites in the cluster. If we click on a site and nothing happens, 
we know that at the level we are looking at, this site does not coincide with another 
site. So every non-trivial cluster corresponds to a screen. By continuing to click 
we always find all n sites. Label the sites with labels from 1 to n and write down 
the labels that are visible in every screen. This produces a nest, a nested set of 
subsets of {1, . . . , n}, i.e. the clusters. If no sites coincide the nest consists of the 
set {1, . . . , n} only. 

The top screen is the first screen we see. It contains all sites, but maybe not all 
sites are separated. We organize things in such a way that every screen contains at 
least 2 non-coinciding sites. That is, in every screen at least two sites are separated. 
We are not interested where in the plane a screen focuses. But we do want to know 
the exact relative position of the sites. We conclude that we consider site sets up 
to scaling and translation. 



7.2.2 Pinpointing sites. 

We want to describe such a family of screens filled by sites. A naive approach is 
by listing the screens together with the coordinates of the sites that occur in each 
screen. This may describe the situation well when the sites do not move, but assume 
that two sites pi and pj that did coincide first in some screen start moving apart. 
This could change the relative positions of the sites in all screens were both pi and 
Pj occur, but it is not at all clear how a change in coordinates in one screen should 
influence coordinates in other screens. 




Figure 7.1: The hook h 1 ^ = [fi] 1 - , af- ) hinged at pi from pj to pk- The little arrow 
indicates the positive direction. 



Instead of marking coordinates of the sites involved, we mark the angle between 
any two sites pi and pj and the hook hfj between any ordered triple (pi,Pj,Pk) of 
sites. A hook h\j consists of an angle a 1 *- and a ratio . Consider Figure [zTl] . 
Suppose that the sites pi and pj, thus the leg piPj have been constructed. If wc 
rotate pj by angle afj with respect to center of rotation pi and multiply the image 
of the rotation by f3jj , we get pk- 

Using these hooks and angles we can fill a screen with sites. Suppose that for a 
set of distinct sites p\, . . . ,p n we know all hooks and angles between the sites. We 
show how to fill a screen with the sites. We put p\ in the origin, as we do not care 
about translation. As we do not care about scaling either, we put pi at distance 1 
of p\. We use the angle ol\2 to assure that the line through p\ and pi has the right 
direction, that is, we set P2 — (cos ai2, sin 0:12). All other sites pi, for i € 3, . . . ,n, 
can be constructed as in Figure 7.1 by a hook h\\ . 
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7.2.3 Degenerate configurations. 

Let CONF n be the configuration space of n distinct sites in the plane and let ^Aff 
be the map that maps an element c G CONF n to the (™) angles and the 6(") hooks 
between the sites. In fact we embed the image space in a slightly bigger space AH n , 
allowing also negative ratios. The space that we are really interested in is the closure 
of the image of CONF n in AH n . We call this space XAH n . Most points x G XAH n 
can be realized as image points of a configuration c G CONF n , that is x — V'Aff ( c ) ■ 
We call a point x G XAH n degenerate if x G V A H„ ( CONF n ) \ "0 AH„ ( CONF n ) . It 
turns out that the degenerate points are exactly those points x G XAH n such that 
at least one ratio coordinate ftf- (x) = 0. 

We go back to the nests and screens. By analyzing the set of all (3f- (x) in some 
point x G XAH n , we can associate a nest C(x) with x. If no (Hjj- (x) = 0, we 
think of x as the image of a configuration where all sites are distinct, and in this 
case we associate the nest {{1, . . . ,n}} with x. And with this nest we associate 
exactly one screen where all n sites occur and are distinct. In the other case at 
least one 8]^ (x) = 0. As we can introduce dfj; = \ Pt ~ Pk \ for distinct sites, we think 

^l] \ I 1 l] \Pi—Po\ 

of djj (x) = as if pi is very close to pk as seen from pj . But this means exactly 
that pi and pk form a cluster with respect to pj. In fact, one can prove that the set 
C(x) consisting basically of all sets CV, = {k | (3'^ ^ 0} U {i, j} defines a nest on 
{1, . . . , n}. So C(x) is the nest that corresponds to x G XAH n . For every cluster in 
C(x) with at least two elements we define a screen as a copy of the plane. 

7.2.4 Factorizing XAH n . 

Fix some x G XAH n . By analyzing the Cjj as explained above, we associate a 
family of screens with x denoted by x-screens. Our next goal is to fill the x- 
screens. The main idea of working with screens is as follows: the site set p\, . . . ,p n 
that we think of as corresponding to some degenerate configuration q G XAH n , 
should be separated completely, somewhere in the hierarchy of x-screens. This 
puts a requirement on those points q G XAH n that we can use to fill the x-screens 
with: the point q should not be more nested than x itself, in the following sense: 
if Plj (x) 7^ 0, co, then also f} 1 ^ (q) ^ 0, oo. Another way to put this: we arc not 
allowed to move down in the stratification with respect to the coincidence of sites. 
We do not need all (™) + 6(") data elements 

{aij(q))i<i<j<n,Kj (s))' k, pairwise distinct indices, 

for filling the x-screens with sites pi(q), ■ ■ ■ ,p n {q), given some suitable q G XAH n . 
We just need enough data elements to pinpoint each site once. By exploiting the 
nested structure of x, we can avoid bad choices in picking the data elements. We do 
not want to be confronted with legs of length • 00 for example. These prerequisites 
result in a factor Dom„(x) of data elements of XAH n that are both suited and 
enough to fill the x-screens, given some suitable q G XAH n . Note that the choice 
of the data elements depends on x, although we use of course the data elements of 
q itself in order to fill the x-screens with sites pi(q), . ■ ■ ,p n (q)- 

7.2.5 The hooked tree. 

Again by using the nested structure of our fixed x G XAH n , we can define a partial 
order on the labels 1, . . . , n. This partial order tells in which order exactly we have 
to add the g-sites pi(q), . . . ,p n {q) in the x-screens. It turns out that with every 
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label i e 1, . . . , n we can associate an x-tag l x {pi). Such an x-tag is a pointer to 
one coordinate of Dom„(i) C XAH n . As any such coordinate is cither an angle 
or a hook, we conclude that every label i corresponds to a geometric construction 
(construct Pi as the end of a line segment of a given angle, or construct as the 
end point of a given hook) that tells how to construct pi . 




1 2 6 3 5 



Figure 7.2: The hooked tree for the nest < {1, 2, 6}, {3, 5} >. 

We store the set of x-clusters in a tree. The vertices of the tree are exactly the 
x-clusters. Two x-clusters are connected by an edge if one cluster is maximal in the 
other. Some edges are labeled by x-tags. This is done in such a way that if going 
up in the tree from leaf i, the first x-tag encountered points exactly to that hook 
or angle that is used to construct pi in the set of x-screens. This tree is called the 
hooked tree. 

7.2.6 Filling screens. 

We discuss the use of the partial order a bit further. Our setup is such that in 
each x-screen separation occurs. Fix some x-screen S. As mentioned above, any 
x-screen S corresponds to some x-cluster Cs. Assume that all x-screens T, such 
that the cluster Ct contains Cs, have already been filled. We say that such screens 
T are above S. The x-screen S should exactly contain the sites with labels in 
the x-cluster Cs after filling. Moreover, the sites corresponding to the maximal 
subclustcrs of S are separated in the x-screen S. Let C\ be the maximal subcluster 
of S with smallest minimal label i and let Ci be the maximal subcluster with second 
smallest minimal label. The first site we put in S is the site pi, that we place at the 
origin. The second site to construct is the site pj, that we place at distance 1 off),. 
The angle that the line segment piPj makes with the positive axis is the so-called 
screen orientation Os of S. It is defined recursively in terms of the screens above 
S. Note that by now we have constructed one leg PiPj in S. All other sites in S 
can be constructed using hooks. 

Suppose, for instance, that we want to construct all sites in some maximal subcluster 
Cfe 7^ C\, C2, where k is the minimal label of Ck and therefore k > i,j. We construct 
Pk by means of a hook hfj . One can think of this hook as a hook on scale, as 
Pi,Pj and pu are separated at the same level. Note that k is automatically minimal 
in the maximal subcluster of Ck that has smallest minimal label. Let m e Ck, 
with m > fc, be minimal in the second maximal subcluster of Ck- We construct p m 
using /ij:™ . This is a different type of hook as pk and p m are separated one step 
deeper in the nest structure. Think of this type of hook as a explosion hook: 
for x itself |/3£™ (x)\ — 0, as pk and p m 'coincide with respect to pi '. If we take 
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an arbitrary point q £ Dom„, then \f3%™ (g)| > 0, and if (q)\ > we could 
say that an explosion occurs, as suddenly p^ and p m do not coincide anymore with 
respect to Pi. 



7.2.7 XAH n is locally the graph of a function. 

One purpose of this setup is to minimize the dimension of the factor Dom„(i), 
that is, we want to use as few data elements of XAH n for filling the x-screens as 



possible. That we cannot do better is shown by Lemma 7.31 that proves that the 
dimension of the factor Dom„(i) equals the dimension of the reduced configuration 
space conf„. Given a set of filled x-screens, we can just read off angles and hooks 



from the screens. Theorem [7.58 proves that filling screens using data elements is 



consistent with reading off data elements, in the sense that we get back the data 
elements that we started with. Note that besides the data elements in the Dom„ (x) 
factor, we can also read off the data elements in the complementary factor Rng„(x) 
of Dom„(i) in the space of all angles and hooks on n points, AH n . This shows that 
XAH n can be written as the graph of some funct ion read : Dom„(x) — > Rng n (x). 



As a consequence, we can prove in Theorem |7.61| that XAH n is a smooth manifold. 
In fact there are some extra requirements on the structure presented to obta in 
smoothness that we did not discuss here. Details can be found in Sections |7.3| to 



7.8 



7.2.8 Connection with Kontsevich-Soibelman. 



In [KS], see also Section 7.1.2 , Kontsevich and Soibclman describe a manifold with 
corners FM2(n) that is closely rela ted to our smooth manifold XAH n . The exact 
relation is given in Theorem 7.65 that describes a map / : FM2{n) — > XAH n 
together with the fibers of /. Due to this close relation we can use the construction 
of the filled x-screens for points x <E FM2(n) as well. This is important in the 
following application. For an overview of the compactifications introduced so far, 



consult Table 7.1 



angles 
mod 27r 



angles 
mod 7r 



clickable 
Voronoi 
diagrams 



FM 2 {n) 



XAHn 



angles 

and 

hooks 



Voronoi 
diagrams 



CDA r . 



CUAn 



angles 

and 

points 



Table 7.1: Overview of compactifications. 



7.2.9 Adding Voronoi diagrams in the x-screens. 

We apply the construction of filled x-screens to the analysis of Voronoi diagrams 
of points sets that may include coinciding points. We proceed by associating a 
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Voronoi diagram to any point x £ FM2(n). First, consider some non-degenerate 
point x £ FM 2 (n). Then by definition, as in the XAH n case, x can be expressed as 
image of some configuration c £ CONF n . As no pf? (x) = for a non-degenerate 
x £ FM 2 (n), there is exactly one x-screen associated with x, the top screen, that 
contains all sites pi(x), . . . ,p n (x) in such a way that all sites are distinct. The 
Voronoi diagram Vfm (x) associated to x in this case is the Voronoi diagram of the 
n sites in the top screen. 

Next assume that a; is a degenerate element of FM 2 (n) . Suppose that we have filled 
the x-screens with the x-sites Pi(x), . . . ,p n (%) as described above. In this case, there 
is more then one x-screen, but it is guaranteed that in any x-screen S exactly those 
clusters are separated that are maximal subclusters of the cluster Cs ■ We define a 
Voronoi diagram Vfm{x) in terms of the x-filled x-screens in two steps. In the first 
step we add the classic Voronoi diagram of the maximal subclusters of Cs to any 
x-screen S. This is called the initialization step. In the so-called completion step we 
recursively add the completed diagrams of all screens below S to S for all x-screens 
S. This is done as follows. A non-trivial cluster is a cluster that consists of at least 
two sites. Any non-trivial maximal cluster M of S corresponds to an x-screen Sm 
that is directly below S. We paste the completed Voronoi diagram Vfm(Sm) of 
Sm into the Voronoi cell of M in the initialized Voronoi diagram V(S). 

This gives us a clickable Voronoi diagram Vfm(x) for any x £ FM 2 . By clicking 
on a non-trivial maximal cluster, if present, a screen appears. This screen contains 
the sites in the cluster and the Voronoi diagram of the sites. Recall the polynomial 
sites from Chapter ||. By computing ratios and angles we can associate an element 
xs £ FM 2 to a set S(t) of polynomial sites. Then the completed Voronoi diagram 
Vfm(T) of the top screen T for xs matches the limit Voronoi diagram of S(t) at 
t = 0. 



7.3 Angles and hooks. 

We start this section by recalling and introducing data elements defined on pairs 
and triples of distinct points in M. 2 . These data elements form the basis of the 
compactification of a suitable quotient of the configuration space of n distinct points 
in the plane that we define later on in this section. 

7.3.1 Data elements of pairs and triples of points. 

Given two distinct points pi and pj , we can determine the angle between the line 
through pi and pj and the positive x-axis. We have done this in Definition |5.1| . 
There we have introduced the directed angle o>ij £ R/27rZ and the undirected angle 
oH £ R/7rZ. We store information on the geometry of three distinct points in hooks, 



see Figure 7.3 



Definition 7.2 

Let pi , pj , and pt be three distinct points in the plane. The angle af- from pj to 
Pk hinged at pi is given by a\j := oiik — ocij mod 2n. The ratio (3fj from pj to pu 
hinged at pi is given by 

oik . \Pi ~ Pk\ r,ik ~ in „ \ 

Pv := I — z — T' % e(0>°°)- 

\Pi Pj | 

The hook hfj from Pj to Pk hinged at Pi equals the ratio together with the angle: 
hfj (ftlj ' a ij )■ The point pi is the hinge point of the hook hf? while the line 
segments PiPj and piPk form its legs . 
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Figure 7.3: The angle a\j from pj to pk hinged at pi. The small arrow indicates 
the positive direction. 



Note that one triple of distinct points gives six hooks. 
Remark 7.3 

An alternative is to consider hooks in the complex plane. Let pi 1 pj, and Pk be 
distinct points in C. Note that cnj = (j>i — Pj)/\pi — Pj\- It is easy to check that, 
with / = V - 1) it holds that 

P^_Pk = «,/«{?. 
Pi ~ Pj 

7.3.2 Geometric interpretation of the hook h\j . 
Notation 7.4 

Let p and q be points in the plane and a £ [— 7T, n) an angle. By rot a ^ p (q) denote 
the image of q under the anti-clockwise rotation around p by angle a. We abbreviate 
rot a! ( ,Q) by rot Q . 

Suppose we are given two distinct points pi and pj , together with a hook h 1 ^ . The 
hook h\j can be seen as a prescription that tells how to construct the point Pk- 
Using afj , we first construct p'p the image of the rotation rot Q ifc p . , applied to pf 



rot 



The half- line that starts in pi and passes through p', contains pk- The ratio /3. 



fixes the distance from pk to pi, thereby fixing pk itself: 

Pk = 0$ ip'j ~ Pi) + Pi- 
Summary: we interpret a hook hfj as a point rotation followed by a vector mul- 
tiplication. As a result the leg piPj is transformed into the leg PiPk ■ 

7.3.3 Negative ratios. 



In this section we give sense to 'negative ratios', as illustrated by Figure [7^. Think 
of the line through the points pi and p'j as a local coordinate axis Ay/ . The origin 
on this axis is pi and the point with coordinate 1 is p'y The position of pk on 
this axis is given by the coordinate /J^ . In this setting, any value of /3j^ in the 
interval (— oo, oo) is meaningful: it just indicates the point with coordinate (3f^ on 
the Ay/-axis. Suppose we allow the full interval (—00,00) for values of (3fj . Then 
there are two different ways to construct a point pk by a hook h] k , given the points 
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Figure 7.4: The local coordinate axis Aiy . 



Pi and pj. The hooks (/^j 5 , ) and (— , af- + n) both result in the same point 
Pk- If identifying — oo with oo, we have by now described a map 

<j> ijk ■ CONF 3 {R 2 ) -> P 1 x R/2ttZ 



(pi,Pj,Pk) (^h,Om), 



under the identification 



+*r). (7.1) 

Recall that CCWf^R 2 ) denotes the configuration space of three distinct points in 
the plane. P 1 denotes the projective line. It can be defined as S 1 with antipodal 
points identified. A point on the projective line can be seen as the coordinate q of 
a third point on an axis with respect to two fixed points ^ 1. A homeomorphism 
from the circl e to the line E U {oo} given by stereographic projection is shown 



in Figure 7.5 A line through P intersects the circle in another point, Q, and hits 
the x-axis in a point R. When P and Q coincide, this line is horizontal, so P is 
mapped to R = oo, cf. also (rJ. 







\^R 









Figure 7.5: Stereographic projection yields a homeomorphism between the circle S 1 
and the projective line KP 1 = RUoo. 



7.3.4 Klein bottle Kjf . 



a 

b b 

a 



Figure 7.6: The Klein bottle is obtained by identifying sides a and b according to 
the arrows. 



The Klein bottle is the quotient space obtained from a rectangle by identifying 
opposite sides, see Figure 7.6; see also [Mu|. In Figure 7.7 we have depicted another 
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rectangle, but this time we add a description of the axes. For the horizontal axis 
we take the projective line P 1 . The vertical axis is ]R/27rZ. Recall the identification 

i ) 



(Mi 



(—Plf ,a]^ + 7r) introduced in 7.1 



Lemma 7.5 

The quotient space P 1 x M/2ttZ / 



■>k defines a Klein bottle. 



Proof. Consider Figure [7.7 
[0,oo 

(i) 
(ii) 
(hi) 



Any class in P x R/2-7rZ/ ~& has a representative in 
x [0, 27r]. On [0,oo] x [0, 2tt], there are three identifications. 
The periodicity of directed angles in ~R/2irZ is indicated by 3>. 
The identification (0, a\j ) ^k (0, af- + tt) is indicated by A. 
If (0, a$ ) ~ k (0, af- + tt) holds on the compact set P 1 x M/2ttZ/ ~ fe , then 

This gives 



also (oo, a\j ) 

the third identification, indicated in Figure 7.7 by A 



(oo, a|y + 7r), as any flf? = corresponds to /3, 



I' or identifying the quotient space we are allowed to cut along a*^ = tt, if we 
eventually paste back again along the same cut. This cut is indicated by oo We do 
cut and paste the bottom half on the top half, applying A and A. This results in 
the cylinder shown on the right of Figure 7.7. On the top of the cylinder, 3> and 
oc are directed clockwise, while on the bottom they are directed counterclockwise. 
This shows that we are in the situation of Figure 7.6. ■ 



T4- 



2tt 



9- - - - 

L 


\ \ A 
i-ee 


L 


* 

X A 
» 




Figure 7.7: The identification (/?- fc , 0$ ] 



What is essential for us is that the Klein bottle is a smooth manifold. Note that 
with every ordered triple k) we have associated a Klein bottle Kf- as the target 
space of 4>ijk- We can write 4>ijk as a composition: <f>iju — Cijk ° h^k- The map 
hijk maps an ordered triple of distinct points to its hook {(if- , aiy ), where > 0. 
The second map, c^, maps a hook (ffij , a * ) to its class [{(3$ , a\j )] under the 
equivalence relation We have: 

<f> ijk : (K 2 ) 3 \A„- fc H fc [0, 00] x K/2ttZ ^ ([-oo s oo]xK/27rZ)/~ fc) 

Proposition 7.6 

TTie map 

: (ho - 00 ! x P/2ttZ)/ - fe -> [0,oo] x K/2ttZ, 
is 2 to 1 if @lj equals or ±00 and 1 to 1 eZse. 
Proof. The inverse images of Cjjfc are as follows: 
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(i) cr.i( [(0, «{* )] ) = {(0, af) ), (0, a {* + tt)}; 

(ii) cr.i ( [(oo, o% )] ) = {(ex., ag ), (oo, a# + tt)}; 

(iii) %I( «« )1 ) = {(°°> )' (°°> + «■)}; 

(iv) cr.i( , )] ) = (/$, ag ), if < $j= < oo; 

(v) cr.i( , )] ) = (-0$ , + tt), if -oo < < 0. 

Remark 7.7 



Consider the image on the right in Figure 7.7. If we cut along 0$ = 0, that is, 



along A, we obtain a Mobius strip with /3y = on the boundary. Cutting along 
= oo, that is, along A, leaves us with two pieces joined along oc. This results 
in a cylinder with (3^ = on one boundary component and fSfj — oo on the other 
component. 



7.3.5 The data map from conf n to AH n . 

In this section we describe a map from a set of distinct points in the plane to all 
ang les and hooks between pairs and triples of those points. Recall that in Definition 



5.4 



we have introduced the reduced configuration space conf fi as the space of n 
distinct points in the plane up to scaling and translations. For a representative c of 
a class [c] G conf„ , write the following data: for every pair of points Pi,Pj S c,i ^ j, 
the undirected angle ay = ciji ; for every triple Pi , pj , pk of points in c, the six hooks 
h-lj ! h-li , h j ji > hfk i ^fcj i an d ■ This gives (™) unordered angles and 6Q) hooks. 
Note that the angles and hooks that we obtain are independent of the choice of a 
representative of a class c G conf„ . 

Definition 7.8 

AH n is the space of hooks and angles on n points: 

AH n := (M/ttZ)(5) x ( P 1 x (R/2ttZ) / ~ fc ) 6 (a) , 



where denotes the identification defined in 7.1 



Remark 7.9 

Being a direct product of circles and Klein bottles, AH n is smooth. 

In the next definition, we introduce a compactification of the reduced configura- 
tion space conf„. 

Definition 7.10 

The data map ipj^jj is the map 

^AH n '■ conf„ ->■ AH n , 

4>AH n '■ conf « (K/ttZ)(*) x ( F 1 x (R/27rZ)/~ fc ) 6 ^, 

, a|| ) | i, j, k pairwise distinct), 

and XAH n is the closure of the image of conf„ in AH n . 

Alternatively we could allow non-negative ratios only, construct a data map from 
conf n to (R/7rZ)(2) x ( [0,oo] x (R/2?rZ) ) 6 (*\ and define XAH n as the closure 
of conf„ in this product. In this approach, however, 'jumps' occur in every hook 
h\j . Indeed, suppose that pj ^ p^ are fixed while Pi moves 'through' p^. Then the 
ratio stays close to zero, but a]j jumps to af- + n at collision, as changes 
direction, while onj stays the same. 



7.4. Nests and screens. 



99 



7.3.6 Example: a smooth compactification of conf 3 . 

Consider three points p±, p 2 , andp 3 in the plane. Their configuration space is given 
by conf 3 = {{pi,P2,Ps) \ Pi ^ P2,Pi ^ Pz,Pi £ P3}- Then conf 3 is mapped by 
^AH ^° a ra ther big space: 

^AH n : COnf 3 ~> S 12 X S 13 X S 23 X K 12 X K 13 X K 21 X #23 X K 31 X #32 ■ 

Indeed, we can make a smooth compactification of conf3 by mapping it to a much 
smaller space: 

Lemma 7.11 

Let -0 be the map ip : conf3 — > Si 2 x .Kff J™™ &2/ (pijP2,P3) 1— ► (0:12, TTien 
^[confa] = 5i2 x K J§ ■ Moreover, S12 x i^if is smooth. 

Proof. We have to show that for every point y G Si 2 x an d every neigh- 
bourhood U of y it holds that f n (/i[con^] 7^ 0. Any point y G S12 x has 
coordinates: 

2/ = (a 12 (y),Pll(y),alUy)). 

If < l/3il(y)l < °°, then (a 12 (y),(3l§(y),a{l(y)) = V(pi,p 2 ,p 3 ), where p x = 
(0, 0), p2 = (cosai2,sinai 2 ), and p 3 = /?ff (cos(ai 2 + a\l ),sin(ai 2 +a\l )). If ratio 
012 (v) = 0) ±00, then in any neighbourhood of y there exists a point z such that 
< \0\l (z)\ < 00, so z G V>[conf 3 ]. ■ 

Later on we prove that XAH[3] is another smooth compactification of conf 3 . Map- 
ping of conf 3 into AH n has the advantage that there are no preferred labels. 



7.4 Nests and screens. 



The main result of this section is Theorem [7.21| that associates a nested set of 
subsets of {1, ... ,n} to an element x G XAH n . This is done by analyzing those 
ratios fif- (x) that are equal to zero. 



7.4.1 Clusters, nests, and screens. 

Definition 7.12 

Let Z be a collection of sets. 

(i) Z is nested iff either A <Z B or B C A or An B = for any two elements 
A.BeZ. 

(ii) A set U is a maximal subset of a set W, w.r.t. some nested set Z iff 
U,W G Z, U CW and there exists no V G Z with U CV CW. 

(Hi) A nest Z on {1 ... n} is a nested set of subsets of {1, . . . , n} that includes 
the set {1, . . . , n} itself and all singleton sets {i}, where i = 1, . . . , n. 

(iv) The elements of a nest are called clusters. The singleton sets are called the 
trivial clusters. 

Example 7.13 

The set Z = {{1, 2, 3, 4, 5, 6}, {1, 2, 6}, {3, 5}, {1}, {2}, {3}, {4}, {5}, {6}} is a nest 
on {1 ... 6}. For simplicity, we often write only the non-obvious elements of Z. The 
obvious elements are the singleton sets and the total set {1, . . . ,n}. We indicate 
this by using the (, ) brackets. In our example we would write Z = ({1, 2, 6}, {3, 5}). 
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Let Z(n) be a nest. Associate a copy Sc of the plane to every non-trivial cluster 
C € Z{n). Every copy is labeled by its corresponding cluster. In this way we get a 
family of screens Sz = {Sc \ C E Z, C non-trivial}, associated to the nest Z. In 
case there is only one non-trivial cluster, we have exactly one screen, iSr„i, where 
stands for {1, ... , n}. In Example 7.13| we have screens 5[6], <Si,2,6j, and 63,5. 



n 



In this section we show how to associate a nest Z to a point x € XAH[n]. In the 
following sections we describe how to fill the family of screens Sz corresponding to 
a nest Z in a meaningful way. 

7.4.2 Properties of ratios. 
Lemma 7.14 

Let Pi, pj, pk, and pi be distinct points in the plane. Then the following relations 
hold: 



(i) 


Pij 


e (0,00); 


(ii) 


pt 


= V/f; 


(in) 


$ 


+ /f >i; 


(iv) 


pi 


oik oik 
' Pij — Pil 


(v) 


Pt 


okl oil 
' Pki — Pij 



■Pli- 

Proof. 

(i) Follows directly from the definition. 

<H» >/t - fcff - «■ 

(iii) This is in fact the triangle inequality: 



pik + 0jk = \Pi ~Pk\ + \Pj -Pk\ = \Pi~Pk\ + \Pj ~Pk\ > 1 



\Pi~Pj\ \Pi~Pj\ \Pi~Pj\ 



f iv \ rH _ \Pi~Pj\ \Pi~Pk\ _ \Pi~Pk\ _ Rlk 

( j PU ~ \ PI -Pl\ ' \ PI -P 3 \ ~ \P*-Pl\ ~ ^ l ' 

(v) Apply the definition. 
Corollary 7.15 

For x G XAH[n] and i, j , k, I distinct labels we have 
(1) \(3^\ € [0,oo]; 

(a) \P%\ = VI^I; 

(in) I/^I + I^I>1; 

(to) 7/max( |/$ |, |/f |) < 00, then \ ■ \/®\ = \. 



Proof. The claims follow from Lemma 7.14 by taking limits, but notice that for 
x E XAH n some /3|y can be negative, see Definition 7.8. ■ 



7.4.3 Separating clusters. 
Remark 7.16 

For the rest of this section, we identify a site pi with its label i. 
Set 

Ca := {Pk I |/$|>0}U{i,j}. (7.2) 

Intuitively, think of the following. Recall the geometric interpretation of a hook 
h % l k given in Section 7.3.2. Suppose we use hf? to draw a point pj in the plane, 
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given the points pi and p k . The ratio \0 1 ^. | gives the distance between the points 
Pi and Pj , compared to the unit distance piPk ■ Imagine an observer situated at p k . 
When |/3*^ | gets very small, pj moves very close to pi. But as long as \(3^ k | > 0, 
the observer at p k still can distinguish pi and Pj. Therefore one could say in this 
setting that Cy consists of all points p k such that pi ^ pj w.r.t. p k . We call Cy the 
separating cluster of Pi and pj . 



Pj 




Pr J Pk 



Figure 7.8: The triangle PiPjPk ■ 



Lemma 7.17 

We have Cij = Cji . 



Proof. Let us start with a heuristic proof: consider Figure 7.8. Note that / = 
\ Pk -pj\,J=\p k -Pi\ and K = \pj-pi\. Therefore \0* | = K/J and \(3% \ = K/I. 
We have to show that K/J > implies that K/I > as well. K/J > only 
excludes that Pi and pj coincide w.r.t. p k . But this is also the only way to make 
K/I small. 

More formally: we have to prove that \(3 l ik \ > implies that |/3j^, | > 0. First of all: 
as \P% | > it follows that \/3j!f \ < oo. Now suppose that \f3f k | = 0. Then also 
IPjl I ' \$j I = \Pkj I = 0- Tnis contradicts the triangle inequality applied to |/3j£j \ 
and \(3f k \. ■ 

Lemma 7.18 

IfPk G Cy and pi Cy, then p x C ik ,C jk . 

Proof. First, p k S Cy implies that |/3j| | < oo. Secondly, as pi $ C\j it follows 
that |=0. So both | and |/3y | are finite, which means that we can apply 



Corollary [7.15| and get \(3f\ = \P$\\P£\ = 



■ik I _ I oij 1 1 oifc 



Lemma 7.19 

Lef i,j and fc &e distinct labels. Then 



Pk £ Cy => Cife C , 
Pk 2" Cy => Cy C Cjj.. 



Proof. If pfc 6 Cy , then Cjfc C Cy by Lemma 7.18. We are left with the case 
Pk & Cij. We prove that in this case Cij c Cjfe. For this purpose, we fix p\ G 
Cy and show that S Cjfe. Since ^ Cy, it follows that |/3^ | = 0. Due to 
G Cy, one has |/3^ | > 0, and in particular \(3f- | is finite. By Coroll 
I = \P% I ' l/3y I = 0. As a consequence |/3|f | > 0, so indeed p, e C 4fc 



7.15 



Lemma 7.20 

IfPa,Pb € Cy, #1671 C ah C Cy. 
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Proof. If pi, G Cia, it follows from Lemma 7.19 that C a b C C\b C Cij. Suppose 
that pb £ C ia . Then Lemma 7.19 gives C ia C Cu,. This implies that p a G C^. 
Therefore, C a h C Cju C C,,-. ■ 



7.4.4 The clusters form a nest. 
Theorem 7.21 

Fix x G XAH n . Let C(x) := {Cy ;}i<i<j< n U {{1}, ■ ■ ■ , {"-}}■ TTien C(x) is a nest 
on n. 



Proof. We have to show that for any clusters C s t an d Cy either C s t n Cy = or 
one of them is contained in the other. Lemma 7.19 deals with the case s = i: it 
shows that the nest condition is fulfilled for Cu and Cy . When C s t n Cy = the 
nest condition is fulfilled as well. So, assume that s, t , i, and j are all distinct and 
that p m G C st n C, 
hold: 



y . We have proved in Lemma 7.19 that the following inclusions 





-*- 


C sm 


n 




n 


c^ 






u 




u 


Cjjn 


* 


Ctm 


open. There 


are 


1 


2 


3 


* C 


D 


c 


* c 


D 





4 

D 
C 



In the first case, and pj are in C s j. So, Cy G C s t by Lemma 7.20. In the second 
case, C st C Cu. Concentrate on the case 3 now. 



C%ra G C s 
Cjm Z> Cj 



I ni 



Pi € C st , 
Pt G Cy. 



If Pj G C sm , then G C s t, so both ^ and pj are in C s t. This implies that 



Cy C C s t ■ When ^ C sm , we can apply Lemma 7.18 as pi G C sm and conclude 



that Cis C Cij. But this implies that p s G Cy. As pt G Cy as well, it follows that 
C s t C Cij. The fourth case follows from the third by relabeling. ■ 



From Theorem 7.21 we obtain an important corollary. 



Corollary 7.22 

If x € XAH[n], then C = C{x) forms a nest. 



For x G XAHn, we call a cluster C G C(x) an x-cluster. By now, the following 
statements are obvious. 

Corollary 7.23 

Let C be a nest. Then 

(i) Cij is the smallest cluster of C that contains both i and j; 
(ii) i and j belong to distinct clusters in the screen Sc ir 

Definition 7.24 

Fix x G XAH n . 
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(i) Let M = . . . , i m }, for m in 2, . . . ,n, be a subset of {1, . . . , n}. The 
separating cluster Cm of M in x is the smallest cluster of C(x) that contains 
,i m . The separating screen Sm is the unique screen that contains exactly 
the elements of M . 

(ii) Let C(x) be as in Theorem 7.2\ . The x-screens are the screens Sc for 

c ec(x). 

Suppose we are given a point x G XAH n . By determining the clusters CV; for x, we 
find the nest C{x) for x. From the nest, we get a list of screens Sc, the x-screens. In 
this section we introduce a data structure that enables us to fill the x-screens, the 
hooked tree. The hooked tree ht(x) encodes both the nest C(x) of x and the relative 
positions of the underlying point set px, . . . ,p n of x. We define ht(x) recursively. 



7.4.5 Definition of hooked tree. 



Definition 7.25 

The hooked tree ht(x) of a point x G XAH n is a rooted tree. Some of its edges 
have tags. The vertices of ht{x) are the clusters of the nest C(x). The root is 
the cluster {1, . . . , n}. The leaves are the clusters {1}, . . . , {n}. Two clusters are 
connected by an edge if one of the clusters is maximal in the other. 

The x-tags refer to hooks and angles in AH n . Fix a cluster S. Order the maximal 
subclusters in S according to the smallest label in each maximal subcluster. Suppose 
S has maximal subclusters c%, . . . , Cfc with minimal labels 

h < h < ■■■ < 3k- 

If S is not the top cluster [n] , there exists a cluster T that contains S as one of its 
maximal subclusters. Suppose that T its maximal subclusters have minimal labels 

h < h < ••• < i m - 

Tag the edges from cluster S to its maximal subclusters as follows. 

Edge to Cj. If S is the top cluster [n] then the edge from S to c\ is tagged top. 

Else: no tag. 
Edge to C2- 

(i) If 5* is the top cluster [n] then the edge to C2 is tagged ct\j 2 . Such tag 
is called an x-type 2. a tag. 

(ii) Else if ji > i\ then tag the edge from S to C2 by h 3 . 1 ^. This is an 
x-type 2.b tag. 

(Hi) Else if ji = i\ then tag the edge to C2 by h^j^. This is an x-type 2.c 
tag. 

Edge to Ct, for t > 3. The edge from S to ct, for t > 3, is tagged hj^j 2 . This is 
an x-type 3 tag. 

The tag oiij in x-tags does not stand for the value aij(x), but rather for a map 
AH n — > R/7rZ. The indices i and j indicate which particular factor of the form 
R/irZ in AH n is meant. 

Definition 7.26 

Fix x G XAH„ . Let q be an arbitrary point in AH n . A q-coordinate hfj (q), respec- 
tively, ctij(q), afj(q), fif? (q) is a x-type 2 hook, respectively, angle, angle, ratio 
iff h\j , otij, a\j or (3^ occurs in x-tags; x-type 3 hooks, angles, and ratios are 
defined in a similar way. 
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7.4.6 The hooked path. 

In this section we associate a unique x-tag to every site pi. Moreover, we show that 
the hooked tree ht(x) induces a partial ordering on the labels 1, . . . , n. 

Definition 7.27 

Fix x € XAH n . 

(i) To any site pi, for i = l,...,n we associate an x-tag l x {pi): it is the first 
x-tag that one encounters when going up in ht(x), starting from the leaf i. 

(ii) The predecessor p x {Pi) of pi with respect to x, for i = 2, . . . ,n is defined as 
follows: 

1. p\ has no predecessor. 



2. If l x {pi) is an x-type 2 tag then i ~ j2 in Definition 7.25 

• In case of x-type 2. a, p x {Pi) = Pi- 

• In case of x-type 2.b, p x (pi) = Pj 1 - 

• In case of x-type 2.c, p x (pi) =Pi 2 - 



3. Ifl x (pi) is an x-type 3 tag then i — j t in Definition 7.21 and p x (pf) = pj 2 . 
(Hi) To a site pi, for i = 1, . . . , n we associate a hooked path L x {pi): 

1. L x ( Pl ) = (). 

2. L x (j>i) equals L x (p x (pi)) with pi prepended. 

If a site pj occurs in L x {pi), then we say that pj is a; -above pi. If Pj = p x (j>i) or 
Pi = Px{Pj)> then we call piPj an x-leg. 

Lemma 7.28 

The predecessor relation gives a partial order on the set pi, . . . ,p n . 

Proof. We have to show that there are no loops. For this purpose we define the 
depth of a site pi as the number of edges from the root in ht(x) to the node 
where the edge that has l x {pi) as its x-tag starts. Consider site pi. Site p\ has 
no predecessor. Assume that i ^ 1 and that pj is its predecessor: pj — p x (pi). 
We claim that either depth(pj) < depth(pj), or depth(pj) = depth(p,) and j < i. 
Indeed, if l x (pi) is of x-type 3, then the last case holds. And if l x (pi) of x-type 2, 
then depth (pi) = depth (pj) + 1. ■ 

Lemma 7.29 

There is a I to 1 correspondence between the sites pi, . . . ,p n and the x-tags. 



Proof. In Definition 7.27 we have associated an x-tag with every site pc the x-tag 
l x (Pi) associated with pi is the first tag one encounters going up from the leaf pi 
in the hooked tree ht(x). To finish the proof, we turn around the association: wc 
show that for every x-tag I there is a unique path to a leaf of the hooked tree such 
that no other x-tag is encountered. If the edge tagged I is itself incident with a leaf 
we are done. So assume that the edge tagged I does not end in a leaf. Then it ends 
in a vertex of the hooked tree corresponding to a cluster S of at least two elements. 



By Definition 7.25 only the edge from S to its maximal subcluster M with minimal 
label ji has no tag. j\ will also be the minimal label of every subcluster of M that 
contains j\ . Therefore there is a unique path without labels to pj ± . ■ 



7.4.7 Dom n (i) and Rng n (x). 



The x-tags define a factor of XAH, 
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Definition 7.30 

Fix x € XAH n . 

(i) The factor Dom„(i) is the factor of AH n that corresponds to the angles 
and hooks in x-tags. The factor of AH n complementary to Dom„(i) is denoted by 
Rng„(x). 

(ii) By tv x , denote the projection of AH n onto Dom„(i); 'tv' stands for tree 
values. The projection of AH n onto Rng„(x) is denoted by ctv x , the complemen- 
tary tree values. 

The following lemma states that dim(Dom„(x) ) = dim(conf„). This shows that 
the hooked tree is optimal in the sense that it refers to as many coordinates of AH n 
as the dimension requires. 

Lemma 7.31 

Let x e XAH n . Then dim( Dom„(i) ) = 2n - 3. 

Proof. The dimension of Dom„(i) is given by the sum of the degrees of freedom 
in the hooks and angles referred to by x-tags. Every hook is defined by one angle 
alj and one ratio , so the dimension of one hook equals 2. The dimension of 



a single angle <xy equals 1. The dimension of 'top' is zero. From Lemma 7.29 it 
follows that the x-tags refer to n — 2 hooks and 1 angle as there are n sites and one 
'top' tag. ■ 



7.4.8 Example of a hooked tree. 
Example 7.32 



Let us consider again, as in Example 7.13 , the nest Z — ({1, 2, 6}, {3, 5}). A re 



alization of Z consists of the three filled screens S^, £1,2,6 and 6*3,5. Recall that 
[6] = {1, . . . , 6}. A possible realizati on is given in Figure |7.9| . The hooked tree 



corresponding to Z is given in Figure 7.10 . Note that the dimension of Dom„(x) 
is indeed 9. 



The first column of Table |7.2| gives all sites in Z. The x-tags corresponding to each 
site are shown in the second column. The third column gives for every x-tag its 
x-type. The fourth column give s tv x {x) coordinates corresponding to the way the 



x-screens are filled in Figure 7.£ . How this filling is done exactly, and the meaning 
of the fifth and sixth column will be explained later in this chapter. 

From the data for the tv x (x) coordinates it follows, for example, that the length of 
the leg 12 is zero with respect to the length of the leg 13. And as the length of the 
leg 16 is of the same order as the length of the leg 12, it must vanish with respect 
to the length of the leg 13 too. Continuing in this way reveals the complete nest 
structure. 



7.5 Filling x-screens. 

In this section we present a method for filling x-screens, given a suitable point 
q S Dom„(x). 

7.5.1 Standard form of x E XAH n . 

First we put a fixed x € XAH n in standard form: if necessary we change some 
representatives in Klein bottles Kf? and/or change the representative for the top 
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Figure 7.9: The screens for the nest ({1, 2, 6}, {3, 5}). 




Figure 7.10: The hooked tree for the nest ({1, 2, 6}, {3, 5}). 



angle ctij 2 (x) modulo ir. We only use the Dom„(i) factor of AH n in order to fill 
the x-screens. Therefore we only have to choose standard representatives for those 
angles and hooks that are referred to in x-tags. 

Definition 7.33 

x E XAH n is represented in standard form iff 

• ctij, 0$ of x -type 2 a\j{x), af? (x) £ [— ir/2,n/2) and 

• ofx-type 3 $j {x) > 0. 

By the following, we do not have to make a sign choice for any x-type 2 ratio j3\j (x). 

Lemma 7.34 

Fix x € XAH n . 

(i) Ifhjlll^ is an x-type 2 hook, then 0fj* 2 (x) = 0. 

(ii) If h-j 1 ^ is an x-type 3 hook, then < (x)\ < oo. 



Proof. Follow the notation of Definition 7.25| . We give a proof for the x-type 2.b 
case. The x-type 2.c case is similar and is omitted. From Definition 7.25 we know 



that there exists an x-screen S such that pj 1 , pj 2 <= S but pi 1 ^ S. The claim follows 
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site x-tag x- type tv x (x) q G Dom„(i) representatives 

1 top top 

3 a 1)3 2. a -7° 69° same 

2 G9j|,al§) 2.c (0,71°) (.35,-91°) (-.35,89°) 

4 (#3, ax!) 3 (1.5,50°) (-1.92,4°) same 

5 (^Ii 5 ,"!!) 2.b (0,44°) (-.34,170°) (.34,-10°) 

6 (P12 , o$) 3 (1,-37°) (.91,-11°) same 



Table 7.2: x-tags and values. 



from the definition of separating cluster, see in particular Equation |^2|. If /3j*j* is 
of x-type 3, then pj 2 e Cj 1 j t , so (x)| > 0. Suppose that {%)\ = 00. Then 
— 0, but this would imply that pj 2 ^ Cj 1 j t . That is a contradiction. ■ 



Changing an x-type 3 ratio /3|y (x) from negative to positive involves of course a 
swap a\j (x) = a\j (x) + it for the corresponding x-type 3 angle: we have to oblige 



the Klein bottle identification 7.1 



Remark 7.35 

We also refer to x-type 3 hooks, angles and ratios as hooks on scale, etc. The 
unique x-type 2. a angle is referred to as top angle. Hooks, angles and ratios of 
x-type 2.b or 2.c are also called explosion hooks, etc. 



An explanation for these names can be found in the geometric interpretation of 



a hook h\j , given in Section |7.3.2[ Think of what changes in the hook h\j as 



we start from the fixed value x € XAH n and move away from this point in some 
general direction inside AH n . If h\j is of x-type 2, then /3j* (x) = 0. Perturbing 
x to x implies in general that |/3|y (x)| > 0. In the geometric interpretation, this 
has the effect of an explosion. If hfj is of x-type 3, then < (x)| < 00 and 
shifting x in general direction does not change this. In our geometric interpretation 
this means that the transformed leg, before and after perturbation, has finite but 
positive length. This explains the name 'hook on scale'. The name 'top angle' refers 
to the place of the x-tag 'top' in the hooked tree ht(x). 

Definition 7.36 

Let x € XAH n be in standard form. A point q € Dom„ (x) is valid with respect to 
x iff the following two conditions on the coordinates of q hold: 



non-orthogonality: no explosion angle is orthogonal to the corresponding angle 
otL (x). The top angle ct\j 2 {q) is not orthogonal to the top angle a\j 2 {x). 

finiteness: \(3^ k (q)\ < 00 for all ratio coordinates 0% {q). 
Remark 7.37 

tv x (x) is valid with respect to x.. 



Proof. The non-orthogonality is obvious, while the finiteness follows from Lemma 

TM. u 
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Figure 7.11: The poset induced by the predecessor relation. 



7.5.2 Example: ordering the sites for the drawing process. 
Example 7.38 



Figure 7.11 depicts the poset induced by the predecessor relation in our running 



example, see also Examples 7.13 and 7.32. This poset tells us, for example, that in 



order to draw pe in the top screen, we should draw pi, p$, and P2 first and in that 
order and, moreover, that we can ignore p4 and p§. So, the hooked path L x (pq) is 
given by L x (p 6 ) = (p 6 ,P2,P3,Pi)- 

This hooked path L x (po) tells us how to draw p$ in the top x-screen using ht(x). 
First pi is placed in the origin. Then p% is drawn using l x {ps) = CK13. As by now 
the leg P1P3 exists, we can apply 1 X {P2) = ^13 in order to add the leg piP2- Finally, 
the leg pipe is formed by applying l x {p 6 ) = h\\ . 



7.5.3 The construction draw x (g). 



Construction [7.42 indicates how to fill x-screens with sites p\ , . . . , p n , given a valid 
point q € Dom„(x). First we define an ordering on the x-screens. 

Definition 7.39 

An x-screen T is above an x-screen S iff the set of labels of S is a subset of the 
set of labels of T . 



The following lemmas relate the x-type of a hook to clusters C in C(x). 

Lemma 7.40 

Let h\j be an x-tag. 

(i) If x-type( hfj ) = 3 then C tj = C ijk = C lk . 

(ii) If x-type(h\j ) = 2 then Cij 7^ Cik and moreover, Cij is directly above Cik- 

Proof. We prove the two claims in the Lemma. 

(i) By construction, an x-type 3 hook involves the minimal labels of the first, 
the second, and another maximal cluster in a fixed cluster. Therefore 



Cij Cijk — C^k . 



(ii) Analyzing the cases in Definition 7.25 shows that i points to the minimal 
label in some cluster S and k points to the second lowest label in S. Moreover, i is 
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a minimal label for one of the maximal clusters of T, the cluster directly above S, 
while j is the minimal label of one of the other maximal clusters in T. Therefore 

Cjj = T 7^ S = Cik- M 

Lemma 7.41 

Let C be a x-cluster. Let p mi and p m2 be the sites with minimal labels in the first, 
resp., the second maximal cluster ofC. Let pk € C, but p^ ^ p mil p m2J and let 
hfj = l x (Pk)- Thenpi,pj € C. 



Proof. Let M ly M 2 , . . . , M s be the maximal subclusters of C. Suppos e tha t p^ £ 
Mi for 1 in {1, ...,s}. Suppose that x-type(/i^' ) = 3. By Definition 7.25 , there 
are two possibilities. Either pk has minimal label in Mi, and in this case I > 2, 
Pi = Pmn and pj = p rri2 or, pkiPi and pj belong all th ree to Mi. We are left with 
the case x-type(hfj ) = 2. According to Definition 7.25 , pk has minimal label in the 
second maximal cluster of some cluster D. As pk ^ p m! , this means that either D 
equals Mi or D is a subcluster of Mi. It follows from Lemma 7.4C that in both cases 
Pi,Pj S C, since these two sites are contained in the cluster directly above D. ■ 



Construction 7.42 ( draw x (q) ). 

Filling a;-screens. 

Input: x-screens, valid q £ Dom„(i). 

Output: x-screens filled with q, prescription draw x (q). 

Step 1. Choosing representatives. 

Choose representatives for the top angle aij 2 (q) and the explosion angles a*{, (q) 
in ht(q). Do so by changing the present value, if necessary, to the value modulo 
2ir most close to the corresponding value aij 2 (x) or a^ k (x) in tv(x ). T his choice 
is unique because of the non-orthogonality condition in Definition 7.36| . See also 



Figure 7.12, which demonstrates that there always exists a value of a(q) within 



distance 7r/2 of a(x) 




Figure 7.12: Non-orthogonality condition. 



Simultaneously, change the sign of an explosion ratio (3fj (q) whenever you change 
the explosion angle by a 1 ^ (q) = af- (q) + ir. Otherwise the Kf~ Klein bottle iden- 



tification 7.1 is not obliged. 

Step 2. Screen orientation. 

We define the screen orientation of an x-screen S = Sc by induction on the 



'above' relation introduced in Definition 7.39| . Assume that the screen orientation 



is defined for all x-screens above S. Let T be the screen directly above S. Following 
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the labeling in Definition 7.25, the screen orientation Os is defined by 



Otij 2 , 






S is a top screen, 
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(7.3) 



Step 3. Drawing the sites. 

The recipe for drawing a site Pi in a screen S — Sc for i G C, is given inductively, by 
induction on the predecessor relation. Below, T denotes the screen directly above 
5. 

l: Site pi is the first site in every path L x (pi) and is drawn in 'the origin' in any 

screen where pi occurs. 
2: Next suppose that i ^ 1 and that all sites p m above pi with m € C have already 

been drawn. 

3: if pi is the first site in the first cluster of C then 
4: pi is drawn in the origin of S. 

5: else if pi is the first site in the second cluster of C then 
6: put pi at (cos Os, sin Os) 
7: else 

8: l :c {Pi) = tijki wnere Pj and pk are above pi. By the induction hypothesis, 
the leg PjPk is already present in S, see also Lemma 7.41 . The leg PjPi is 
constructed as follows: 
9: Rotate PjPk over o?^ k around pj. 
10: Scale the result by a factor /3j l k with respect to pj . 
11: To summarize, 



Pi = Pfk Tot a % (Pk ~P])+ Pj ■ 



(7.4) 



12: By induction, the length of the leg PjPk is finite. The finiteness condition in 
Definition 7.36 assures that the resulting leg PjPi is finite as well. 



Remark 7.43 

Note that the legs piPj constructed in draw x (q) while filling the screens are exactly 
the IE-legs. 

Some drawings are better than others. 
Definition 7.44 

A valid q € Dom„(i) is accepted with respect to x iff in any x-screen S filled by 
draw x (g), q-sites that belong to distinct maximal x-clusters of S do not coincide. 



Remark 7.45 

Suppose that q e Dom„(i) is valid and accepted. If type x {(3f- ) = 3, then 
Definition 7.46 

Every valid and accepted q € Dom„(i) gives a filing of the x-screens consisting of 
71 sites. The set of all filled x-screens that can be obtained in this way is denoted by 
Scr n (x). 



7.6. Reading off data elements from Riled screens. 



Ill 
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■6 2 \ 



6 
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5* 



Figure 7.13: re-screens filled with some q £ Dom„(i). 



7.5.4 Example: filling x-screens. 
Example 7.47 



In F i gure 7.15 we h ave filled the ir-screens of our running example, see Exam ples 
7.1S| , 7.32] and 7.3£ , with the q- values present in the the fifth column in Table |7.2| . 



The chosen representatives for the q- values are given in the last column. According 



to Step 2 in Construction 7.42, we have changed the representatives for /1J3 and 



/i|f . Let us look more carefully at the relative positions of the sites p\, p%, and p§ 
in the top screen SVgi and the screen S\ t 2,6- It seems like the complete cluster 1, 2, 6 
has been rotated over n w.r.t. site p\ in the bottom screen. In order to find an 
explanation, we compute the coordinates of pi in both and 51.2,6- The hooked 
path L x (p2) equals (p2,pz,pi). In both screens p\ — (0,0). In the top screen, 
Ps = (cos Os [6] , sin Os [6] ) , wher e O s [R] = "13- We determine the position of P2 in 
the top screen, using Equation |7.4| : 

P2(S[e]) = Pil rot^i (p 3 ) = fill (cos a 13 +a\l, sin a 13 + a\l). 

In 51.2,6, however, pi is the first site in the second cluster, so 

P2(*Sl,2,6) = ( cosC, Si,2,6> sinC 's'i,2,6)- 

According to Step 2 in Construction 7.42] , 



Csi.2,6 = °s l6] + a\l = "13 + all > 

as i\ = 1, «2 = 3 and ji — 1, ja = 2. We conclude that the only difference between 
the argument of p 3 in S'jg] and 5*1,2,6 can be caused by /?i| : this ratio is used in 
drawing p2 in S[q] but not used in drawing p2 in Si, 2, 6- Indeed, sign(/3|| ) = — 1, 
which causes a n shift in arg(p2) from the top screen to the screen £1,2,6- 



7.6 Reading off data elements from filled screens. 

Let q 6 Dom„(i). Suppose we have filled the cc-screens with g-sites by applying 
dr&w x (q). In this section we tell exactly how to read off angles and hooks coor- 
dinates from the q-filled cc-screens. More precisely, we explain how to determine 
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every angle ixy between two sites pi (q) and pj (q) and every hook h A between three 
sites Pi(q), Pj(q), and Pk{q)- Basically we read off lacking coordinates in a suitable 
x-screen. Recall the notions of separating clusters and separating screens that were 



introduced in Definition 7.24. These notions are used for fixing the ir-screens used 



for determining the angles and hooks coordinates. 
Remark 7.48 

The separating screen Sij(x) is the only cc-screen in which it is guaranteed that 
Pi(q) 7^ Pj{q) for q being valid and accepted. 



7.6.1 Determining hooks and angles from draw^g). 
Definition 7.49 

Let q £ Dom„(:r) be valid and accepted. From drecw x (q) we determine a point 
q = readj; (q) £ AH n as follows. 

(i) oiij{q) £ R/27rZ is the angle of the directed line through pt(q) and Pj(q) in 
the q-filled screen Sij(x). 
( H) a t (?) = a ik (q) ~ aij (q) . 

(Hi) Read off all sites coordinates in screen Sijk- If Pi(q) ^ Pj(q) in Sijk, then 
where 

Else, P$(q)=l/fti(q). 

Note that by definition of separating screen, not all three sites pi, pj, andpfc coincide 
in Sijk- Therefore pi (q) — Pj(q) in Sijk implies that pi{q) ^Pk(q)- 

pk 



(Pk(q) -Pi(q)) ■ (Pj(q) - Pi{q)) 

\p'M-pM\ 2 



(7.5) 



rot ai* (<i)(pj(q) -Pi(q)) +Pi(q)- 



(7.6) 




Figure 7.14: Determining . 



Remark 7.50 

We give a geometric explanation for Definition 7.49. It is rather straightforward to 
read off Oiij{q) and a\j (q), so let us focus on 0r? (q). Go to the separating screen 
Sijk ( x ) ■ Write down the coordinates of pi (q) , pj (q) , and pk (q) in this screen. First 
consider the case where Pi{q) ^ Pk{q)- If proceeding through the definition in the 



right order, we know by now a\j (q). See also Figure 7.14. Use o& (q) in order to 



rotate the leg pj — pi onto the line lik through pi and pk, thereby fixing a point 
p'j. Finally determine the coordinate of p^ on the ij'-axis with respect to the 'unit 
vector' pj — pi as follows. Let x = pk — Pi, y = p'j — Pi, and 9 the angle between 
x and y. From the definition of dot product x ■ y = \x\\y \ cos#, it follows that the 
quantity ||| cos 9 we look for is given by (x ■ y)/\y\ 2 . Note that in our case 9 always 
equals either or n. 



7.6. Reading off data elements from Riled screens. 
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Proposition 7.51 

The components of read^. (q) are infinitely often differentiable (C°° ) on the set of 
valid and accepted q S Dom„(a:). 



Proof. The main problem is the continuity of Ye&d x (q), so this is stressed in the 
proof. To make the proof complete, replace every occurrence of 'continuous' by 
C°°. 

(i) AH n is smooth as a direct product of circles and Klein bottles. Dom„(i) 
is a factor of AH n and is a product of circles and Klein bottles as well. read x (g) 
determines a point in AH n . 

(ii) We restrict ourselves to valid q £ Dom„(x). This means that we have 
imposed two conditions on q: the non-orthogonality condition, which is is an open 
condition as it assures that one point of a circle S 1 is avoided. The finiteness 
condition is open as well, as one point, infinity, of P 1 is avoided. 

(iii) For valid q <E Dom n (:c), the construction draw I ((j) is defined. The construc- 
tion consists of three steps. First of all representatives are chosen for the top angle 
and explosion angles. The chosen representatives depend continuously on q. No 
jumps can occur due to the non-orthogonality condition. In step 2 screen orienta- 
tions Os are defined. The formulas presented in Equation 7.3, depend continuously 
on the chosen representatives for q. In the third step, points are drawn in the screen. 
The positions of these points depend cont inuously on the screen orientation Os and 
the q coordinates, compare Equation 7.4 . 

(iv) Filled screens that are not accepted are thrown away. A filled screen is not 
accepted when certain clustering between sites occurs. Avoiding clustering is yet 
another open condition. 

(v) For valid and accepted q E Dom„ (i), the point read-^g) £ AH n is de fined. 
The formulas for the coordinates of read^g) are given in Definition ftA^ . All 
equations depend continuously on the position of the sites in accepted draw^ (q) . ■ 



7.6.2 Example: reading off hooked tree elements. 



Example 7.52 

Applying D efinit ion 7.49, we read off the values of the angle and hooks in Dom„(x) 
from Figure 7.13 These values are presented in Table 7J3. Next to the data elements 
to read off, the read off screen is given. Compare the values in Table 7J? with the 
values of q € Dom„(i) in Table 7J1 used above to fill the x-screens. Note that 



these two sets of values are indeed the same but only up to identification in the 
appropriate Klein bottle. 



angle screen value angle value ratio screen value 



«1,3 


S[ 6 ] 


69° 












OL\,2 


5*1,2,6 


157° 




88° 


Pli 


5 [6] 


-.35 


Q!l,4 


«%] 


-107° 




-176° 


Pli 


S[ 6 ] 


1.92 


«3,5 


5*3,5 


-123° 


a 35 
"31 


-12° 


Pi! 


5[ 6 ] 


.34 


Oilfi 


5*1,2,6 


146° 


rylS 
"12 


-11° 


Pit 


5*1,2,6 


.91 



Table 7.3: Read off values. 
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7.7 Consistency Theorem. 

In this section we show that our construction of filling screens by means of hooked 
tree data is consistent with the procedure of reading off data from the filled screens. 



7.7.1 Managing x-type 3 hooks. 
Lemma 7.53 

Let q € Dom„(i) be valid and accepted and put q = read x (g) £ AH n . Suppose that 
type x (h%) = 3. Then 



4 (9) = 





Figure 7.15: (q) > and (q) < 0. 



Proof. Follow the labeling in Definition 7.25: i = ji, j = j'2, and k — j t . For 
convenience we set j% = 1, j2 = 2 and j'j = t. By Lemma 7.40, all evaluations take 
place in x-screen S = S12 = 5i2t = Sit- For an illustration of the proof, see Figure 
7.15. As pi — (0,0) and P2 = (cos Os, sinOg) it follows that 



ait (9) - "12(g) = arg(p t ) - arg(p 2 ) = arg(p t ) - O s . 



If /3» (q) > 0, then arg(p t ) = O s + a\\ (q), so a» (g) = a}* (<?). 
hand, #| (<?) < 0, then arg(p t ) - ^„ _l. ™it f„\ _l ^ o„ „,i* 
required. 



+ all (q) + 7T. So, aj| (g) 



If, on the other 
= "12 (?) + 7T as 



Lemma 7.54 

Let q G Dom„(i) be valid and accepted. Set q 
type x (f3%) = 3. Then 

ft) \^m = \^{i)\- 



= read x (g) G AH n . Suppose that 



(ii) 



ffi (9) 



> 0. 



Proof. As in the proof of Lemma 7.53 set i = 1, j = 2, k = t and conclude from 
Lemma 7.4C that the read off screen is given by screen S = S12 = Sit- We use 



Equation 75 for determining (<?), while remembering that p\ = (0,0). 



pt -p 2 
p'2 ■ p'2 



(7.7) 



p' 2 is given by 

P2 



r ot a itp 2 = (cos[O s +a{l(q)],sin[O s 



0-1 



(?)])■ 



7.7. Consistency Theorem. 
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Therefore \p' 2 \ 



so j3\ 2 (q) reduces to (3) 



Pt • p' 2 . Observe that p 2 = 



(cos Os, sin Os) and apply Lemma |7.53 . If (3\ 2 (q) > 0, then a\ 2 (q) — a{ 2 (g) and 



by evaluating Equation 7.7 it follows that (3\ 2 (q) = (3\ 2 (q). If, on the other hand, 



(q) < 0, then a{ 2 (q) = a\ 2 (q) + tt, so p t = —fi\ 2 {q)p' 2 , which implies that 
0u (?) — — fii2 (?)> as claimed. ■ 



Example 7.55 

In the hooked tree in our running example, see Examples 7.13| , 7.32) , 7.38, 7.47 and 
two x-type 3 hooks occur, and /ijj ■ The q- values, that are those values 
that where used to fill the cc-screens, and the read x (g)-v alue s, the values that where 
read off from the filled x-screens, are collected in Table 7.4 . 



q coordinate reaci^g) coordinate 



( - 1.92,4° ) 
h\l (.91,-11°) 



"13 
,16 



( 1.92,-176° 
( .91,-11° 



Table 7.4: Swap of representatives for type^ 3 hooks. 



Note that j3\^ (q) < 0. According to Lemmas 7.53 and 7.54 this implies that the 



representative in swaps. Indeed, /3|| 



changes sign, while a]^ (q) 



'13 



Ratio j3\ 2 (q) is positive so here nothing should change and this is reflected in the 
values in the table. 



7.7.2 Managing x-type 2 hooks. 

Next lemma shows that in case of a x-type 2 hook, the representative in the Klein 
bottle stays the same. 



Lemma 7.56 

Let q e Dom„(i) be valid and accepted. Put q 
x-type(h$ ) = 2. Then 



read T (<7) G AH n . Suppose that 



(q),ctf (q)) 



(fl? («), <*?< (<?)) 



Proof. Throughout the proof, we label as in Definition 7.25 . The screen orienta- 
tions 0s for the different cases are defined in Equation 7.3. Set k — j 2 and i = j±. 
We distinguish three possibilities for the hinge point pj 1 . 

(i) ji — i\. In this case, j — i 2 . First we show that a£ 3 £ (q ) = oP 3 £ (q). In 
order to determine oa^iq), we need the coordinates Pi 1 {T) = (0,0) and Pi 2 (T) — 
(cos0T,sin0T) in T = Si 1 i 2 . It follows that a^i^q) = arg(p i2 ) = Or- For 
determining j 2 (q) = Oij t j 2 (q) we compute in screen S the coordinates of pj 1 (S) = 
(0,0) and: 



Phi s ) = (cos 0s, sin c 



(cos 0-z 



a^!,sin0 T 



M2 ■ 



We get that c££ (q) 



T + alf(q) 



On 



a- 



(?)■ 



Next we determine /3^ 2 2 (q). For this purpose we compute the coordinates of p^, 
Pi 2 , p' i2 , and pj 2 in Si 1 i 2 j 2 — T. We gave the first two above, the latter two are given 
by 



(cos Ot + o 1 ^, sin O t 
B lin p' 

' l\l 2 ri 2 
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By substituting those coordinates in Equation 7.5 we obtain that Plli 2 (q ) 
PZil (?)' as claimed. 



J 2 




Figure 7.16: /3g£ (q) = (q). In the picture, (q) < 0. 



(ii) j\ — %2- In this case, j = i\. Again we start by showing that a^if io) = 
a) 2 J 2 (q). For an illustration of what happens in screen T — Si 1 i 2 , consult Figure 



?. f ?. 1 \ J- / -L A ~ ____ „ L „. z , 

7.l4 As in the case ]\ = i\ it holds that on x i 2 (q) = Ot- Therefore, a^^g) 
Ot + 7T- Angle aj 3j - 2 (g) = <x,- U2 (g) is determined from pj x (S) = (0, 0) and pj 2 (S) 
(cosOs,smOs) in S = 5j i:)2 = 5: screen orientation O5 is given by O5 = Ot 



7r. Therefore a.^ J2 (g) = arg(pj 



q t + a Jin 



7T. SO 



On 



TT-Ol 



1132 (q). 



For det ermining /3| 2 ^ (g), which occurs in the x-screen T, we apply Equations 7.5 
and 7.6 directly: 



(Pj2 • (gil -Pi 2 ) 

(Pj2 -P12) ■ ip'n -Pi 3 ), 

^ (9). 



(iii) ji = i f , t = 3, . . . m. In this case j = i^. We start by analyzing ot? ^ . 
For determining a itil (q) we compute in T — coordinates p^iT) — (0,0), 

Pi 2 {T) = (cos Ot, sin Ot), and 



(g)(cos(0 T + a££ (g)), S in(0 T + a*% («))). 



Therefore, 



"Mi® 



Or + a^(g), > 0, 

0T + <i<(g)+7r, < 0. 



We use the coordinates, Pi t (S) = pj 1 (S) = (0,0) and Pj,jS) = (cos 0s, sinOg) in 



S = S„ 



H32 



Sj 1 j 2 for determining o>i t j 2 (q). From Equation |7.3| it follows that 



O s = 



Ot + <2 + (g) > 0, 

OT + «:;:;(g) + a^(g), /^(g) < 0. 



7.8. XAH n is smooth. 
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After subtracting it follows that aj*^ (q) = a£*j° (q). We only have to show that 
(#) = ^i'ii (^)- We ^° so ^ filling in coordinates in Equation L5, The relevant 



coordinates in screen Si. 

Ph 

Pit 



Pn 



Therefore 



T are as follows. 
= (0,0), 

= (cos0t,sui0t), 

= (?>ot a *xH (pi 2 -Pi^+Ph, 

= ft* (g)rot a<lit (p i2 ), 
= rot i tj2 (p i± -p it ) +p it , 
= rot (-pij +p it , 
= /tlfrot^^-PiJ- 
■ ' {-Pi*) + Pi 



Pi t , 



H32 



<T (9)rot Q , t32 (-p lt ) • rot aH32 (-p h ) 
i t i i Hn 

rot a ,„ 2 (-p it ) • rot a ,„2 (~Pi t ) 



(ff). 



Example 7.57 

In our running example we have two type^ hooks, and . The reader is 
invited to check that the representatives computed for read x (q) in and K^f 
that are the same as those for q itself. See in particular Example [7.47 and Example 
7J5|. 



7.7.3 Consistency theorem for hooked tree elements. 
Theorem 7.58 

Fix x £ XAH n . Let q £ Dom„(i) be valid and accepted. Then 



Dom„(x) 



is the identity map. 



Scr„(a;) 
drecw x (q) 



ALL 



reads (?) 



Dom„(i), 



Pro of. We analyze the three types of data elements pres ent in ht(x), see Definition 
[7.25| . The type^ 3 hooks are dealt with in Lemma 7.53 and the type x 2 hooks in 
Lemma 7.56) . We are left with the type^ 2. a angle between the top screen and its 
second maximal cluster. In Definition 7.25 this angle is labeled a\j 2 . The separating 
screen of p\ and pj 2 is the top screen S\ n ] , so we have to read off in this screen. The 
coordinates ofpi and pj 2 in S[ n ] are pi = (0,0) and pj 2 = (cosaij 2 (q),smaij 2 (q)). 
The read off value a\j 2 (q) equals arg(pj 2 ) = a\j 2 (q). This completes the proof of 
the theorem. ■ 



7.8 XAH n is smooth. 

Lemma 7.59 

Fix x £ XAH n . Let q £ XAH n be such that: 
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ft) tv x {q) is valid and accepted; 
(ii) q = V'AffJM) f or some i c q] e conf„ 
Then q = re&d x (tv x (q)). 



Proof. The consistency theorem, Theorem 7.58 , proves the claim for those q coor- 
dinates that belong to Dom„ (x) . In the rest of the proof, we analyze mod tt 
and h% . Write q := re&d x (tv x (q)). 



First we show that a a b(q) — oi a b(q) for a, b <E 1, . . . , n, with a ^= b. Definition 7.49 
tells that we have to determine a a b(q) in the screen S a b- Let j\ and ji be the first 
label in the first cluster and the first label in the second cluster of screen S a b- Then 



Pj 1 and pj 2 are the sites drawn first and second in Construction 7.42. If a =/= Ji, J2, 
then p a was constructed out of the x-leg Pj 1 j 2 by repeatedly applying hooks from 
ht(x) in the order indicated by the hooked path L x (p a ). The same holds for pb- For 
hooks h\j that are referred to in ht(x) one has hf- (q) = h 1 ^ (q). As a consequence, 
the positions of p a and pb relative to pj 1 and pj 2 are the same in S a b as in any 
configuration c q € [c q ]. 

The angle of the directed line passing first through pj 1 and then through pj 2 in 
S a b is exactly the screen orientation Os ab = Osu ^ ■ Writing out any screen ori- 



entation according to the recursion in Equation 7.3 gives a list of data elements 



corresponding exactly to the hooked path L x (p D ). As all angles occurring in the 



screen orientation definition are coordinates of Dom„(i), Theorem 7.58 applies 
again. At any type 2 angle afj however a tt switch can occur: the corresponding 
type 2 ratio /3$ is not used in the definition. Moreover, representatives might swap 



in Step 1 of Construction 7.42. It follows that 



O s ■ 

°J\32 



a h32 (<?) + 



To show that hf~ (q) — hfj (q) is easy now. The proof consists of three observations: 
(i) Any (/-filled x-screen 5* is just a copy of the original configuration, due to 
the Consistency Theorem. But, it also satisfies the following: points not belonging 
to cluster C(S) corresponding to screen S are omitted; points in S are scaled such 
that the distance from pj t t o pj 2 equals one; points in S are possibly rotated over 
tt according to Const ructi on |7.42| , Step 1 . 



(ii) By Definition 7.49 , a 1 ^ (q) — ctik(q) — onj{q). From the proof of part 1. of 

(q) equals a\!j 



this lemma it follows that the angle ofi (q) equals af~ (q) up to tt 



(iii) In case ctf- (q) = a-j" (q) + tt this 7r-switch is automatically corrected by the 



sign of /3lj , as afj (q) is used in Equation 7.6 



Lemma 7.60 

Fix x S XAH n . Let U C AH n be the set of points in AH n such that tv x (u) is valid 
and accepted for all u £ U . Then 

XAH n r\U = graph(ctv x (rend x )). (7.8) 

Proof. We prove the two inclusions. First let y S XAH n n U. Then its projection 
q y := tv x {y) G Dom„(i) is valid and accepted with respect to x. There always exist 
points z e XAH n r\U of the form z — ip 4 ([c z ]), for [c 2 ] € conf„, that are arbitrary 



close to y. For such points z expression d7.8| ) holds by Lemma 7.59. So these points 
can be expressed as the image of a configuration consisting of non-coinciding points. 
Therefore these points z belong to the graph of graph(cto a: (read a ;)). This implies 
that y itself belongs to the graph as well, as y can be constructed as limit point of 
points of the form z = 4>AH (I c ^])' f° r con Li in XAH n . 



7.9. Voronoi diagrams in the x-screens. 
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Secondly, suppose we start with a point g q on the graph. This implies that q = 
tv x (g q ) £ Dom„(x) is valid and accepted. As both conditions (valid, accepted) are 
open conditions, there exists e with < e << 1, such that if we replace all g-ratios 
that vanish by e that then q(e) is still valid and accepted. But this means that 
dracw x (q(e)) produces a filling of the x-screens in which all sites in the top screen 
are distinct. Writing down the coordinates gives a class [c g ( e )] £ conf„. Then 

^AflJ[ c *)l) = {q(e);ctv x (read x (q(e)))). 
The claim follows by observing that XAH n is the closure of tpAH ( coni "n) m AH n . 



Theorem 7.61 

XAH n is a submanifold of AH n 



Proof. Being a submanifold of a manifold can be determined locally, see [BG] 



section 6.2. Lemma 7.60 shows that XAH n can locally be written as the graph 



of a function, and this is one of the characterizations of submanifold. See [BG] 
Theorem 2.1.2. notably part (iv). I 



7.9 Voronoi diagrams in the x-screens. 

So far in this chapter we have introduced the space XAH n in order to model config- 
urations of coinciding sites in the plane. We have seen that we can associate with 
every x £ XAH n a family of x-screens. Moreover, for almost all q £ XAH n , we can 
fill the x-screens with q-sites so that for any two q sites Pi(q) and Pj(q) there exists 
an x-screen where Pi(q) and Pj(q) are distinct. 

In this section we add Voronoi diagrams to the x-screens. As in any x-screen at least 
two sites are distinct, we could apply the ordinary definition of Voronoi diagram. 
But there is an ambiguity in the screen model, as we take angles between pairs 
of points in R/7rZ and allow negative ratios 0f k . 

Example 7.62 

Recall Example 7.47. This example demonstrates that q-clusters can get rotated 



over 7r in distinct x-screens. 



7.9.1 Connection with Kontsevich-Soibelman space. 

We nail down the clusters in such a way that they cannot get swapped over 7r in 
distinct x-screens. We do this by taking angles between pairs of points in M/27rZ 
and allowing ratios (3^ £ [0,oo] only. That is, we consider the map 

^EDAH n : conf » -> (K/2ttZ)(S) x ([0,oo] x R/2ttZ) 6 ('), 
[(p u ..., Pn )] I ► ((Oii)l<i<4<n, ,o# ))• 

where i, j, and k are pairwise distinct indices. Denote by EDAH n the product 
(R/2ttZ)(S) x ([0,oo] x R/2ttZ) 6 ^) . Recall the manifold with corners FM^n) in- 
troduced by Kontsevich and Soibelman, see K| and described in Section 7.1.2. 



Proposition 7.63 

Let XEDAH n denote the closure of^EDAH ( coni n) in EDAH n . Then 

X EDAH n = FM 2 (n). 
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Proof. FM2 is defined as the closure of the image of conf„ under the map 

[{pi,...,Pm}] >-> {{ a ij)l<i<j<n,0ij), 



where i, j, and k are pairwise distinct indices and f3jj 
the claim follows directly. 



/■'< 


—Ph\ 


\pl 


-Pj\ 



As 1 





3 


1 • 




• 

2 



Figure 7.17: The Voronoi diagram of pi(t), P2(t) and p^if). On the left for negative 
t, on the right for positive t. 



Example 7.64 

Consider Pl {t) = £(1,0), p 2 (*) = t(-|, and p 3 (t) = *(— |, — |V3). We com- 

pare the angles associated to pi, P2, and P3 in XAH[3] and FM2(3) for positive and 
negative t. First we consider XAff[3]. The undirected angles 012 = — J, 013 = f 
and 023 = § do not change at t = 0. In FM2(3), however, angles ai2, CK13 and a 2 3 
change by it at £ = 0. Think of this as moving into another corner in the manifold 
with corners. In Figure 7.17 the Voronoi diagram of pi(t), P2(t) and Pa(t) is shown. 
On the left for any negative value of t, on the right for any positive value of t. Note 
the jump in the diagram at t = 0. It demonstrates the need of the directed angles 
in the screen model if we want to use it for displaying Voronoi diagrams. 



a = ao 



a = «o + 7r 





FM 2 {n) 



If 



a — a?o 




P = 



XAH n 



(3 = 



Figure 7.18: The manifold with corners FM2(n) above the smooth manifold XAH n , 
with a, (3, and (3' in Dom„(i). That slice is presented where all coordinates are 
fixed except for (3 and (3'. 



Theorem 7.65 



7.9. Voronoi diagrams in the x-screens. 
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See also Figure 7.1l\ . Consider the 



map 



I ■ FM 2 (n) -> XAH n , 

(K-)i<i<i<„,(/3t,^)) - (m 1 < l<J < n ),m^ai';)}), 

where [.] denotes the class of [((if* , a\j )] in the Klein bottle Kff with respect to 
the equivalence relation Fix x E XAH n . Let m be the number of (iff (x) in 

Dom„(i) such that (if* (x) = 0. Then 



\r\x)\ 



im+1 



Proof. Fix a class [(/3J* (x),afj (x))] such that (3ff (x) = 0. Any such class has two 



elements, (0, af- ) and (0, afi- +%). By Lemma 7.60| , XAH n can be written as the 



graph of a function on Dom„(i). Therefore those /3y components that are not 
in Dom„(x) are a function of those components that are in Dom„(i). Moreover, 
read x (x) even determines the representative for those /3y components that are 
not in Dom„(x). For any [3$ S Dom„(i), the two representatives correspond to 
two distinct pieces of the boundary in FM^ri). We conclude that we can use the 
signvector of those (3ff (x) € Dom„ (x) with (iff (x) — as a coordinate system 
pointing to the distinct corners in the fiber of the map /. The additional 1 in 2 m+1 
comes from the two representatives of the unique -component in Dom„(i). ■ 

Example 7.66 

Consider x € XAH3 given by 

S12 ai 3 a 23 (/3i2,ail) (021 ^l) (Pii,<4i) 
f I f (i,o) (0,-ff) (0,§) 



Theorem 7.21 associates to x the nest ({1, 2, 3}, {2, 3}). The factor Dorr^x) con- 
sists of the two components c?i2(x) and h^i (x). The standard form of x restricted 
to Dom3(x) is given by the representatives 012(2;) = ^ and a^f (x) = — ff, com- 



pare Definition |7.33| . The fiber / _1 (x) consists exactly of those points in FM 2 (3) 
that are mapped on the same standard form in XAH n . We obtain these points by 
analyzing the four distinct representatives for cti2 and o^f . Choose for example 
— ana a 21 — 



Q'12 = — ^ and a|f = — ff- It follows from the nest structure that 0113 = CC12, 
while a 23 = al? + Q!2i = -¥ + ? = -f- The four P oints A, ■ ■ ■ , A in FM 2 (3) 



obtained in this way with coordinates (0:12, 0:13, 023, /3if , /?23 , /?32 )j are given by 

fx = (f,f,f ,1,0,0), 
/2 = (f , f , — f , 1) 0,0), 
h = (-f>-f ,-f,l,0,0), 
/4 = ,1,0,0). 

Remark 7.67 

In Chapter |5| we have encountered a similar situation in comparing the spaces CUA3 
and C-DA3 of undirected and directed angles on three points. Recall particularly 



Examples 5.30 and 5.31 



We can define and fill x-sreens in FM2(n) as well. Fix a point x € FM2(n). The 
set of vanishing ratios (3\j (x) for x defines a family of x-clusters and corresponding 



x-screens by applying the methods of Section 7.4. A stripped version of Construe 
tion 



7.42 fills the x-screens, given some q € Dom„(x). In short, we discuss the 
adaptations in Construction 7.42] . 



Ad Step 1. We do not have representatives in EDAH n , so we skip Step 1. 
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Ad Step 2. The definition of Os still holds, but negative ratios 01* do not occur 
in EDAH n . 

Ad Step 3. In Step 3 nothing changes. 



7.9.2 Adding Voronoi diagrams in the x-screens. 

In this section we define a Voronoi diagram Vfm (x) for x € FM2 (n) . 
Definition 7.68 

Let x £ FM2(n). The Voronoi diagram Vfm{x) is defined in terms of the x-screens 
filled by draw x {x). Fix an x-screen T filled by draw x {x). Let C be the cluster 
corresponding to T. Its maximal subclusters are denoted by C\,... ,C m . Suppose 
that the coordinates of C\, . . . , C m in T are given by t%, . . . , t m . 

Initialization. V(T) is defined as the Voronoi diagram V({ti, . . . ,t m }). 

Completion. This is a recursive step from the screens displaying single sites up 
to the top screen. Assume that the maximal subscreens Si, ... , S m of T have 
been completed. Denote the completed diagram in Si by Vfm (Si). After the 
initialization step, any cluster Ci lives in a Voronoi cell V(U) in T. Define 

V FM (U) = V(ti)nV FM (Si), 

and 

Vfm(T) = |J V FM (U). 

i—l...m 

The Voronoi diagram Vf m (x) is the rooted tree of all Vf m (T) , where T runs over all 
x-screens. The ordering on Vfm(T) is inherited from the ordering on the x-screens. 






11 





11 




Figure 7.19: Constructing Vfm(x)'. on the left, initialization; on the right, comple- 
tion. 



Example 7.69 

Recall Example 4.42] . In this example we consider an ordered set of sites in R[£] x M.[t] 
given by 

S(t) = ((-6,4),(-2,-6),(-l,-2-3t),(0,-3t),(3,5),(6,-3), 

(-1 - 3i, -2 + t), (-2t, -2t), (-2t, 2t), (2t, 0), [2t, 2t), (-1 + 2t, -2)). 



7.10. Conclusion. 



123 



For any pair of sites Pi(t) and Pj(t) we have defined an angle or direction aij(t) at 
t = 0. In a similar way we can define a ratio 



_ u„M*)-Pfc(*)l 



#*(()) = lim 



no |pj(t) — pj-(t)| ' 

This establishes a map \ from a set 5(f) of n sites in R[i] x R[t] to a point x := 
x(5(0)) G FM2. We construct the Voronoi diagram Vfm(x) by applying Definition 



■6J . The nest of x is given by ({3, 7, 12}, {4, 8, 9, 10, 11}). Figure |7.19| shows on 
the left the three filled screens together with the Voronoi diagrams added in the 
initialization step. On the right the diagram Vfm (%) is shown. It is constructed by 
adding the two diagrams corresponding to clusters 3, 7, 12 and 4, 8, 9, 10, 11 in the 
cells of those clusters in the top screen. 

Example 7.70 

We analyze an example where the maximal depth of a site equals 2. We start with 
the ordered set of sites in R[t] x R[t] given by 5(f) = {pi,P2,P3,P4}, where 

Pi = (0,0), 

p 2 = f (cos |,sin|), 

P3 = P2 +t 2 (cos|,sin|), 

Pi = P3 + i 3 (cosf,sin : §). 



As in Example 7.69, one can construct a point y — x(5(0)) G FM2. The nest of y 
is given by (234, 34). In Figure 7.20] the tree of filled screens for y is shown: on the 



left after the initialization step but before the completion step of Definition 7.68 
on the right after the completion step. In this example one can see clearly how 
the Voronoi diagram Vfm(T) of the top screen T is built up recursively out of the 
screens corresponding with the subclusters of {1, 2, 3, 4}. 



Remark 7.71 

Note the compatibility with the construction of limit Voronoi diagram for sites in 
R[t] x R[i] as constructed in Chapter 3. Compare also Examples 7.69 and 7.70. 



7.10 Conclusion. 

7.10.1 An easy model to remember. 

In this chapter we have presented a real version of the Fulton-MacPherson compact- 
ification X [n] . An advantage of our approach is that we do not need any machinery 
from algebraic geometry neither in constructing our compactification space XAH n 
nor in proving that it is a smooth manifold. Instead of algebraic blowups we have 
used the angles between two points and the hooks between three points. This idea 
was already proposed by Kontsevich and Soibelman. We have shown how to adapt 
their approach so that the resulting modified compactification space is a smooth 
manifold instead of just a manifold with corners. Moreover, we have shown that 
an explicit analysis of the ratios on triples of points gives a combinatorial counter- 
part of Fulton-MacPherson's description of degenerated configurations in terms of 
screens. 

Because of its nature, the construction and description of the space XAH n can serve 
as a illustrative example to the theory of configuration spaces that has attracted a 
lot of interest in recent years. The setup in terms of screens can serve as a general 
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Figure 7.20: Constructing Vfm(u)' on the left, initialization; on the right, comple- 
tion. 



framework for studying properties of point sets in the plane. Especially for point 
sets that contain almost coinciding points, a description by angles and hooks can be 
a more robust alternative for storing the relative positions of the points. Another 
advantage of working with hooks and angles is that translations and scalings are 
already eliminated. That is, one concentrates on point configuration up to affine 
transformations. 



[FM], Theorem 2 states that the Fulton-MacPherson compactification X[n] equals 
the closure of CONF n in the product of those B1a(A s ) for S C {l,...,n} of 
cardinality 2 and 3. Note that our model also supports this 'three is enough' motto: 
we only consider properties of pairs of points (the angles) and triples of points (the 
hooks). 



7.10.2 Relation to earlier chapters. 

A first application is given by studying limits of Voronoi diagrams: we have intro- 
duced a method to associate a collection of filled x-screens to a point x E XAH n . In 
the final sections it turned out that this method can be ported to the manifold with 
corners FM 2 (n), introduced by Kontsevich-Soibclman. Just because of its corners, 
the screen structure associated to a point x in FM2(n) is well suited for modeling 
a possibly degenerated Voronoi diagram for x. Recall the description of Voronoi 
diagrams of polynomial sites, presented in Chapter ^: we have by now two notions 
of degenerated Voronoi diagram that match both combinatorially and according to 
shape with the notion of classic Voronoi diagram. 

From a set S(t) of n distinct sites {pi(t), . . . ,p n (t)} in R[t] x K[t] we can construct 
a point s € XAH n at t = 0, compare Example [7.69 . A construction in the other 



direction can be made as well. Suppose we are given some x G XAH n . And suppose 
the cc-screens are filled by x-sites pi(x), . . . ,p n (x). Define the depth of a screen S 



7.10. Conclusion. 
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as: 



depth(S) = 1 + # screens above S. 
Let pj be the site directly above pi with respect to the predecessor relation defined 



in Definition 7.27. The polynomial site Pi(t) representing pi (x) is defined as: 



Pi(t) = (0,0), 

p t (t) = Pj {t) +t dc P th( - s ^ (coordinates^; (a:)) in Sy). 

Mapping a set S(t) of polynomial sites to a point s £ XAH n and back results in 
this way in a normal form for S(t). 

In Chapters || and |6| we have studied the compactification CDA n . This compactifi- 



cation is the closure of the graph of the directed angle map, cf. Definition 5.11. One 



of the drawbacks of considering only directed angles is that collinear configurations 



can not be reconstructed from just these angles, compare Section 5.2. This prob- 
lem is solved in the compactification XAH n by adding the hooks between triples of 
points. 



In Chapter |6[ Theorem 6.11 we prove essentially that two data sets 7 n ,?7„ £ CDA n 



that are Euclidean close have Voronoi diagrams V(j n ),V(S n ) that are Hausdorff 
close. As there is an obvious continuous map (forget all ratios) from FM2(n) to 
CDA n , this continuity result also holds for the top screen of the Voronoi diagram 
Vfm(x). More precisely: suppose we have some x £ FM.2{n). Fill the top screen 
S\ n ] by x. The coordinates of the sites pi(x), . . . ,p n (x) can be read off from S\ n ], 
while the directed angles between pairs of sites are part of x. These data together 
determine a point 7„ £ CDA n . This establishes a map / : FM2 — * CDA n . We have 
the following situation: 

, / x / . Voronoi map . . 

M 2 (fi) — > CDA n — > diagram space, 

As both / and the 'Voronoi map' are continuous, the composition is continuous as 
well. In this short analysis, we have left out the incorporation of the compactness 
condition on the Voronoi map. This is left as an exercise to the reader. 
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